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Abstract 

o 

^ Massive gravity has seen a resurgence of interest due to recent progress which 

has overcome its traditional problems, yielding an avenue for addressing important 
open questions such as the cosmological constant naturalness problem. The possibility 
of a massive graviton has been studied on and off for the past 70 years. During 
this time, curiosities such as the vDVZ discontinuity and the Boulware-Deser ghost 
^ were uncovered. We re-derive these results in a pedagogical manner, and develop the 

in 



Stiikelberg formalism to discuss them from the modern effective field theory viewpoint. 
We review recent progress of the last decade, including the dissolution of the vDVZ 
ly-^ discontinuity via the Vainshtein screening mechanism, the existence of a consistent 

effective field theory with a stable hierarchy between the graviton mass and the cutoff, 
and the existence of particular interactions which raise the maximal effective field 
^ theory cutoff and remove the ghosts. In addition, we review some peculiarities of 

X ...sive paviton, on curved ,p«e. novel theories in three dimensions, and examples 

^ of the emergence of a massive graviton from extra-dimensions and brane worlds. 
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1 Introduction 



Our goal is to explore what happens when one tries to give the graviton a mass. This is a 
modification of gravity, so we first discuss what gravity is and what it means to modify it. 

1.1 General relativity is massless spin 2 

General relativity (GR) [1] is by now widely accepted as the correct theory of gravity at low 
energies or large distances. The discovery of GR was in many ways ahead of its time. It was 
a leap of insight, from the equivalence principle and general coordinate invariance, to a fully 
non-linear theory governing the dynamics of spacetime itself. It provided a solution, one 
more elaborate than necessary, to the problem of reconciling the insights of special relativity 
with the non-relativistic action at a distance of newtonian gravity. 

Had it not been for Einstein's intuition and years of hard work, general relativity 
would likely have been discovered anyway, but its discovery may have had to wait several 
more decades, until developments in field theory in the 1940's and 50's primed the culture. 
But in this hypothetical world without Einstein, the path of discovery would likely have been 
very different, and in many ways more logical. 

This logical path starts with the approach to field theory espoused in the first volume 
of Weinberg's field theory text [2] . Degrees of freedom in flat four dimensional spacetime are 
particles, classified by their spin. These degrees of freedom are carried by fields. If we wish 
to describe long range macroscopic forces, only bosonic fields will do, since fermionic fields 
cannot build up classical coherent states. By the spin statistics theorem, these bosonic fields 
must be of integer spin s = 0,1,2,3, etc. A field, ip, which carries a particle of mass m, 
will satisfy the Klein-Gordon equation (□ — 'm'^)ilj = 0, whose solution a distance r from a 
localized source goes like ~ ^e"'"''. Long range forces, those without exponential suppression, 
must therefore be described by massless fields, m = 0. 

Massless particles are characterized by how they transform under rotations transverse 
to their direction of motion. The transformation rule for bosons is characterized by an 
integer h > 0, which we call the helicity. For h = 0, such massless particles can be carried 
most simply by a scalar field 0. For a scalar field, any sort of interaction terms consistent 
with Lorentz invariance can be added, and so there are a plethora of possible self-consistent 
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interacting theories of spin particles. 

For helicities s > 1, the field must carry a gauge symmetry if we are to write interactions 
with manifest Lorentz symmetry and locality. For helicity 1, if we choose a vector field to 
carry the particle, its action is fixed to be the Maxwell action, so even without Maxwell, we 
could have discovered electromagnetism via these arguments. If we now ask for consistent 
self interactions of such massless particles, we are led to the problem of deforming the action 
(and possibly the form of the gauge transformations), in such a way that the linear form of 
the gauge transformations is preserved. These requirements are enough to lead us essentially 
uniquely to the non-abelian gauge theories, two of which describe the strong and weak forces 
[3]. 

Moving on to helicity 2, the required gauge symmetry is linearized general coordinate 
invariance. Asking for consistent self interactions leads essentially uniquely to GR and full 
general coordinate invariance [4, 5, 6, 7, 8, 9, 10] (see also chapter 13 of [2], which shows how 
helicity 2 implies the equivalence principle). For helicity > 3, the story ends, because there 
are no self interactions that can be written [11] (see also chapter 13 of [2], which shows that 
the scattering amplitudes for helicity > 3 particles vanish). 

This path is straightforward, starting from the principles of special relativity (Lorentz 
invariance), to the classification of particles and fields that describe them, and finally to their 
possible interactions. The path Einstein followed, on the other hand, is a leap of insight and 
has logical gaps; the equivalence principle and general coordinate invariance, though they 
suggest GR, do not lead uniquely to GR. 

General coordinate invariance is a gauge symmetry, and gauge symmetries are redun- 
dancies of description, not fundamental properties. In any system with gauge symmetry, 
one can always fix the gauge and eliminate the gauge symmetry, without breaking the phys- 
ical global symmetries (such as Lorentz invariance) or changing the physics of the system 
in any way. One often hears that gauge symmetry is fundamental, in electromagnetism for 
example, but the more correct statement is that gauge symmetry in electromagnetism is 
necessary only if one demands the convenience of linearly realized Lorentz symmetry and 
locality. Fixing a gauge will not change the physics, but the price paid is that the Lorentz 
symmetries and locality are not manifest. 

On the other hand, starting from a system without gauge invariance, it is always 
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possible to introduce gauge symmetry by putting in redundant variables. Often this can be 
very useful for studying a system, and can elucidate properties which are otherwise difficult 
to see. This is the essence of the Stiikelberg trick, which we will make use of extensively 
in our study of massive gravity. In fact, as we will see, this trick can be used to make any 
lagrangian invariant under general coordinate diffeomorhpisms, the same group under which 
GR is invariant. Thus general coordinate invariance cannot be the defining feature of GR. 

Similarly, the principle of equivalence, which demands that all mass and energy grav- 
itate with the same strength, is not unique to GR. It can be satisfied even in scalar field 
theories, if one chooses the interactions properly. For example, this can be achieved by it- 
eratively coupling a canonical massless scalar to its own energy momentum tensor. Such 
a theory in fact solves all the problems Einstein set out to solve; it provides a universally 
attractive force which conforms to the principles of special relativity, reduces to newtonian 
gravity in the non-relativistic limit, and satisfies the equivalence principle^. By introduc- 
ing diffeomorphism invariance via the Stiikelberg trick, it can even be made to satisfy the 
principle of general coordinate invariance. 

The real underlying principle of GR has nothing to do with coordinate invariance 
or equivalence principles or geometry, rather it is the statement: general relativity is the 
theory of a non-trivially interacting massless helicity 2 particle. The other properties are 
consequences of this statement, and the implication cannot be reversed. 

As a quantum theory, GR is not UV complete. It must be treated as an effective field 
theory valid at energies up to a cutoff at the Planck mass, Mp, beyond which unknown 
high energy effects will correct the Einstein-Hilbert action. As we will discuss in Section 
6.1, for a given background such as the spherical solution around a heavy source of mass M 
such as the Sun, GR has three distinct regimes. There is a classical linear regime, where 
both non-linear effects and quantum effects can be ignored. This is the regime in which r 
is greater than the Schwarzschild radius, r > r5 ~ jp-. For M the mass of the Sun, we 
have ~ 1 km, so the classical linear approximation is good nearly everywhere in the solar 
system. There is the quantum regime r < very near the singularity of the black hole, 
where the effective field theory description breaks down. Most importantly, there is a well 

■^This theory is sometimes known as the Einstein-Fokker theory, first introduced in 1913 by Nordstrom 
[12, 13], and later in a different form [14, 15]. It was even studied by Einstein when he was searching for a 
relativistic theory of gravity that embodied the equivalence principle [16]. 
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separated middle ground, a classical non-linear regime, < r < rs, where non-linearities 
can be summed up without worrying about quantum corrections, the regime which can be 
used to make controlled statements about what is going on inside a black hole. One of 
the challenges of adding a mass to the graviton, or any modification of gravity, is to retain 
calculable yet interesting regimes such as this. 

1.2 Modifying general relativity 

A theory of massive gravity is a theory which propagates a massive spin 2 particle. The 
most straightforward way to construct such a theory is to simply add a mass term to the 
Einstein-Hilbert action, giving the graviton a mass, m, in such a way that GR is recovered 
as m — 7- 0. This is a modification of gravity, a deformation away from the elegant theory of 
Einstein. Since GR is the essentially unique theory of a massless spin 2 degree of freedom, 
it should be remembered that modifying gravity means changing its degrees of freedom. 

Despite the universal consensus that GR is a beautiful and accurate theory, there has 
in recent years arisen a small industry of physicists working to modify it and test these 
modifications. When asked to cite their motivation, they more often than not point to 
supernova data [17, 18] which shows that the universe has recently started accelerating in 
its expansion. If GR is correct, there must exist some dark energy density, p ~ 10~^^ g/cm^. 
The simplest interpretation is that there is a constant term. A, in the Einstein-Hilbert action, 
which would give p ~ MpA. To give the correct vacuum energy, this constant has to take the 
small value A/Mp ~ 10^^^, whereas arguments from quantum field theory suggest a value 
much larger, up to order unity [19]. It is therefore tempting to speculate that perhaps GR is 
wrong, and instead of a dark energy component, gravity is modified in the infrared [20, 21], 
in such a way as to produce an accelerating universe from nothing. Indeed may modifications 
can be cooked up which produce these so-called self accelerating solutions. For example, one 
well-studied modification is to replace the Einstein-Hilbert lagrangian with F{R), a general 
function of the Ricci scalar [22, 23], which can lead to self-accelerating solutions [24, 25]. 
This modification is equivalent to adding an additional scalar degree of freedom. 

These cosmological reasons for studying modifications to gravity are often criticized on 
the grounds that they can only take us so far; the small value of the cosmological acceleration 
relative to the Planck mass must come from somewhere, and the best these modifications 



7 



can do is to shift the fine-tuning into other parameters (see [26] for an illustration in the 
F{R) I scalar-tensor case) . 

While it is true the small number must come from somewhere, there remains the 
hope that it can be put somewhere which is technically natural, i.e. stable to quantum 
corrections. Some small parameters, such as the ratio of the Higgs mass to the Planck mass 
in the standard model, are not technically natural, whereas others, like small fermion masses, 
are technically natural, because their small values are stable under quantum corrections. A 
rule of thumb is that a small parameter is technically natural if there is a symmetry that 
appears as the small parameter is set to zero. When this is the case, symmetry protects a 
zero value of the small parameter from quantum corrections. This means corrections due 
to the small parameter must be proportional to the parameter itself. In the case of small 
fermion masses, it is chiral symmetry that appears, whereas in the case of the Higgs mass 
and the cosmo logical constant, there is no obvious symmetry that appears. 

Of course there is no logical inconsistency with having small parameters, technically 
natural or not, and nature may explain them anthropically [27], or may just employ them 
without reason. But as practical working physicists, we hope that it is the case that a small 
parameter is technically natural, because then there is a hope that perhaps some classical 
mechanism can be found that drives the parameter towards zero, or otherwise explains its 
small value. If it is not technically natural, any such mechanism will be much harder to find 
because it must know about the quantum corrections in order to compensate them. 

One does not need a cosmological constant problem, however, to justify studying mod- 
ifications to GR. There are few better ways to learn about a structure, whether it's a sports 
car, a computer program or a theory, than to attempt to modify it. With a rigid theory 
such as GR, there is a level of appreciation that can only be achieved by witnessing how 
easily things can go badly with the slightest modification. In addition, deforming a known 
structure is one of the best ways to go about discovering new structures, structures which 
may have unforeseen applications. 

One principle that comes to play is the continuity of physical predictions of a theory 
in the parameters of the theory. Surely, we should not be able to say experimentally, given 
our finite experimental precision, that a parameter of nature is exactly mathematically zero 
and not just very small. If we deform GR by a small parameter, the predictions of the 



8 



deformed theory should be very close to GR, to the extent that the deformation parameter 
is small. It follows that any undesirable pathologies associated with the deformation should 
cure themselves as the parameter is set to zero. Thus, we uncover a mechanism by which 
such pathologies can be cured, a mechanism which may have applications in other areas. 

Massive gravity is a well developed case study in the infrared modification of gravity, 
where all of these points are nicely illustrated. Purely from the consideration of degrees of 
freedom, it is a natural modification to consider, since it amounts to simply giving a mass to 
the particle which is already present in GR. In another sense, it is less minimal than F{R) 
or scalar tensor theory, which adds a single scalar degree of freedom, because to reach the 5 
polarizations of the massive graviton we must add at least 3 degrees of freedom beyond the 
2 of the massless graviton. 

With regard to the cosmological constant problem, there is the possibihty of a tech- 
nically natural explanation. The deformation parameter is m, the graviton mass, and GR 
should be restored as m — > 0. The force mediated by a massive graviton has a Yukawa 
profile ~ ^e""^'', which drops off from that of a massless graviton at distances r > ^, so one 
could hope to explain the acceleration of the universe without dark energy by choosing the 
graviton mass to be of order the Hubble constant m ^ H. Of course, this does not eliminate 
the small cosmological constant, which reappears as the ratio m/Mp. But there is now hope 
that this is a technically natural choice, because deformation by a mass term breaks the 
gauge symmetry of GR, which is restored in the limit m — > 0. As we will see, a small m 
is indeed protected from quantum corrections (though as we will see, there are other issues 
that prevent this, at our current stage of understanding, from being a completely satisfactory 
reahzation of a technically natural cosmological constant). 

There are also interesting lessons to be learned regarding the continuity of physical 
predictions. The addition of a mass term is a brutality upon the beautiful structure of 
GR, and does not go unpunished. Various pathologies appear, which arc representative of 
common pathologies associated with any infrared modification of gravity. These include 
strong classical non-linearities, ghost-like instabilities, and a very low cutoff, or region of 
trustability, for the resulting quantum effective theory. In short, modifying the infrared 
often messes up the UV. New mechanisms also come into play, because the extra degrees of 
freedom carried by the massive graviton must somehow decouple themselves as m — > to 
restore the physics of GR. 



9 



The study of the massless hmit leads to the discovery of the Vainshtein mechanism, by 
which these extra degrees of freedom hide themselves at short distances using non-linearities. 
This mechanism has already proven useful for model builders who have long-range scalars, 
such as moduli from the extra dimensions of string theory, that they want to shield from 
local experiments that would otherwise rule them out. 

1.3 History and outline 

The possibility of a graviton mass has been studied off and on since 1939, when Fierz and 
Pauli [28] first wrote the action describing a free massive graviton. Following this, not much 
occurred until the early 1970's, when there was a flurry of renewed interest in quantum 
field theory. The linear theory coupled to a source was studied by van Dam, Veltman, and 
Zakharov [29, 30], who discovered the curious fact that the theory makes predictions different 
from those of linear GR even in the limit as the graviton mass goes to zero. For example, 
massive gravity in the m — )■ limit gives a prediction for light bending that is off by 25 
percent from the GR prediction. The linear theory violates the principle of continuity of 
the physics in the parameters of the theory. This is known as the vDVZ discontinuity. The 
discontinuity was soon traced to the fact that not all of the degrees of freedom introduced 
by the graviton mass decouple as the mass goes to zero. The massive graviton has 5 spin 
states, which in the massless limit become the 2 helicity states of a massless graviton, 2 
helicity states of a massless vector, and a single massless scalar. The scalar is essentially 
the longitudinal graviton, and it maintains a finite coupling to the trace of the source stress 
tensor even in the massless limit. In other words, the massless limit of a massive graviton is 
not a massless graviton, but rather a massless graviton plus a coupled scalar, and the scalar 
is responsible for the vDVZ discontinuity. 

If the linear theory is accurate, then the vDVZ discontinuity represents a true physical 
discontinuity in predictions, violating our intuition that physics should be continuous in 
the parameters. Measuring the light bending in this theory would be a way to show that 
the graviton mass is mathematically zero rather than just very small. However, the linear 
theory is only the start of a complete non-linear theory, coupled to all the particles of the 
standard model. The possible non-linearities of a real theory were studied several years later 
by Vainshtein [31], who found that the nonlinearities of the theory become stronger and 
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stronger as the mass of the graviton shrinks. What he found was that around any massive 
source of mass M, such as the Sun, there is a new length scale known as the Vainshtein radius, 
Ty ~ f j • At distances r < ry, non-linearities begin to dominate and the predictions 
of the linear theory cannot be trusted. The Vainshtein radius goes to infinity as m — )■ 0, so 
there is no radius at which the linear approximation tells us something trustworthy about the 
massless limit. This opens the possibility that the non-linear effects cure the discontinuity. 
To have some values in mind, if we take M the mass of the Sun and m a very small value, 
say the Hubble constant m ~ 10^'^^ eV, the scale at which we might want to modify gravity 
to explain the cosmological constant, we have ry ~ 10^® km, about the size of the Milky 
Way. 

Later the same year, Boulware and Deser [32] studied some specific fully non-linear 
massive gravity theories and showed that they possess a ghost-like instability. Whereas the 
linear theory has 5 degrees of freedom, the non-linear theories they studied turned out to 
have 6, and the extra degree of freedom manifests itself around non-trivial backgrounds as 
a scalar field with a wrong sign kinetic term, known as the Boulware- Deser ghost. 

Meanwhile, the ideas of effective field theory were being developed, and it was real- 
ized that a non-renormalizable theory, even one with apparent instabilities such as massive 
gravity, can be made sense of as an effective field theory, valid only at energies below some 
ultraviolet cutoff scale A. In 2003, Arkani-Hamed, Georgi and Schwartz [33] brought to 
attention a method of restoring gauge invariance to massive gravity in a way which makes it 
very simple to see what the effective field theory properties are. They showed that massive 
gravity generically has a maximum UV cutoff of A5 = (Mpm'^Y^^ . For Hubble scale gravi- 
ton mass, this is a length scale A5 ^ ~ 10^^ km. This is a very small cutoff, parametrically 

smaller than the Planck mass, and goes to zero as m — )■ 0. Around a massive source, the 

/ \ 1/3 

quantum effects become important at the radius rg = i-^) xT, which is parametrically 
larger than the Vainshtein radius at which non-linearities enter. For the Sun, rg ~ lO^'* km. 
Without finding a UV completion or some other re-summation, there is no sense in which 
we can trust the solution inside this radius, and the usefulness of massive gravity is limited. 
In particular, since the whole non-linear regime is below this radius, there is no hope to 
examine the continuity of physical quantities in m and explore the Vainshtein mechanism in 
a controlled way. On the other hand, it can be seen that the mass of the Boulware-Deser 
ghost drops below the cutoff only when r < rg, so the ghost is not really in the effective 
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theory at all and can be consistently excluded. 

Putting aside the issue of quantum corrections, there has been continued study of the 
Vainshtein mechanism in a purely classical context. It has been shown that classical non- 
linearities do indeed restore continuity with GR in certain circumstances. In fact, the ghost 
degree of freedom can play an essential role in this, by providing a repulsive force in the 
non-linear region to counteract the attractive force of the longitudinal scalar mode. 

By adding higher order graviton self-interactions with appropriately tuned coefficients, 
it is in fact possible to raise the UV cutoff of the theory to A3 = (Mpm^)^/^, corresponding 
to roughly Aj^ ~ 10^ km. In 2010, the complete action of this theory in a certain decoupling 
limit was worked out by de Rham and Gabadadze [34], and they show that, remarkably, 
it is free of the Boulware-Deser ghost. Recently, it was shown that the complete theory is 
free of the Boulware-Deser ghost. This A3 theory is the best hope of realizing a useful and 
interesting massive gravity theory. 

The subject of massive gravity also naturally arises in extra-dimensional setups. In a 
Kaluza-Klein scenario such as GR in 5d compactified on a circle, the higher Kaluza-Klein 
modes are massive gravitons. Brane world setups such as the Dvali-Gabadadze-Porrati 
(DGP) model [35] give more intricate gravitons with resonance masses. The study of such 
models has complemented the study of pure 4d massive gravity and has pointed towards 
new research directions. 

The major outstanding question is whether it is possible to UV extend the effective 
field theory of massive gravity to the Planck scale, and what this UV extension may look 
like. This would provide a solution to the problem of making the small cosmological constant 
technically natural, and is bound to be an interesting theory its own right (the analogous 
question applied to massive vector bosons leads to the discovery of the Higgs mechanism 
and spontaneous symmetry breaking). In the case of massive gravity, there are indications 
that a UV completion may not have a local Lorentz invariant form. Another long shot, if 
UV completion can be found, would be to take the m — )■ limit of the completion and hope 
to obtain a UV completion to ordinary GR. 

As this review is focused on the theoretical aspects of Lorentz invariant massive gravity, 
we will not have much to say about the large literature on Lorentz-violating massive gravity. 
We will also not say much about the experimental search for a graviton mass, or what the 
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most likely signals and search modes would be. There has been much work on these areas, 
and each could be the topic of a separate review. 

Conventions: Often we will work in an arbitrary number of dimensions, just because 
it is easy to do so. In this case, D signifies the number of spacetime dimension and we stick 
to -D > 3. d signifies the number of space dimensions, d = D — 1. We use the mostly plus 
metric signature convention, rj^^i, = (—,+,+,+,■■■). Tensors are symmetrized and anti- 
symmetrized with unit weight, i.e T(^,,) = ^ {T^^ + T^^), T[^^] = ^ {T^^ - Tj,^). The reduced 
4d Planck mass is Mp = (g^g)i/2 ~ 2.43 x 10^^ GeV. Conventions for the curvature tensors, 
covariant derivatives and Lie derivatives are those of Carroll [36]. 



2 The free Fierz-Pauli action 



We start by displaying an action for a single massive spin 2 particle in fiat space, carried by 
a symmetric tensor field /i^j,, 

S = j d''x-^dxh^,d^h^'' + d^K^d''h^^-d,M^''d,h + ^dxhd^h-^m\h^,h^''-h^). (2.1) 

This is known as the Fierz-Pauli action [28]. Our point of view will be to take this action 
as given and then show that it describes a massive spin 2. There are, however, some (less 
than thorough) ways of motivating this action. To start with, the action above contains all 
possible contractions of two powers of h, with up to two derivatives. The two derivative 
terms, those which survive when m = 0, are chosen to match exactly those obtained by 
linearizing the Einstein-Hilbert action. The m = terms describe a massless helicity 2 
graviton and have the gauge symmetry 

Sh^u = d^i, + d,i^, (2.2) 

for a spacetime dependent gauge parameter ^^(x). This symmetry fixes all the coefficients 
of the two- derivative part of (2.1), up to an overall coefficient. The mass term, however, 
violates this gauge symmetry. The relative coefficient of —1 between the h"^ and h^^^h^'^ 
contractions is called the Fierz-Pauli tuning, and it not enforced by any known symmetry. 

However, the only thing that needs to be said about this action is that it describes 
a single massive spin 2 degree of freedom of mass m. We will show this explicitly in what 



13 



follows. Any deviation from the form (2.1) and the action will no longer describe a single 
massive spin 2. For example, violating the Fierz-Pauli tuning in the mass term by changing 
to -Im^ihf,^'' - (1 - a)h^) for a gives an action describing a scalar ghost (a scalar 
with negative kinetic energy) of mass = ^^^f^^^, in addition to the massive spin 2. For 
small a, the ghost mass squared goes like ~ ^. It goes to infinity as the Fierz-Pauli tuning 
is approached, rendering it non-dynamical. Violating the tuning in the kinetic terms will 
similarly alter the number of degrees of freedom, see [37] for a general analysis. 

There is a method of constructing lagrangians such as (2.1) to describe any given spin. 
See for example the first few chapters of [2], the classic papers on higher spin lagrangians 
[38, 39], and the reviews [40, 41]. 



2.1 Hamiltonian and degree of freedom count 

We will begin our study of the Fierz-Pauli action (2.1) by casting it into hamiltonian form 
and counting the number of degrees of freedom. We will show that it propagates — i 

degrees of freedom in D dimensions (5 degrees of freedom for D = A), the right number for 
a massive spin 2 particle. 

We start by Legendre transforming (2.1) only with respect to the spatial components 
hij. The canonical momenta are'^ 



TT, 



- ^ 

ohij 



hij — hkkSij — 2(9(i/ij)o + 2dkhQk5ij. (2.3) 



Inverting for the velocities, we have 

hij = TTij - ^ _ ^ T^kkSij + 2d(ihj)Q. (2.4) 

In terms of these hamiltonian variables, the Fierz-Pauli action (2.1) becomes 
S = / - H + 2fto. (9...) + rnV, + ft„„ (v=A.. - 3^,h, - m^ft..) . (2.5) 



■^Note that canonical momenta can change under integrations by parts of the time derivatives. We have 
fixed this ambiguity by integrating by parts such as to remove all derivatives from h^i and ft-oo- 
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where 



2 '^ 2D -2 '' 

+ -^dkhijdkhij — dihjkdjhik + dihijdjh^k — -^dihjjdihf^k + -m^{hijhij — h\ 



(2.6) 



First consider the case m = 0. The time-hke components h^i and /iqo appear hnearly 
multiphed by terms with no time derivatives. We can interpret them as Lagrange multiphers 
enforcing the constraints djUij = and V'^ha — didjhij = 0. It is straightforward to check 
that these are first class constraints, and that the hamiltonian (2.6) is first class. Thus (2.5) 
is a first class gauge system. For D = 4, the hij and TTjj each have 6 components, because 
they are 3x3 symmetric tensors, so together they span a 12 dimensional (for each space 
point) phase space. We have 4 constraints (at each space point), leaving an 8 dimensional 
constraint surface. The constraints then generate 4 gauge invariances, so the gauge orbits are 
4 dimensional, and the gauge invariant quotient by the orbits is 4 dimensional (see [42] for an 
introduction to constrained hamiltonian systems, gauge theories, and the terminology used 
here). These are the two polarizations of the massless graviton, along with their conjugate 
momenta. 

In the case m 7^ 0, the hoi are no longer Lagrange multipliers. Instead, they appear 
quadratically and are auxiliary variables. Their equations of motion yield 

hoi = — KdjTCij, (2.7) 

which can be plugged back into the action (2.5) to give 

s = / d-x - n + A„„ - - , (2,8) 



where 



"H = ^TTij - ^jyZ2^'^' ^dkhijdkhij - dihjkdjhik + dihijdjhkk - ^dihjjdihkk 

1 o,. . .0. 1 



+ ^m'^{hijhij - hi) + {djTTijf . 

III 

(2.9) 



The component /iqo remains a Lagrange multiplier enforcing a single constraint C 



—V'^hii + didjhij + rn^hii = 0, but the hamiltonian is no longer first class. One secondary 
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constraint arises from the Poisson bracket with the hamiltonian H = J d'^x H, namely 
{H, C}pB = j^m'^7!-ii + didjUij. The resulting set of two constraints is second class, so there 
is no longer any gauge freedom. For D = 4 the 12 dimensional phase space has 2 constraints 
for a total of 10 degrees of freedom, corresponding to the 5 polarizations of the massive 
graviton and their conjugate momenta. 

Note that the Fierz-Pauli tuning is crucial to the appearance of /iqo as a Lagrange 
multiplier. If the tuning is violated, then /iqo appears quadratically and is an auxiliary 
variable, and no longer enforces a constraint. There are then no constraints, and the full 12 
degrees of freedom in the phase space are active. The extra 2 degrees of freedom are the 
scalar ghost and its conjugate momentum. 



2.2 Free solutions and graviton mode functions 

We now proceed to find the space of solutions of (2.1), and show that it transforms as 
a massive spin 2 representation of the Lorentz group, showing that the action propagates 
precisely one massive graviton. The equations of motion from (2.1) are 

S S 

= nhf,^ - dxdf,h\ - dxd^h^ + r]^^dxd^h^'' + d^d^h - r]^^Uh - m^{h^^ - r]^^h) = 0. 

(2.10) 

Acting on (2.10) with d'^, we find, assuming ^ 0, the constraint d^h^^ — d^h. 
Plugging this back into the equations of motion, we find D/i^i, — d^d^h — vn?{hy^y — ri^yK) = 0. 
Taking the trace of this we find h = 0, which in turn implies d^h^i,=0. This, along with 
h = applied to the equations of motion (2.10), gives (□ — m'^)h^y = 0. 

Thus the equations of motion (2.10) imply the three equations, 

(□ - m^)/^^, = 0, a^V = 0, h = 0. (2.11) 

Conversely, it is straightforward to see that these three equations imply the equations of 
motion (2.10), so (2.11) and (2.10) are equivalent. The form (2.11) makes it easy to count 
the degrees of freedom as well. For D = 4, the first of (2.11) is an evolution equation for the 
10 components of the symmetric tensor /i^^,, and the last two are constraint equations on the 
initial conditions and velocities of /i^j^. The last determines the trace completely, killing one 
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real space degree of freedom. The second gives 4 initial value constraints, and the vanishing 
of its time derivative, i.e. demanding that it be preserved in time, implies 4 more initial 
value constraints, thus killing 4 real space degrees of freedom. In total, we are left with the 
5 real space degrees of freedom of a 4 dimensional spin 2 particle, in agreement with the 
hamiltonian analysis of Section 2.1. 

The first equation in (2.11) is the standard Klein-Gordon equation, with the general 
solution 

h^'^ix) = I {h^''{p)e'P-' + h^''*{p)e-^^-^) . (2.12) 



Here p are the spatial momenta, oOp = ^Jp"^ + rn^, and the D-momenta are on shell so 
that = {cOp, p). 

Next we expand the Fourier coefficients h^'^{p) over some basis of symmetric tensors, 
indexed by A, 

/i'^'^(p) = ap,Ae'^'^(p,A). (2.13) 

We will fix the momentum dependence of the basis elements e^'^{Pi A) by choosing some 
basis e^'^(k. A) at the standard momentum = (m, 0, 0, 0, . . .) and then acting with some 
standard boost"^ L{p), which takes k into p, p^ = L{p)'^j^k'^. This standard boost will choose 
for us the basis at p, relative to that at k. Thus we have 

e^'^ip, A) = LipYM^^^^K A). (2.15) 

Imposing the conditions d^h^'^ = and /i = on (2.12) then reduces to imposing 

k^e'^-'ik, A) = 0, v^.e^'^ik, A) = 0. (2.16) 
^We choose the standard boost to be 

i^-(p) = % + ^(7-l)pV, 

IpI 

L\ip) = L°,{p) = ^^Vl^, 

L\{P)^1, (2.14) 

where 

7 



TO TO 

is the usual relativistic 7. See chapter 2 of [2] for discussions of this standard boost and general representation 
theory of the Poincare group. 
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The first says that e'*'^(k, A) is purely spatial, i.e. e^'^(k, A) = 0. The second says that it 
is traceless, so that ?j(k. A) = also. Thus the basis need only be a basis of symmetric 
traceless spatial tensors, A = 1, . . . , ^^^til _ i demand that the basis be orthonormal, 

e^'^(k,A)e;,(k,A') = 5,,,. (2.17) 

This basis forms the symmetric traceless representation of the rotation group SO{d), 
which is the little group for a massive particle in D dimensions. If W^^, is a spatial rotation, 
we have 

R^^M'^re^'iK A') = <e^^(k. A'), (2.18) 

where R^'^ is the symmetric traceless tensor representation of R^^,. We are free to use any 
other basis e''^(k. A), related to the e^''(k. A) by 

€^"'{k,X) = B\€'"'{k,X'), (2.19) 

where B is any unitary matrix. 

Given a particular spatial direction, with unit vector A;*, there is an SO{d — 1) subgroup 
of the little group SO{d) which leaves A;* invariant, and the symmetric traceless rep of SO{d) 
breaks up into three reps of SO{d — 1), a scalar, a vector, and a symmetric traceless tensor. 
The scalar mode is called the longitudinal graviton, and has spatial components 



After a large boost in the A;* direction, it goes like ei ~ /m^. As we will see later, in the 
massless limit, or large boost limit, this mode is carried by a scalar field, which generally 
becomes strongly coupled once interactions are taken into account. The vector modes have 
spatial components 

^i, = ^kHt (2.21) 

and after a large boost in the /c* direction, they go like ~ p/m. In the massless limit, these 
modes are carried by a vector field, which decouples from conserved sources. The remaining 
linearly independent modes are symmetric traceless tensors with no components in the 
directions, and form the symmetric traceless mode of SO{d— 1). They are invariant under a 
boost in the /c' direction, and in the massless limit, they are carried by a massless graviton. 
In the massless hmit, we should therefore expect that the extra degrees of freedom of the 
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massive graviton should organize themselves into a massless vector and a massless scalar. 
We will see later explicitly how this comes about at the lagrangian level. 

Upon boosting to p, the polarization tensors satisfy the following properties: they are 
transverse to and traceless, 

p^e'^'^(p,A)=0, r/^.e^^(p,A) = 0, (2.22) 

and they satisfy orthogonality and completeness relations 

e^^(p,A)e;,(p,A') = 5AA', (2.23) 
^ e^''(p, A)e*"^(p, A) = i(p/^"P'^/^ + pM/3p^") _ _i_pM^pa/3^ (2.24) 

A 

where P"^ = rj^^ + The right hand side of the completeness relation (2.24) is the 

projector onto the symmetric and transverse traceless subspace of tensors, i.e. the identity 
on this space. We also have the following symmetry properties in p, which can be deduced 
from the form of the standard boost (2.14), 

e^^(-p,A) = e^^(p,A), z, j = 1, 2, . . . , (2.25) 
e°^(-p,A) = -eO*(p,A), z = l,2,...,rf (2.26) 
e°°(-p,A) = e°°(p,A). (2.27) 

The general solution to (2.10) thus reads 

h^^{x) = [ Yl ^p^^'^iP^ ^y'"" + «;,Ae*"^(P, A)e-^-. (2.28) 



A 



The solution is a general linear combination of the following mode functions and their 
conjugates 

«^:,(x)^-=:i==e^'^(p,A)e^^-, A = l,2,...,ci. (2.29) 
A/(27r)''2a;p 

These are the solutions representing gravitons, and they have the following Poincare trans- 
formation properties 

<A(^-«)=<:A(^)e-'"-", (2-30) 



AVAV<A'(A'^a:) = ./^W^(A,J>)a'a<;,,(x) 

V 



(2.31) 
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where W{A,p) = L ^{Ap)AL{p) is the Wigner rotation, and W{A,p)x'x is its spin 2 rep, 
R^^ — 7- {B^^RB)x'x.^ Thus the gravitons are spin 2 solutions. 

In terms of the modes, the general solution reads 

h^'^ix) = [d'pJ2 Ka<:a(^) + «;,A<A (^)) • (2-32) 
The inner (symplectic) product on the space of solutions to the equations of motion is. 



{h,h') = i I (T^ h'"'*{x)doh'^^{x) 



(2.33) 



The u functions are orthonormal with respect to this product, 

Ka,V,a') = 5''(p-p')5aa', (2.34) 

K,A,«kA') = -'5'(p-p')'5AA', (2.35) 

(np,A,^i;,A') = 0, (2.36) 

and we can use the product to extract the a and a* coefficients from any solution h^i^lx), 

ap,A = (Mp,A, h), (2.37) 

<A = -K,A,/^)- (2.38) 



^We show the Lorentz transformation property as follows 

= [L{Ap) {L-\Ap)AL{p))]^^, [L{Ap) {L-'{Ap)AL{p))]\, e{k,\r'^' ^^^'^ 
= [L{Ap)W {Kp)T^, [L{Ap)W{A,p)]\,e{kAY''^'e^^^-\ 
The little group element is a spatial rotation. For any spatial rotation R^^, we have 

i?^^,i?'^,,e^'^'(fc,A) = i?^^,i?'',,i?^;e-^''''(k, A') ^ B^^.e^'^ {k, X") 
= {B-'RB)\ef"'{k,X'). 

Plugging back into the above, 

^WiA,p)^y''iAp,X')e^^P--, 
where W is the spin 2 representation of the httle group in a basis rotated by B, = B^^RB. 

20 



In the quantum theory, the a and a* become creation and annihilation operators which 
satisfy the usual commutations relations and produce massive spin 2 states. The fields 
hij and their canonical momenta 7?,^, constructed from the a and a*, will then automatically 
satisfy the Dirac algebra and constraints of the hamiltonian analysis of Section 2.1, providing 
a quantization of the system. Once interactions are taken into account, external lines of 
Feynman diagrams will get a factor of the mode functions (2.29). 

2.3 Propagator 

Integrating by parts, we can rewrite the Fierz-Pauli action (2.1) as 

S = Jd^'x hi^,0^^'''f'Kp, (2.39) 

where 

O^^^p = [vtv l - r/'-r^^^) (□ - m^) - 2a(^9(„r/% + d^^d'^r^^^ + d^d^r^^'', (2.40) 
is a second order differential operator satisfying 

In terms of this operator, the equations of motion (2.10) can be written simply as = 
0^"''''^ha,3 = 0. 

To derive the propagator, we go to momentum space, 
O'^l^id ^ip) = - (vtv'^l - r'^T]^^^ ip' + m^) + 2p('^p(„r/'^)) - p^'^r/^^ - Po^pp^". (2.42) 
The propagator is the operator T^ap^aX with the same symmetries (2.41) which satisfies 

0^^'"^2^a/3,.A = ^(5^55: + «)- (2.43) 

The right side is the identity operator on the space of symmetric tensors. 
Solving (2.43), we find 



T) - ~^ 



2 [PaaPpX + PaxP/Sa) — -^Pa/sPaX 



(2.44) 



where Pa/s = r]a/5 + 
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In the interacting quantum theory, internal hues with momentum p will be assigned 

2 

this propagator, which for large p behaves as ~ This growth with p means we cannot 
apply standard power counting arguments (like those of chapter 12 of [2]) to deduce the 
renormalizability properties or strong coupling scales of a theory. We will see later how to 
overcome this difficulty by rewriting the theory in a way in which all propagators go like 
~ ^ at high energy. 

Massless propagator 

The massive graviton propagator (2.44) may be compared to the propagator for the case 
m = 0. For m = 0, the action becomes 

Sm=o = j d^'x hi^^e^^^^'^^Kp, (2.45) 

where the kinetic operator is 

= <^"'Jm=o = "^1 - ^'''^"/^) ° - 25^'^("^1) + ^'^''^o.p + d^dp^^r (2.46) 
This operator has the symmetries (2.41). Acting on a symmetric tensor Z^j^^ it reads 

e^^'"'^^^;^ = nZ^"^ - T^^'^UZ - 2d'^''daZ^^'' + d^d^Z + r/^^S^S/jZ"^. (2.47) 

The m = action has the gauge symmetry (2.2), and the operator (2.46) is not invert- 
ible. Acting with it results in a tensor which is automatically transverse, and it annihilates 
anything which is pure gauge 

(e^'^'^^Z^^j) = 0, e''"'''^ (a,e/3 + 9^3^) = 0. (2.48) 

To find a propagator, we must fix the gauge freedom. We choose the Lorenz gauge 
(also called harmonic, or de Bonder gauge), 

d^h^u - \d,h = 0. (2.49) 

We can reach this gauge by making a gauge transformation with chosen to satisfy = 
- [d^h^^ - \d^h) . This condit ion fixes the gauge only up to gauge transformations with 
parameter satisfying = 0. In this gauge, the equations of motion simphfy to 

□V - \n^.unh = 0. (2.50) 
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The solutions to this equation which also satisfy the gauge condition (2.49) are the Lorenz 
gauge solutions to the original equations of motion. 

To the lagrangian of (2.45) we add the following gauge fixing term 

2 



GF 



-dnh 
2 



(2.51) 



Quantum mechanically, this results from the Fadeev-Popov gauge fixing procedure. Clas- 
sically, we may justify it on the grounds that the equations of motion obtained from the 
action plus the gauge fixing term are the same as the gauge fixed equations of motion (2.50). 
The gauge condition itself, however, is not obtained as an equation of motion, and must be 
imposed separately. We have 



£ + £, 



GF 



-h^^Dh'"' - -hDh, 
2 ^ 4 



(2.52) 



whose equations of motion are (2.50). 

We can write the gauge fixed lagrangian as £ + £gf = \h^uO^'^'°'^hai3, where 



QIJ.U, 



,al3 



□ 



2 2 



Going to momentum space and inverting, we obtain the propagator. 



V. 



al3,crX 



—I 
p2 



1 / N 1 



(2.53) 



(2.54) 



which satisfies the relation (2.43) with O in place of O. This propagator grows like ~ ^ 
at high energy. Comparing the massive and massless propagators, (2.54) and (2.44), and 
ignoring for a second the terms in (2.44) which are singular as m — )■ 0, there is a difference 
in coefficient for the last term, even as m — 0. For D = A, it is 1/2 vs. 1/3. This is the first 
sign of a discontinuity in the m — )■ limit. 



3 Linear response to sources 



We now add a fixed external symmetric source T^^^{x) to the action (2.1) , 

S = j d''x-^dxh^,d^h^''+d^h,xd''h^^^-d^h^''d,h+^dxhd^h-^m\h^,h^''-h^ 

(3.1) 
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Here k = Mp ^ is the coupling strength to the source^. 
The equations of motion are now sourced by Tfj^^, 

(3.2) 

In the case m = 0, acting on the left with d'^ gives identically zero, so we must have the 
conservation condition d^T^y = if there is to be a solution. For m 7^ 0, there is no such 
condition. 



3.1 General solution to the sourced equations 

We now find the retarded solution of (3.2), to which the homogeneous solutions of (2.2) can 
be added to obtain the general solution. Acting on the equations of motion (3.2) with d^, 
we find, 

K 

(3.3) 



Plugging this back into (3.2), we find 

Where ddT is short for the double divergence d^d^T^'^ . Taking the trace of this we find 



(3.4) 



Applying this to (3.3), we find 



K 



T A d'^T ,y 

m"^ D — I 



m^{D - 1) 

which when applied along with (3.4) to the equations of motion, gives 

1 



(3.5) 



{d' -m')h 



— K 



T, 



D - 1 



(3.6) 



+ - 



K 



1 



D 



-h,, + (Z}-2) 



ddT 



2 sc 



"The normalizations chosen here are in accord with the general relativity definition T^'^ = -^j^j^ 
well as the normalization Sg^^, — 2Kh^^. 



as 
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Thus we have seen that the equations of motion (3.2) imply the following three equations, 



— K 



T 



+ - 



dadu 



T 



1 



K 



K 



h 



—d,T+—d^T, 
1 



D - 1 

K D-2 



r^,. + (D-2)^ 



ddT 



dyddT, 



-T 



K D 



ddT. 



(3.7) 



Conversely, it is straightforward to see that these three equations imply the equations of 
motion (3.2). 



Taking the first equation of (3.7) and tracing, we find (□— 



h + 



K D-2 



0. Under the assumption that (9^ — m^)/ = ^ / = for any function /, the third equa- 
tion is implied. This will be the case with good boundary conditions, such as the retarded 
boundary conditions we impose when we are interested in the classical response to sources. 
The second equation of (3.7) can also be shown to follow under this assumption, so that 
we may obtain the solution by Fourier transforming only the first equation of (3.7). This 
solution can also be obtained by applying the propagator (2.44) to the Fourier transform of 
the source. 

Despite the absence of gauge symmetry, we will often be interested in sources which 
are conserved anyway, d^T'^'^ = 0. When the source is conserved, and under the assumptions 
in the paragraph above, we are left with just the equation. 



ddT 



{d^ -m^)h 



— K 



T 

D 



The general solution for a conserved source is then. 



(2^ 



^tpx _ 



p2 -|- J77,2 



D - 1 



dndy 



T 



(3.8) 



im? 



T{p) 



(3.9) 



where T^'^{p) is the Fourier transform of the source, T^'^{p) = J d^x e '^'^^ T^'^{x). To get 
the retarded field, we should integrate above the poles in the p^ plane. 
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3.2 Solution for a point source 



We now specialize to four dimensions so that n = 1/Mp, and we consider as source the stress 
tensor of a mass M point particle at rest at the origin 



Note that this source is conserved. For this source, the general solution (3.9) gives 



(3.10) 



hooix) 
hoi{x) 
hij{x) 



2M 



P 



3MpJ (27r)3 p2 + m2 
0, 



M 



1 



3Mp J (27r)3 p2 + m2 



+ 



PiPj 



(3.11) 



Using the formulae 



(27r)3 
d'p 
(27r)3 



p2 + 



1 e"""^ 
47r r ' 



p2 + 



1 e"'" 
An r 



(27r)3 

— (1 + mr)6ij j(3 + 3mr + m?r'^)xiXj 



p2 _|_ ^2 



(3.12) 



where r = y/xjXi, we have 

, , , 2M 1 e-"^^ 

hoi{x) 
hij{x) 



3Mp 4:71 r 
0, 

M 1 e""^ 
3Mp47r~ 



2^2 



1 + mr + m r 



- 



(3 + 3mr + m r )xiX 



(3.13) 



Note the Yukawa suppression factors e characteristic of a massive field. 

For future reference, it will be convenient to record these expressions in spherical coordi- 
nates for the spatial variables. Using the formula [F{r)5ij + G{r)xiXj\ dx^dx^ = {F{r) + r'^G{r)) dr'^+ 
F{r)r'^dVL^ to get to spherical coordinates we find 



h^^dx'^dx'' = -B{r)df + C{r)dr' + A{r)r^dn^ 



(3.14) 
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where 

2M 1 e""''' 

B{r) = 
C{r) = 



3Mp 47r r ' 
2M 1 e-™'^l + mr 



3Mp Att r m?r'^ ' 
Ml e"'"'' 1 + mr + m^r^ 



In the hmit r <^l/m these reduce to 

^ 2M 1 

±5 r = 



3Mp 47rr 
Cir) = 1 



Corrections are suppressed by powers of mr. 



Solution for the massless graviton 



For the purposes of comparison, we will compute the point source solution for the massless 
case as well. We choose the Lorenz gauge (2.49). In this gauge, the equations of motion 
simplify to 

□/i„,,y - -ViuyOh = -KT„,,y. (3.17) 



Taking the trace, we find Oh 



nh 



j^kT, and upon substituting back, we get 
1 



T, 



(3.18) 



D-2 

This equation, along with the Lorenz gauge condition (2.49), is equivalent to the original 
equation of motion in Lorenz gauge. 

Taking d'^ on (3.17) and on its trace, using conservation of T^^^, and comparing, we 
have D{d'^h^^ — ^d^h) = 0, so that the Lorentz condition is automatically satisfied when 
boundary conditions are satisfied with the property that □/ = ^ / = for any function 
/, as is the case when we impose retarded boundary conditions. We can then solve 3.17 by 
Fourier transforming. 



W7 



D 



1 

D-2 



(3.19) 
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where T^^{p) = J d^x e T'^^(x), is the Fourier transform of the source. To get the 
retarded field, we should integrate above the poles in the p° plane. 

Now we specialize to -D = 4, and we consider as a source the point particle of mass M 
at the origin (3.10). For this source, the general solution (3.19) gives 

M f d'p 1 M 1 



hmix) = —— / -—76 



""^■^^ 2MpJ (27r)3^ p2 2Mp47rr' 
hoi{x) = 0, 



h,Jx) = ^ [ /ii^e^P"^5,, = ^7^5i,-. (3.20) 
^ 2Mp J (27r)3 p2 2Mp47rr ^ ^ ^ 

For later reference, we record this result in spherical spatial coordinates as well. Us- 
ing the formula [F{r)6ij + G{r)xiXj]dx'^dx^ = {F{r) + r'^G{r)) dr"^ + F{r)r'^dQ'^ to get to 
spherical coordinates we find 

h^^dx^dx" = -B{r)dt^ + C{r)dr'^ + A{r)r'^dn'^, (3.21) 

where 

Ml 

Bir) = 



2Mp47ir 
C(r) = 

...Ml , , 

A(r) = . 3.22 



3.3 The vDVZ discontinuity 

We would now like to extract some physical predictions from the point source solution. Let's 
assume we have a test particle moving in this field, and that this test particle responds to 
h^j_u in the same way that a test particle in general relativity responds to the metric deviation 
^QijLu = -^h^^. We know from the textbooks (see for example chapter 7 of [36]) that if h^^, 
takes the form 2hoo/Mp = —2(f), 2hij/Mp = —2ilj6ij, h^i = for some functions 0(r) and 
ip{r), then the newtonian potential experienced by the particle is given by 0(r). Furthermore, 
if ■ip{r) = 'j(f){r) for some constant 7, called the PPN parameter, and if (/)(r) = — ^ for some 
constant k, then the angle for the bending of light at impact parameter b around the heavy 
source is given by a = 2(1 + 7)/^. Looking at (3.20), the massless graviton gives us the 



28 



values 



GM GM 

ip = , massless graviton, (3.23) 



using = SnG. The PPN parameter is therefore 7 = 1 and the magnitude of the hght 
bending angle for light incident at impact parameter b is 

AGM , , 

a = — ; — , massless graviton. (3.24) 
b 

For the massive case, the metric (3.13) is not quite in the right form to read off the 
newtonian potential and light bending. To simplify things, we notice that while the massive 
gravity action is not gauge invariant, we have assumed that the coupling to the test particle 
is that of GR, so this coupling is gauge invariant. Thus we are free to make a gauge trans- 
formation on the solution h^,^, and there will be no effect on the test particle. To simplify 
the metric (3.13), we go back to (3.11) and notice that the term in hij is pure gauge, so 
we can ignore this term. Thus our metric is gauge equivalent to the metric 

, , , 2M 1 e"™^ 

SMpA-K r 
hoi{x) = 0, 

M 1 e"'"'' 

hiAx) = Sii. (3.25) 

We then have, in the small mass limit, 

AGM 2GM ^ 

= , ip = Oij. massive graviton, (3.26) 

3 r 3 r 

These are the same values as obtained for u = Brans-Dicke theory. The newtonian potential 
is larger than for the massless case. The PPN parameter is 7 = ^, and the magnitude of the 
light bending angle for light incident at impact parameter b is the same as in the massless 
case, 

a = — - — , massive graviton. (3.27) 

If we like, we can make the newtonian potential agree with GR by scaling G — > ^G. Then 
the light bending would then change to a = ^^^y^, off by 25 percent from GR. 

What this all means is that linearized massive gravity, even in the limit of zero mass, 
gives predictions which are order one different from linearized GR. If nature were described 
by either one or the other of these theories, we would, by making a finite measurement. 
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be able to tell whether the graviton mass is mathematically zero or not, in violation of our 
intuition that the physics of nature should be continuous in its parameters. This is the vDVZ 
(van Dam, Veltman, Zakharov) discontinuity [29, 30] (see also [43, 44]). It is present in other 
physical predictions as well, such as the emission of gravitational radiation [45]. 

4 The Stiickelberg trick 

We have seen that there is a discontinuity in the physical predictions of linear massless 
gravity and the massless limit of linear massive gravity, known as the vDVZ discontinuity. 
In this section, we will expose the origin of this discontinuity. We will see explicitly that 
the correct massless limit of massive gravity is not massless gravity, but rather massless 
gravity plus extra degrees of freedom, as expected since the gauge symmetry which kills the 
extra degrees of freedom only appears when the mass is strictly zero. The extra degrees of 
freedom are a massless vector, and a massless scalar which couples to the trace of the energy 
momentum tensor. This extra scalar coupling is responsible for the vDVZ discontinuity. 

Taking m — > straight away in the lagrangian (3.1) does not yield not a smooth limit, 
because degrees of freedom are lost. To find the correct limit, the trick is to introduce new 
fields and gauge symmetries into the massive theory in a way that does not alter the theory. 
This is the Stiikelberg trick. Once this is done, a limit can be found in which no degrees of 
freedom are gained or lost. 

4.1 Vector example 

To introduce the idea, we consider a simpler case, the theory of a massive photon coupled 
to a (not necessarily conserved) source J^, 

S = Jd^'x - \f^uF^' - ^m'A^A^ + A^J^ (4.1) 

where F^y = d^j^Ay — dyA^. The mass term breaks the would-be gauge invariance, 5A^ = 9^A, 
and for = 4 this theory describes the 3 degrees of freedom of a massive spin 1 particle. 
Recall that the propagator for a massive vector is ^2^^ {Vt^v + which goes like ~ ^ 

for large momenta, invalidating the usual power counting arguments. 
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As it stands, the limit m — )■ of the lagrangian (4.1) is not a smooth hmit because 
we lose a degree of freedom - for m = we have Maxwell electromagnetism which in 
D = A propagates only 2 degrees of freedom, the two polarizations of a massless helicity 1 
particle. Also, the limit does not exist unless the source is conserved, as this is a consistency 
requirement in the massless case. 

The Stiickelberg trick consists of introducing a new scalar field 0, in such a way that 
the new action has gauge symmetry but is still dynamically equivalent to the original action. 
It will expose a different m — )■ limit which is smooth, in that no degrees of freedom are 
gained or lost. We introduce a field, 0, by making the replacement 

A^^A^ + d^<j), (4.2) 

following the pattern of the gauge symmetry we want to introduce [46]. This is emphati- 
cally not a change of field variables. It is not a decomposition of into transverse and 
longitudinal parts (A^ is not meant to identically satisfy c^^A'" = after the replacement), 
and it is not a gauge transformation (the lagrangian (4.1) isn't gauge invariant). Rather, 
this is creating a new lagrangian from the old one, by the addition of a new field 0. F^^^, 
is invariant under this replacement, since the replacement looks like a gauge transformation 
and Ffj_^ is gauge invariant. The only thing that changes is the mass term and the coupling 
to the source, 

S = jd^x - \f^.F^^ - \m\A^ + 5^0)2 + A^J^ - <pd^jr (4.3) 

We have integrated by parts in the coupling to the source. The new action now has the 
gauge symmetry 

6A^ = d^A, 6(t> = -A. (4.4) 

By fixing the gauge = 0, called the unitary gauge (a gauge condition for which it is 
permissible to substitute back into the action, because the potentially lost equation is 
implied by the divergence of the A^ equation) we recover the original massive lagrangian 
(4.1), which means (4.3) and (4.1) are equivalent theories. They both describe the three 
degrees of freedom of a massive spin 1 in D = 4. The new lagrangian (4.3) does the job 
using more fields and gauge symmetry. 

The Stiikelberg trick is a terrific illustration of the fact that gauge symmetry is a 
complete sham. It represents nothing more than a redundancy of description. We can take 
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any theory and make it a gauge theory by introducing redundant variables. Conversely, 
given any gauge theory, we can always eliminate the gauge symmetry by eliminating the 
redundant degrees of freedom. The catch is that removing the redundancy is not always a 
smart thing to do. For example, in Maxwell electromagnetism it is impossible to remove 
the redundancy and at the same time preserve manifest Lorentz invariance and locality. Of 
course, electromagnetism with gauge redundancy removed is still electromagnetism, so it 
is still Lorentz invariant and local, it is just not manifestly so. With the Stiikelberg trick 
presented here, on the other hand, we are adding and removing extra gauge symmetry in a 
rather simple way, which does not mess with the manifest Lorentz invariance and locality. 

We see from (4.3) that cf) has a kinetic term, in addition to cross terms. Rescaling 
(h — )■ —d) in order to normalize the kinetic term, we have 

S = jd^'x - \f^uF^^ - \n^A^A^ - mA^d^(t> - ]^d^ct>d^(t> + - (4-5) 
and the gauge symmetry reads 

5A^ = d^K, 6(f) = -mA. (4.6) 

Consider now the m — )■ limit. Note that if the current is not conserved (or its 
divergence does not go to zero with at least a power of m [47]), then the scalar becomes 
strongly coupled to the divergence of the source and the limit does not exist. Assuming a 
conserved source, the lagrangian becomes in the limit 



£ = -^F.^F^-' - ^9^09^0 + A^J^ (4.7) 



and the gauge symmetry is 

SA^ = d^A, 6(t> = 0. (4.8) 

It is now clear that the number of degrees of freedom is preserved in the limit. For D = 4 
two of the three degrees of freedom go into the massless vector, and one goes into the scalar. 

In the limit, the vector decouples from the scalar, and we are left with a massless gauge 
vector interacting with the source, as well as a completely decoupled free scalar. This m — t- 
limit is a different limit than the non-smooth limit we would have obtained by taking m — 
straight away in (4.1). We have scaled — )■ ^0 in order to canonically normalize the scalar 
kinetic term, so we are actually using a new scalar (pnew = m(j)oid which does not scale with 
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m, so the smooth hmit we are taking is to scale the old scalar degree of freedom up as we 
scale m down, in such a way that the new scalar degree of freedom remains preserved. 

Rather than unitary gauge, we can instead fix a Lorentz-like gauge for the action (4.3), 

df.A^' + m0 = 0. (4.9) 

This gauge fixes the gauge freedom up to a residual gauge parameter satisfying (□ — m^)A = 
0. We can add the gauge fixing term 

-^{d^A" + m4>f . (4.10) 

As in the massless case, quantum mechanically this term results from the Fadeev-Popov 
gauge fixing procedure. Classically, we may justify it on the grounds that the equations of 
motion obtained from the action plus the gauge fixing term are the same as the gauge fixed 
equations of motion (the gauge condition itself, however, is not obtained as an equation of 
motion, and must be imposed separately). Adding the gauge fixing term diagonalized the 
lagrangian, 

S + Sgf= [ d^'x Ia,{D - m')A^ + ^</>(n - m^)4> + A^J^ - -(^9^ J^ (4.11) 
and the propagators for Ay, and are respectively 



2^ 2' 2^ 2' (4-12) 

pZ _|_ pZ _|_ jy^Z 

which go like ~ ^ at high momenta. Thus we have managed to restore the good high energy 
behavior of the propagators. 

It is possible to find the gauge invariant mode functions for and (f), which can then 
be compared to the unitary gauge mode functions of the massive photon. In the massless 
limit, there is a direct correspondence; (p is gauge invariant and becomes the longitudinal 
photon, the A^ has the usual Maxwell gauge symmetry and its gauge invariant transverse 
modes are exactly the transverse modes of the massive photon. 

4.2 Graviton Stiikelberg and origin of the vDVZ discontinuity 

Now consider massive gravity, 

S = jd^'x Cm=o - Im^hy^h^" - h'') + ^hyuT^", (4.13) 
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where £m=o is the lagrangian of the massless graviton. We want to preserve the gauge 
symmetry 5h^y = d^^y + d^^^ present in the m = case, so we introduce a Stiickelberg field 
patterned after the gauge symmetry, 

h^^ ^h^^ + d^,A^ + d^A^. (4.14) 

The £m=o term remains invariant because it is gauge invariant and (4. 14) looks like a gauge 
transformation, so all that changes is the mass term, 

S = j d^x Cm=o - \m\h^,h>^'' - h^) - ^m^F^^F'^-' - 2m' {h^^d^A'' - hd^A^) 

+Kh^,T^'' - 2KA^d,T^'', (4. 15) 

where we have integrated by parts in the last term, and where F^j_y = dfj,Ay — duA^. 
There is now a gauge symmetry 

5K. = d^i, + d,i^, 5A^ = -i^, (4.16) 

and fixing the gauge = recovers the original massive gravity action (as in the vector 
case, this is a gauge condition for which it is permissible to substitute back into the action, 
because the potentially lost A^ equation is implied by the divergence of the h^u equation). 
At this point, we might consider scaling m^/^ normalize the vector kinetic term, 

then take the m — )■ limit. In this limit, we would end up with a massless graviton and 
a massless photon, for a total of 4 degrees of freedom (in 4 dimensions). So at this point, 
m — )■ is still not a smooth limit, since we would be losing one of the original 5 degrees of 
freedom. 

We have to go one step further and introduce a scalar gauge symmetry, by introducing 
another Stiickelberg field 0, 

A^^A, + d,(l). (4.17) 

The action (4.15) now becomes 

2 ^ ^ ; 2 

- 2m' {h^^d^'A" - hd^A^") - 2m' ( ^^^9^0 - hd'<i)) + nh^^T^"^ 

- 2KA^d^T^"' + 2K(f)ddT, (4.18) 
where ddT = dndpT'^^ and we have integrated by parts in the last term. 



S= I d^'x Cm=o - ^m'{h^,h^^ - h') - -m'F.^F'^'' 
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There are now two gauge symmetries 



5V = d^i, + d,i^, 5A^ = -^^ (4.19) 
5A^ = d^A, 50 = -A. (4.20) 

By fixing the gauge = we recover the lagrangian (4.15) . 

Suppose we now rescale — )■ :^A^, — )■ ^0, under which the action becomes 



- 2m {h^^d^A'' - hd^A") - 2 - hd^cp) + Kh^^T'"' 

- —KA.d^T^"' + ^K<i)ddT, (4.21) 

and the gauge transformations become 

SA^ = d^A, (50 = -mA, (4.22) 

where we have absorbed one factor on m into the gauge parameter A. 

Now take the m — > hmit. (If the source is not conserved and the divergences do not 
go to zero fast enough with m [47] , then and become strongly coupled to the divergence 
of the source, so we now assume the source is conserved.) In this limit, the theory now takes 
the form 

S = jd^'x Cm=o - ^F^uF^" - 2 ( V^^^'</> - hd^<P) + ^h^,T^\ (4.23) 

we will see that this has all 5 degrees of freedom; a scalar tensor vector theory where the 
vector is completely decoupled but the scalar is kinetically mixed with the tensor. 

To see this, we will un-mix the scalar and tensor, at the expense of the minimal coupling 
to T^"^, by a field redefinition. Consider the change 

V = h'^, + Trr^^,, (4.24) 

where tt is any scalar. This is the linearization of a conformal transformation. The change 
in the massless spin 2 part is (no integration by parts here) 



Cm=o{h) = Cm=o{h') + {D-2) 



(4.25) 
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This is simply the hnearization of the effect of a conformal transformation on the Einstein- 
Hilbert action. 

By taking vr = in the transformation (4.24), we can arrange to cancel all the 

off-diagonal h(j) terms in the lagrangian (4.23), trading them in for a (p kinetic term. The 
lagrangian (4.23) now takes the form, 

S = Jd^'x C^=o{h') - \f,.F^'' - 2^^a^0a'^0 + kH'^^T^^^ + (4.26) 
and the gauge transformations read 

= d^i, + d,i,, M^ = (4.27) 

5A^ = d^A, 50 = 0. (4.28) 

There are now (for D = 4) manifestly five degrees of freedom, two in a canonical massless 
graviton, two in a canonical massless vector, and one in a canonical massless scalar''. 

Note however, that the coupling of the scalar to the trace of the stress tensor survives 
the 171 = limit. We have exposed the origin of the vDVZ discontinuity. The extra scalar 
degree of freedom, since it couples to the trace of the stress tensor, does not affect the 
bending of light (for which T = 0), but it does affect the newtonian potential. This extra 
scalar potential exactly accounts for the discrepancy between the massless limit of massive 
gravity and massless gravity. 

As a side note, one can see from this Stiickelberg trick that violating the Fierz-Pauli 
tuning for the mass term leads to a ghost. Any deviation from this form, and the Stiickelberg 
scalar will acquire a kinetic term with four derivatives ~ (D^)^, indicating that it carries 
two degrees of freedom, one of which is ghostlike [48, 49]. The Fierz-Pauli tuning is required 
to exactly cancel these terms, up to total derivative. 

Returning to the action for m ^ (and a not necessarily conserved source), we now 
know to apply the transformation hf^^ = /i'^^ + j^4>rifj_^, which yields, 

^Ordinarily the Maxwell term should come with a 1/4 and the scalar kinetic term with a 1/2, but we 
leave different factors here just to avoid unwieldiness. 
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+ Kh'T^^'' + -^K<j)T - -KA^d.T'^'' + ^K(l)ddT. (4.29) 
^ D — 2 m 

The gauge symmetry reads 

2 

D~2 

6A^ = df,A, 6(f) = -mA. (4.31) 



^K'^ = d^i, + d,i,, + ——mkr^^,, 5A^ = -mi^ (4.30) 



We can go to a Lorentz-like gauge, by imposing the gauge conditions [50, 51] 

d''h'^,-]^d,h' + mA, = Q, (4.32) 

d^A>^ + m {hi' + 2^^0^ = 0. (4.33) 

The first condition fixes the symmetry up to a residual transformation satisfying (□ — 
m^)^^ = 0. It is invariant under A transformations, so it fixes none of this symmetry. The 
second condition fixes the A symmetry up to a residual transformation satisfying (□— m^)A = 
0. It is invariant under transformations, so it fixes none of this symmetry. We add two 
corresponding gauge fixing terms to the action, resulting from either Fadeev-Popov gauge 
fixing or classical gauge fixing, 

^GFi = j d''^ - (d'h'^, - ^d,h' + mA^^ , (4.34) 

5gf2 = jd''^ - (^M^^ + ^ (^^' + 2§^</') ) = 0- (4-35) 
These have the effect of diagonalizing the action, 

^ + ^GFi + ^GF2 = j d''xhi'^,{n-m')h'^'' -^-h' {n-m^)h' 

+A^ (□ - m^) A^^ + 2^^0 (□ - m") 

+Kh'Tf''' + y^K<pT - -KA^d^T'^" + ^n4>ddT. 
^ D — 2 m 

(4.36) 
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The propagators of h'^^^, and are now, respectively, 

1 -17]^^ D 



p2 _j_ ^2 



1 / N 1 



(4.37) 



which all behave as ~ ^ for high momenta, so we may now apply standard power counting 
arguments. 

With some amount of work, it is possible to find the gauge invariant mode functions 
for h'^^, A^ and 0, which can then be compared to the unitary gauge mode functions of 
Section 2.2. In the massless limit, there is a direct correspondence; (f) is gauge invariant 
and its one degree of freedom is exactly the longitudinal mode (2.20), the A^ has the usual 
Maxwell gauge symmetry and its gauge invariant transverse modes are exactly the vector 
modes (2.21), and finally the /i^^ has the usual massless gravity gauge symmetry and its 
gauge invariant transverse modes are exactly the transverse modes of the massive graviton. 



4.3 Mass terms as filters and degravitation 

There is a way of interpreting the graviton mass as a kind of high pass filter, through which 
sources must pass before the graviton sees them. For a short wavelength source, the mass 
term does not have much effect, but for a long wavelength source (such as the cosmological 
constant), the mass term acts to screen it, potentially explaining how the observed cosmic 
acceleration could be small despite a large underlying cosmological constant [52]. 

First we will see how this works in the case of the massive vector. Return to the action 
(4.5), with a conserved source, before taking the m — )■ limit, 

S = j d^'x - ^F^uF^" - ^m'A^A^ - mA^d^ - ^d^c^d'^cP + A^r. (4.38) 
The equation of motion is 

□0 + m9-A = O. (4.39) 

We would now like to integrate out cj). Quantum mechanically we would integrate it out of 
the path integral. Classically we would eliminate it with its own equation of motion. Solving 
the equation of motion involves solving a differential equation, so the result is non-local, 

777 

= --5-A (4.40) 
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Plugging back into (4.38), we obtain a non-local lagrangian 

S = jd^'x - \f,, [l - ^) F^'^ + J^ (4.41) 

where we have used F^^^F^'^ = —2A^^A'^ — 2d ■ A^d ■ A, arrived at after integration by 
parts. The lagrangian (4.41) is now a manifestly gauge invariant but non-local lagrangian for 
a massive vector. The non-locality results from having integrated out the dynamical scalar 
mode. The equation of motion from (4.41) is 

— jd,F^'' = -r. (4.42) 

This is simply Maxwell electromagnetism, where the source is seen thorough a filter ^1 — 
For high momenta p ^ m, the filter is ~ 1 so the theory looks like ordinary electromag- 
netism. But for p <^ m, the filter becomes very small, so the source appears weakened. We 
can think of this as a high-pass filter, where m is the filter scale. 

Applied to gravity, the hope is to explain the small observed value of the cosmological 
constant. The cosmological constant, being a constant, is essentially a very long wavelength 
source. Gravity equipped with a high pass filter would not respond to a large bare cosmolog- 
ical constant, making the observed effective value appear much smaller, while leaving smaller 
wavelength sources unsuppressed. This mechanism is known as degravitation [53, 54, 52, 55]. 

This filtering is essentially just the Yukawa suppression e"™'^ that comes in with massive 
particles, so we should be able to cast the massive graviton into a filtered form. Look again 
at the action (4.15) with a conserved source, before introducing the Stiikelberg scalar, 

S = jd^'x Cm=o - \m\h^,W'' - e) - hn'F^.F^" - 2m' {h^^d'^A" - hd^A'^) + /^V^^'- 

(4.43) 

Now consider the following action containing an additional scalar field N, 

I, .„., 1 



S = I d^x Cm=o + rn? 



^Kuh^" + -^h' + A^UA^' + N{h - N) 



- A^{d^h-2d''h^, + 2d^N) 
The field is an auxiliary field. Its equation of motion is 



+ >^KuT^r (4.44) 



N=hi + (4.45) 
39 



which when plugged into (4.44) yields (4.43). Thus the two actions are equivalent, and 
(4.44) is another action describing the massive graviton. Here, however, there is no gauge 
symmetry acting on the scalar; is gauge invariant^. 

Instead of eliminating the scalar, we can eliminate the vector using its equations of 
motion, 

= ^ Qs^/i - 9^ V + d,N^ ■ (4.47) 

Plugging back into (4.44) gives 

S = Jd^'x ^ V (l - ^) ^^"'"^/ia/? - 2N^ {d^d.h^^ - Uh) + /^VT^^ (4.48) 

where E^^^j^ is the second order differential operator for the massless graviton (2.46). Now, 
to diagonalize the action, make a conformal transformation 

after which (4.48) becomes 

(4.50) 

Finally, making the field redefinition N' = ^^^^N to render the coupling to the source local, 

S = Jd^x \h'^^ [l - ^) S'^'^^-^K, + 2§^iV'(n - m')N' + Kh'^^T^'^ + ^^iV'T. 

(4.51) 

Thus a massive graviton is equivalent to a filtered graviton coupled to T^^, and a scalar with 
mass m coupled with gravitational strength to the trace T. The scalar is the longitudinal 
mode responsible for the vDVZ discontinuity. 

It is not hard to see that a linear massive graviton screens a constant source. Looking 
at the equations of motion (3.2) where the source is a cosmological constant T^y oc rj^y, 



^For another form of the massive gravity action, we can take N' ^ N — 9^v4^ in (4.15), which gives 
S = d^x £,,1=0 + '71^ 



l^.h^" + ^h^ - ^F^.F'^'' + N'ih - N') - 2A'^ (d^h - d'^h^y) 



+ nh^yT^"" ~ 2KA^dyT^'r (4.46) 
The field N' now takes the value N' — and is no longer gauge invariant. 
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and taking the double divergence, we find d^d^h^j^i, — D/i = 0, which is the statement that 
the hnearized Ricci scalar vanishes, so a cosmological constant produces no curvature. If 
degravitation can be made to work cosmologically, then this provides an interesting take 
on the cosmological constant problem. Of course the smallness of the cosmological constant 
reappears in the ratio m/Mp, but as we will see, in massive gravity a small mass is technically 
natural. There are other obstacles as well, and promising avenues towards overcoming them, 
and we will have more to say about these things while studying the non-linear theory. 



5 Massive gravitons on curved spaces 

We now study some new features that emerge when the Fierz-Pauli action is put onto a 
curved space. One new feature is the existence of partially massless theories. These are 
theories with a scalar gauge symmetry that propagate 4 degrees of freedom in D — A. 
Another is the absence of the vDVZ discontinuity in curved space. 



5.1 Fierz-Pauli gravitons on curved space and partially massless 
theories 



We now study the linear action for a massive graviton propagating on a fixed curved back- 
ground with metric g^„. As in the flat space case, the massless part of the action will be the 
Einstein-Hilbert action with a cosmological constant, ^y/^{R — 2A), expanded to second 
order in the metric perturbation Sg^i, — 2Ac/i^j,, about a solution g^^,^. The solution must be 
an Einstein space, satisfying 



R 



R 



9nu, A 



D-2 



D"'^ ' V 2D 

Appending the Fierz-Pauli mass term, we have the action 



R. 



(5.i; 



(5.2) 



Here the metric, covariant derivatives and constant curvature R are those of the background. 
Notice the term, proportional to R, that kind of looks like a mass term, but does not have 
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the Fierz-Pauli tuning. There's some representation theory behind this [56], and a long 
discussion about what it means for a particle to be "massless" in a curved space time [57], 
but at the end of the day, (5.2) is the desired generalization of the flat space Fierz-Pauli 
action, which, for most choices of m^, propagates 5 degrees of freedom in D = 4. See 
[58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72] for some other aspects of massive 
gravity on curved space. 

For some choices of m^, (5.2) propagates fewer degrees of freedom. For m = 0, the 
action has the gauge symmetry 

Sh^, = V^C, + V.^M, (5.3) 

and the action propagates 2 degrees of freedom in D = A. As we will see momentarily, for 
R = '^^"'^'' m^, m 7^ 0, the action has a scalar gauge symmetry, and propagates 4 degrees 
of freedom in D = 4. For all other values of and R, it has no gauge symmetry and 
propagates 5 degrees of freedom in D = 4. This is summarized in Figure 1. 

We introduce a Stiickelberg field, A^, patterned after the m = gauge symmetry, 

V ^ V + + V^A^. (5.4) 

The Cm=o term remains invariant, the source term does not change because we will assume 
covariant conservation of T'^'^, so all that changes is the mass term. 



S = I d^x Cm=0 + 



- ^m'ih.^h'^'' - h^) 



- ^n'F^^F^" + '^m'RA^A^ - 2m^ ( VV^A'^ - hV ^A^') + t^h^^T^" 

(5.5) 

where F^^ = d^Ay — dyA^ = V^A^ — V^A^, and we have used the relation V f^A^V^A^^ = 
(V^A^)^ — R^yA^A'^ to see that there is now a term that looks like a mass for the vector, 
proportional to the background curvature. There is now a gauge symmetry 

5V = V^^- + V.e^' M^ = -e^, (5.6) 
and fixing the gauge = recovers the original action (5.2). 

Introducing the Stiikelberg scalar and its associated gauge symmetry, 

A^^A^ + V^(t), M^ = V/,A, (50 = -A, (5.7) 
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we have 



S 



1 



— -m 



+ 



2m^R 



D 



(5.8) 



Under the conformal transformation 



V = h'txv + T^9nv, (5.9) 
where tt is any scalar, the change in the massless part is (no integration by parts here) 



-9 



D 



{D - 2) (v^irWh' - V^ttV./i'^'^ + ~ l)V^7rV'^7r^ 

(5.10) 



Applying this in the case tt = ■j^m?4' yields, 



-m 



D 



2m^ {h'^^V^A'' - h'Vf.A'') + 2w? 



- 2m' 



D - 1 

-i 

D-2 

D-2 D V D-2 



-m 



R 
D 



{2(f)V|,A^' + h'(l)) 



D - 1 



R 



+ Kh'^^T^^ + ^—^m^Kct^T 



(5.11) 



The gauge symmetry reads 



5h'. 



d.A, 6(P = -A. 



(5.12) 

(5.13) 



Note that for the special value 



R = 

£1-2 
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(5.14) 



the dependence on completely cancels out of (5.11). Setting unitary gauge = 0, and 
given the replacements (5.4), (5.7) and the conformal transformation (5.9), this implies that 
the original lagrangian (5.2) with the mass (5.14) has the gauge symmetry 

Shf,^ = V^VuX + jj^^m^Xg/Mu, (5.15) 

where X{x) is a scalar gauge parameter. The theory at the value (5.14) is called partially 
massless [73, 74, 75, 76, 77]. Due to the gauge symmetry (5.15), this theory propagates one 
fewer degree of freedom than usual so for D = A it carries four degrees of freedom rather 
than five. Consistency demands that the trace of the stress tensor vanish for this theory (if 
it is conserved). In addition, it marks a boundary in the R,ni? plane between stable and 
unstable theories, see Figure 1. 



R 









R=6m^ 


2 DOF (stable) ^ 


^ I^OF (stable) 






5 DOF (unstable) \\ 




\\X W 


^ X- m 




5 DOF (stable) 


4 DOF (unstable) 





Figure 1: Degrees of freedom and their stability for values in the R,m? plane for massive gravity 
on an Einstein space (shown for L) = 4, the other dimensions follow similarly). The line R = 
6m?, 7^ is where a scalar gauge symmetry appears, reducing the number of degrees of freedom 
by one. The line = is where the vector gauge symmetries appear, reducing the number of 
degrees of freedom by three. 
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5.2 Absence of the vDVZ discontinuity on curved space 



To study the fate of the vDVZ discontinuity, we now take a massless hmit, while preserving 
the number of degrees of freedom. However, there are many paths in the i?, plane and 
correspondingly, different ways to take the massless limit. 

For example, let's take the m — )■ limit for fixed but non-zero R. Here we will see that 
the vDVZ discontinuity is absent [73, 78, 79, 80]. First we go to canonical normalization 
for the vector by taking — )■ Then we notice that we can immediately take the 

m — )■ limit, without the need to introduce the second Stiickelberg field 0. This is because 
a mass term for the vector is present in this limit, so no degrees of freedom are lost. Thus 
our limiting action is 



2 + ^ 



(5.16) 



The massive vector completely decouples from the stress tensor, so there is no vDVZ dis- 
continuity. Notice that the vector is a tachyon in dS space but healthy in AdS, consistent 
with the stability regions shown in Figure 1. These regions can all be investigated in similar 
fashion. Finally, note also that the i? — )■ and m — limits do not commute. 



6 Non-linear interactions 

Up to this point, we have only studied the linear theory of massive gravity, which is deter- 
mined by the requirement that it propagate only one massive spin 2 degree of freedom. We 
now turn to the study of the possible interactions and non-linearities for massive gravity. 

6.1 General relativity 

We start by reviewing the story of non-linearities in GR. We will then repeat it for massive 
gravity, to see exactly where things differ. General relativity is the theory of a dynamical 
metric g^;,y^ with the action 

S = ^ [d^'x ^R. (6.1) 
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The action is invariant under (pullback) diffeomorphism gauge symmetries f^ix), 

9f^-(^) ^ d^d^^''^ ^^^^^^ ■ ^^-^^ 
Infinitesimally, for f^{x) = + ^'^{x), this reads 

^Qfiu = ^^gf^u = ^t,^u + ^ui^,, (6.3) 

where is the gauge parameter, is the Lie derivative, and indices are lowered by the 
metric. 

The field equation for the metric is 

- \r9^.u = 0, (6.4) 

and the most symmetric solution is fiat space g^j.^ = rj^i,. 

To see that this is a theory of an interacting massless spin 2 field, we expand the action 
around the fiat space solution r]^^, 

To second order in hf^^ the action is 

(6.5) 

where indices on /i^,^ are raised and traced with the fiat background metric rj^y and we have 
ignored total derivatives. After canonical normalization, h^y = 2nh^y, this linear action for 
GR is exactly that of the m = spin two particle in Minkowski space (2.1). 

If we continue the expansion around fiat space to higher non-linear order in h^y, we 
get a slew of interaction terms, all with two derivatives and increasing powers of /i, and 
coefficients all precisely fixed so that the result is diffeomorphism invariant and sums up to 
(6.1). Schematically, 

S = j d^X d^h^ + Kd^h"^ + ■■■ + + . . . (Q_Q) 

The higher and higher powers of h^y are suppressed by appropriate powers of k. The action 
is expanded in powers of nh and the linearized expansion is valid when nh <til. 
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When expanding around the background, we must put all of the gauge symmetry into 
hfj,iy, SO that the transformation rule is 

K.{x) ^ d^rd,f7]^p - 7]^, + d^rd.fKp (fix)) . (6.7) 

For infinitesimal transformations = + S,^, this is 

5hf,i, = d^,iu + dyi^ + C^h^y. (6.8) 

This is the all orders expression in h^i, for the infinitesimal gauge symmetry, which shows 
that it gets modified at higher order from its linear form (2.2). 

This argument can be turned around. We can start with the massless linear graviton 
action and ask what higher power interaction terms can be added. The possible terms can be 
arranged in powers of the derivatives, and lower derivatives will be more important at lower 
energies. Starting with two derivatives, we ask what terms of the form (6.6) can be added. 
We must add higher order terms in such a way that the linear gauge invariance is preserved, 
though the form of the gauge transformations may be altered at higher order in h. These 
requirements are strong enough to force the interactions to be those obtained from expanding 
the Einstein-Hilbert term [4, 5, 6, 7, 8, 9, 10]. Looking at it from this direction, an amazing 
thing happens. The linear action with which we started has a non-dynamical background 
metric 7]^^, but after adding all the interactions of the Einstein-Hilbert term, the change of 
variables h^^ — > g^^i, — rj^y completely eliminates the background metric from the action. The 
fully interacting Einstein-Hilbert action turns out to be background independent. This will 
not be the case once we add a mass term, though we will still be able to introduce gauge 
invariance through the Stiikelberg trick. 



Zero derivative interactions mean curved space 

If we ask for interactions terms with fewer than two derivatives, the only option is zero 
derivatives, and diffeomorphism invariance forces them to sum up to a cosmological constant 
^—g = + O (/i^). This contains a term linear in h, which means there is a tadpole and 
h = is not a solution to the equations of motion, so we are not expanded around a solution 
[81]. Instead we may consider GR with a cosmological constant A, 

I rf^x v^(i?-2A). (6.9) 



47 



The equations of motion G^i^ + h.g^p = 0, implies that the background solution gj^J is an 
Einstein space, 

R,. = ^9,., ^=(^^)^- (6.10) 

Expanding around the background gfj_u = gj^J + hf^,, to quadratic order, we have the linearized 
action 



(fx \f--g 



1 

D \ ^ 2 



f (total d), 



where the covariant derivatives, curvature and metric determinant out front are those of the 
background. After canonical normalization, h^^ = 2K,h^y, this linear action is exactly that 
of the m = spin 2 particle in an Einstein space (5.2) we used in Section 5. 

Upon expanding around the background, we must put all of the gauge transformations 
into /i^i/, so that the transformation rule is 

V(^) ^ d.rdj^g^^l (/(x)) - g^^J + d.rd.fKp (fix)) . (6.11) 

For infinitesimal transformations = x^ + S,^, this is 

(^V = ^^i- + V.^M + A V' (6-12) 

where the covariant derivatives are of the background. This is an all orders expression in h^i, 
for the infinitesimal gauge symmetry. To linear order, this reproduces the massless curved 
space gauge symmetry (5.3). As in the fiat space case, this argument may be reversed. 
The only possible interactions for a massless graviton propagating on an Einstein space (an 
Einstein space is the only space on which a free graviton can consistently propagate [70]) 
should be those of Einstein gravity with a cosmological constant. 



Spherical solutions 

Returning now to a fiat background A = 0, and setting D = 4, we attempt to find spherically 
symmetric static solutions to the equations of motion R^j^y — \Rgfj.y = 0, using an expansion 
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in powers of non-linearity, a method we will repeat for the non-linear massive graviton. The 
most general spherically symmetric static metric can be written 

Qi^^dx^dx" = -B{r)dt^ + C{r)dr'^ + A{r)r'^dn'^. (6.13) 

The most general gauge transformation which preserves this ansatz is a reparametrization of 
the radial coordinate r. We can use this to set the gauge A{r) = C{r), bringing the metric 
into the form 

Qf.^dxi'dx'' = -B{r)de + C(r) [dr^ + r^d^'^] . (6.14) 

The linear expansion of this around flat space will be seen to correspond to the Lorenz gauge 
choice. Plugging this ansatz into the equations of motion, we get the following from the tt 
equation and rr equation respectively, 

3r (C)^ - 4C (2C" + rC") = 0, (6.15) 
AB'C^ + 2{2B + rB')C'C + Br{C'f = 0. (6.16) 

The 66 equation, (which is the same as the (jxj) equation by spherical symmetry) turns out to 
be redundant. It is implied by the tt and rr equations (this happens because of a Noether 
identity resulting from the radial re-parametrization gauge invariance). 

We start by doing a linear expansion of these equations around the flat space solution 

Bo{r) = 1, Co(r) = 1. (6.17) 

We do this by the method of linearizing a non-linear differential equation about a solution. 
We introduce the expansion 

B{r) = So(r) + eSi(r) + e'S2(r) + --- , (6.18) 

C{r) = Co{r) + eC,{r) + e'C2{r) + ■ ■ ■ , 

where e will be a parameter that counts the order of non-linearity. We proceed by plugging 
into the equations of motion and collecting like powers of e. The O{0) part gives = 
because Bo,Cq,Ao are solutions to the full non-linear equations. At each higher order in e 
we will obtain a linear equation that lets us solve for the next term in terms of the solutions 
to previous terms. 
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At 0{e) we obtain 



2C 

C:' + -^ = 0, B[ + C[ = 0. (6.19) 

r 



There are three arbitrary constants in the general solution. Demanding that Bi and Ci go 
to zero as r — > oo, so that the solution is asymptotically flat, fixes two. The other constant 
remains unfixed, and represents the mass of the solution^. We choose it to reproduce the 
solution (3.22) we got from the propagator. We have then, 

c. = ?^. (6.21) 

At O(e^) we obtain another set of differential equations 

3G^M^ 2C' 



^2 ^// 



C'i = (6.22) 



7/^2 j\/r2 

+ B'^ + C'^ = 0. (6.23) 



Again there are three arbitrary constants in the general solution. Demanding that B2 and 
C2 go to zero as r — )■ 00 again fixes two. The third appears as the coefficient of a ^ term, 
and we set it to zero so that the second order term does not compete with the first order 



^If we had used the gauge A = 1, these would be 1st order equations and there would be only two 
constants to fix, i.e. the order of the equation seems to depend on the gauge. The reason for this is that the 
Lorenz gauge does not completely fix the gauge. Under an active diffeomorphism f{x), the metric transforms 
as gfjiv{x) -> dfif^duf^ g\a{f{x)). Under a radial reparametrization f(r) this becomes 

B{r)^B{f{r)), C{r) ^C{f{r)), A(r) ^ A(f(r))^. (6.20) 

or r 

Choosing the gauge A = 1 amounts to solving A(f(r))^^ = 1 for f, and this is an algebraic equation so the 
solution is unique and this choice completely fixes the gauge. Choosing the gauge A = C amounts to solving 
fC{r{r)) = A(f(r))^, which is a differential equation for r. Thus the solution is not unique because there 
is an integration constant. The transformations that preserve the gauge choice solve ^ = with solution 
f — kr, i.e. constant scalings of r. This appears as an extra boundary condition that must be specified, 
because we must fix the scaling. 

In addition, our ansatz is also invariant under time scaling, t = kt, under which B{r) — !■ kB{r). This 
represents another unfixed gauge symmetry. We generally fix this and the radial scaling by demanding that 
A, i? — >■ 1 as r — > 00. Then the only boundary condition is the mass. 
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term as r — oo. We can continue in this way to any order, and we obtain the expansion 

2GM f GM \ 

B{r)-l = 1 + ■■■ , 6.24 

r \ r J 

^, . ^ 2GM / 3GM \ 

C(r)-1 = ——h+ — + ...]. (6.25) 

The dots represent higher powers in the non-hnearity parameter e. We see that the non- 
hnearity expansion is an expansion in the parameter rg/r, where 

rs = 2GM, (6.26) 

is the Schwarzschild radius. Thus the Schwarzschild radius rs ~ M/Mp represents the radius 
at which non-hnearities become important. This scale can also be estimated straight from 
the lagrangian (6.6). The non-linear terms are suppressed relative to the linear terms by 
powers of the factor h/Mp. The linear solution is ~ so h/Mp becomes order one 

when r ~ M/Mp ~ rs- 

In GR, the linearity expansion can be easily summed to all orders by solving the original 
equations exactly, 

B{r) = ): T:2 . C{r) = (1 + ^)^ 



\^ ^ GM) 



This is the Schwarzschild solution, in Lorenz gauge. 



GR as a quantum effective field theory 

We can understand the previous results from an effective field theory viewpoint, and in the 
process check that the black hole solution is still valid despite quantum corrections. Pure 
Einstein gravity in D = 4 is not renormalizable. It contains couplings with negative mass 
dimension carrying the scale Mp. Thus it must be treated an effective field theory with 
cutoff at most Mp [82]. This can also be seen from scattering amplitudes; by dimensional 
analysis the 2 — )■ 2 graviton scattering amplitude at energy E goes like which becomes 
order one and violates unitarity at an energy E ~ Mp. 

Since we have an effective theory, we expect quantum mechanically the presence of a 
plethora of other operators in the effective action, suppressed by appropriate powers of Mp 
and order one coefficients. Higher derivatives term, those beyond two derivatives, will be 
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associated with higher order effects in powers of some energy scale over the cutoff. By gauge 
invariance, all operators with two derivatives sum up to y/—gR and correct the Planck 
mass, naively by order one. However, we can generate operators with higher numbers of 
derivatives^*^, suppressed by appropriate powers of the Planck scale, for example, 

' d^h\ -^d^h\ ^d^h\ ... (6.27) 



Ml 'Ml ' Ml 

By gauge invariance, they must sum up to higher order curvature scalars, multiplied by 
appropriate powers of Mp, for instance, schematically, 

^ —gRV'R ~ -^,d'h' + ^d'h^ + ■ ■ ■ , (6.29) 



My ^ Ml Ml 

^ ^9R' - ^d'h' + ^d'h^ + ■ ■ ■ (6.30) 



My Ml Ml 

These corrections include terms which are second order in the fields, but higher order in 
the derivatives, which naively lead to new degrees of freedom, some of which may be ghosts 
or tachyons. One might worry why these terms are generated here, however the masses 
of these ghosts and tachyons is always near or above the cutoff Mp, so they need not be 
considered part of the effective theory, since the unknown UV completion may cure them. In 
line with this logic, they must not be re-summed into the propagator (this would be stepping 
outside the Mp expansion), but rather treated as vertices in the effective theory. 

The important observation is that all these higher terms are suppressed relative to any 
term in the Einstein-Hilbert part by powers of the derivatives 

d I , , 

-w- ~ -v^- 6.31 

Mp Mpr ^ ' 

Thus, at distances r ^ more than a Planck length from the central singularity of our 
spherical solution, quantum effects are negligible. Only when approaching within a Planck 
length of the center does quantum gravity become important. The regimes of GR are shown 
in Figure 2. 

^''Notc that quantum corrections will also generate the terms with no derivatives, the cosmological constant. 
We can consistently declare that these are zero at the expense of a fine tuning. This is the usual cosmological 
constant problem [f9]. 
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An important fact about GR is that there exists this parametrically large middle regime 
in which the theory becomes non-linear and yet quantum effects are still small. This is the 
region inside the horizon r = rs but farther than a Planck length from the singularity. In this 
region, we can re-sum the linear expansion by solving the full classical Einstein equations, 
ignoring the higher derivative quantum corrections, and trust the results. This is the reason 
why we know what will happen inside a black hole, but we do not know what will happen 
near the singularity. As we will see, this crucial separation of scales, in which the scale of 
non-linearity is well separated from the quantum scale, does not always occur in massive 
gravity. It only occurs if the parameters of the interactions are tuned in a certain way. 



Quantum Classical 



Non-linear Linear 



J_ M 
Mp ^ mJ 



Figure 2: Regimes for GR. 



6.2 Massive general relativity 



We now turn to non-linearities in massive gravity. What we want in a full theory of massive 
gravity is some non-linear theory whose linear expansion around some background is the 
massive Fierz-Pauli theory 2.1. Unlike in GR, where the gauge invariance constrains the full 
theory to be Einstein gravity, the extension for massive gravity is not unique. In fact, there 
is no obvious symmetry to preserve, so any interaction terms whatsoever are allowed. 

The first extension we might consider would be to deform GR by simply adding the 
Fierz-Pauli term to the full non-linear GR action, that is, choosing the only non-linear 
interactions to be those of GR, 



S 



1 

2^ 



-gR) - v/r^^m2(^(°)^-^7^°^"^ ( - h^aKp) 



(6.32) 



Here there are several subtleties. Unlike GR, the lagrangian now explicitly depends on a 
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fixed metric g^^} , which we will call the absolute metric, on which the linear massive graviton 
propagates. We have h^^ = g^y — gj^J as before. The mass term is unchanged from its 
linear version, so the indices on h^^, are raised and traced with the absolute metric. The 
presence of this absolute metric in the mass term breaks the diffeomorphism invariance of 
the Einstein-Hilbert term. Note that there is no way to introduce a mass term using only 
the full metric (7^,^, since tracing it with itself just gives a constant, so the non-dynamical 
absolute metric is required to create the traces and contractions. 

Varying with respect to g^^, we obtain the equations of motion 

v^(i?^^ - iRgn + {g^'^^'^g^'^^'h^, - g^'^'^'K.g^^^^^) = 0. (6.33) 

Indices on R^^, are raised with the full metric, and those on h^^, with the absolute metric. We 
see that if the absolute metric gj^J satisfies the Einstein equations (6.4), then g/j^^ = gjjj , i.e. 
hfiu = 0, is a solution. When dealing with massive gravity and more complicated non-linear 
solutions thereof, there can be at times two different background structures. On the one 
hand, there is the absolute metric, the structure which breaks explicitly the diffeomorphism 
invariance. On the other hand, there is the background metric, which is a solution to the 
full non-linear equations, about which we may expand the action. Often, the solution metric 
we are expanding around will be the same as the absolute metric, but if we were expanding 
around a different solution, say a black hole, there would be two distinct structures, the 
black hole solution metric and the absolute metric. 

If we add matter to the theory and agree to use only minimal coupling to the metric 
gfj^u, then the absolute metric does not directly influence the matter. It is the geodesies and 
lengths as measured by the solution metric that we care about. Unlike in GR, if we have a 
solution metric, we cannot perform a diffeomorphism on it to obtain a second solution to the 
same theory. What we obtain instead is a solution to a different massive gravity theory, one 
whose absolute metric is related to the original absolute metric by the same diffeomorphism. 

Going to more general interactions beyond (6.32), our main interest will be in adding 
interactions terms with no derivatives, since these are most important at low energies. The 
most general such potential which reduces to Fierz-Pauli at quadratic order involves adding 
terms cubic and higher in h^^, in all possible ways 



{V^gR)~^nf\m'U{g'^'\h) 



(6.34) 
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where the interaction potential U is the most general one that reduces to Fierz-Pauli at 
linear order, 

Uig'^'l h) = U,{g^'\ h) + U,{g^''\ h) + U,{g^''\h) + U.ig^'^l h) + ■ ■ ■ , (6.35) 



Us{g^'\h) 



[h^]-[h]\ 

+Ci [h'] + C2 [h^] [h] + C3 [hf 



(6.36) 
(6.37) 



uM'\h) = +Dr[h^]+D2[e][h] + D,[h'Y + D,[e][hf + D,[h]\ (6.38) 
U^ig^'K h) = +Fi [h'] + F2 [h'] [h] + F3 [h'] [h]' + F4 [e] [h'] + F, [h'Ylh] 



The square bracket indicates a trace, with indices raised with i.e. [h] = 

[h'^] = g^^^'^^'ha/sg^^'^^'^h,^^, etc. The coefficients Ci, C2, etc. are generic coefficients. Note that 
the coefficients in Un{g^^\ h) for n > D are redundant by one, because there is a combination 
of the various contractions, the characteristic polynomial C™{h) (see Appendix A), which 
vanishes identically. Thus one of the coefficients in Un{g^'^\ h) ioi n > D (or any one linear 
combination) can be set to zero. 

If we like, we can re-organize the terms in the potential by raising and lowering with 
the full metric g^'^ rather than the absolute metric g^^'>^^'^ ^ 



+F,[h'] [hf + F,[h]\ 



(6.39) 



S 



1 

2^ 




(6.40) 



where 



V{g, h) = V2{g, h) + 1/3(^7, h) + V^{g, h) + V,{g,h) + ■ ■ ■ , 



(6.41) 



V2{g,h) 
V3{g,h) 
V,{g,h) 
V,{g,h) 



{h') - {h)\ 

+Cl(/^=^)+C2(/^')(/^)+C3(/^)^ 

+ d2{h'){h) + d^{hY + d,{h^){hy + d,{h)\ 
+fi{h') + f2{h'){h) + Mh'){hf + h{h'){h^) + h{h'')\h) 



(6.42) 
(6.43) 
(6.44) 



(6.45) 
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where the angled brackets are traces with the indices raised with respect to g^'^. It does not 
matter whether we use potential (6.34) with indices raised by the potential (6.40) 

with indices raised by g'^". The two carry the same information and we can easily relate 
the coefficients of the two by expanding the inverse full metric and the full determinant in 
powers of h^^ raised with the absolute metric, 

1. 1 / 1 



(6.47) 



-g = V-^7(°) 



2 4 V ^ 2 



The following is useful for this purpose, 



r) = E(-l)M ' )[^'^"]- (6-48) 



While the zero derivative interaction terms we have written in (6.34) are general, the 
two derivative terms are not, since we have demanded they sum up to Einstein-Hilbert. The 
potential has broken the diffeomorphism invariance, so there is no symmetry reason for the 
two derivative interaction terms to take the Einstein-Hilbert form. We could deviate from it 
if we like, but we will see later that there are good reasons why it is better not to. We may 
also conceivably add general interactions with more than two derivatives, but we will omit 
these for the same reasons we omit them in GR, because they are associated with higher 
order effective field theory effects which we hope will be small in suitable regimes. 



6.3 Spherical solutions and the Vainshtein radius 

We will now look at static spherical solutions, doing for massive gravity what we did for 
GR in Section 6.1. We specialize to four dimensions, and for definiteness we pick the action 
(6.32) with the minimal mass term. We attempt to find spherically symmetric solutions to 
the equations of motion 6.33, in the case where the absolute metric is flat Minkowski in 
spherical coordinates, 

g^Sdx^'dx" = -dt" + dr^ + r'^d^f. 
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We consider a spherically symmetric static ansatz for the dynamical metric^^ 

g^.^dx^dx'' = -B{r)df + C{r)dr^ + A{rydn\ (6.49) 

Plugging this ansatz into the equations of motion, we get the following from the tt equation, 
rr equation and 66 equation (which is the same as the 00 equation by spherical symmetry) 
respectively, 

4BC^m\^A^ + (2B{C - 3)C^m^r^ - AVA^BC (C - rC")) A^ 
+2^m3C {2C^ - 2r (3A' + rA!') C + r'^AC) A + C^A^BCr'^ {Af = 0, 

(6.50) 

4{B + rB') A^ + {2PA'B' - AB {C - rA')) A + Br^ {A')^ 2{2A + B- 3)m2 _ 

(6.51) 

-2B'^C^m^rA^ - 2B'^C\B + C- 3)m^rA^ 

-VA'^BC (2C'B^ + {rB'C - 2C {B' + rB")) B + Cr (5')^) A^ 

+BVA^BC {CrA'B' + B {AC A' - rC'A' + 2CrA")) A - B^CVA^BCr {A')^ = 0. 

(6.52) 

In the massless case, A{r) could be removed by a coordinate gauge transformation, and the 
last equation was redundant - it was a consequence of the first two. With non-zero m, there 
is no diffeomorphism invariance, so no such coordinate change can be made, and the last 
equation is independent. 

As we did in Section 6.1 for GR, we expand these equations around the flat space 
solution 

5o(r) = l, Co(r) = l, Ao{r) = l. (6.53) 
We introduce the expansion 

B{r) = Bo{r) + eBi{r) + e^B2{r) + --- , (6.54) 

C{r) = Co(r) + eCi(r) + e2C2(r) + -- - , 

A{r) = Ao{r) + eAi{r) + e'^A2{r) + --- . 

^^In general, when there are two metrics staticity and spherical symmetry are not enough to put both in 
diagonal form. An r dependent off diagonal drdt term can remain in one of them. We will not seek such 
off-diagonal metrics, and will limit ourselves to the diagonal ansatz. 
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Plugging into the equations of motion and collecting like powers of e, the O{0) part gives 
= because Bq,Co, Aq are solutions to the full non-linear equations. At each higher order 
in epsilon we will obtain a linear equation that lets us solve for the next term. At 0{e) we 
obtain 



2 (mV^ -l)Ai + (mV^ + 2) Ci + 2r {-3A\ + C[ - rA'l) = 0, (6.55) 
4B,,„^+(i-m^)A. + ^Mii^il^^0. (6.56) 
rAim'^ + rBim^ + rC^m^ - 2A[ -B[ + C[- rA[ - rB'[ = 0. (6.57) 



One way to solve these equations is as follows. Algebraically solve them simultaneously 
for Ai,A[,A'( in terms of -Bi's and Ci's and their derivatives. Then write the equations 
£Ai{B,C) = A[{B,C) and £A'{B,C) = A"{B,C). Solve these two equations for Ci and 
C[ in terms of i?i's its derivatives. Then write -^Ci{B) = C[{B). What is left is 

- SrBim'^ + 6B[ + SrB'l = 0. (6.58) 

There are two integration constants in the solution to (6.58), one is left arbitrary and the 
other must be sent to zero to prevent the solutions from blowing up at infinity. We then 
recursively determine Ci and Ai. Thus the whole solution is determined by two pieces of 
initial data^^. 



The solution is 



8GMe-'"'^ 

Bi{r) = , 6.59 

6 r 



8GMe— 1 + mr 

Ci r = —, 6.60 

6 r m'^r^ 

, , , 4GM e"'"'^ 1 + mr + 

Mr) = —- , (6.61) 



where we have chosen the integration constant so that we agree with the solution (3.15) 
obtained from the Green's function. 

We can now proceed to O(e^). Going through the same procedure, we find for the 



^^Naively, it is a second order equation in Ai and Bi, first order in Ci and we might think this requires 5 
initial conditions, but in fact it is a degenerate system, and there are second class constraints bringing the 
required boundary data to 2. 
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solution, when mr <C 1, 

8 GM f 1 GM \ 

CM-1 ^ -?^(l-14^....), ,6.63) 

4 GM f GM \ 
A{r)-l = ^ 1-4— ^ + --- . 6.64 

The dots represent higher powers in the non-hnearity parameter e. We see that the non- 
hnearity expansion is an expansion in the parameter ry/r, where 

gmVI^ 

is known as the Vainshtein radius. As the mass m approaches 0, ry grows, and hence the 
radius beyond which the solution can be trusted gets pushed out to infinity. As argued by 
Vainshtein [31], this perturbation expansion breaks down, and says nothing about the true 
non-linear behavior of massive gravity in the massless limit. Thus there is reason to hope 
that the vDVZ discontinuity is merely an artifact of linear perturbation theory, and that the 
true non-linear solutions show a smooth limit [31, 83, 84, 85]. 

One might hope that a smooth limit could be seen by setting up an alternative expan- 
sion in the mass m^. We take a solution to the massless equations, the ordinary Schwarzschild 
solution, with metric coefficients Bo,Go,Ao, and then plug an expansion 

B{r) = Bo{r) +m^Bi{r) +m^B2{r) + ■■■ , 
C{r) = Co{r) + m^Ciir) + m^C^ir) + ■ • • , 

A{r) = Ao{r) +m^Ai{r) +m^A2{r) + ■■■ , (6.66) 

into the equations of motion. Collecting powers of m yields a new perturbation equation at 
each order, but in this case the equations are generally non-linear. Even the equation we 
obtain at 0{m?) for the first correction to Schwarzschild is non-linear, so working with this 
expansion is much more difficult than working with the linearized expansion. 

The linearity expansion is valid is the region r ^ ry- If general relativity is restored at 
distances near the source, the mass expansion should be valid in the opposite regime r ^ ry, 
and the full solutions should interpolate between the two expansions. There have been 
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several extensive numerical studies of the full non-linear solutions. At first, troubles were 
encountered trying to find a complete and satisfactory solution that interpolates between the 
two regimes [86]. Later, the problem was revisited with more sophisticated methods, both 
in the decoupling limit [87], and more extensively in the full theory [88, 89, 90], with the 
final result being that the non-linearities can in fact work to restore continuity with GR. We 
will see later the mechanism by which this occurs. Other solutions, including some analytic 
solutions in various cases, are claimed in [91, 92, 93, 94, 95, 96]. 



6.4 Non-linear hamiltonian and the Boulware-Deser mode 

We now go on to study the hamiltonian of the non-linear massive gravity action (6.32) with 
fiat absolute metric r]^'^, 



2k^ 



(6.67) 



We saw in Section 2.1 that the free theory carries five degrees of freedom in D = 4, due 
to the fact that the time components /iqo appeared as a Lagrange multiplier in the action. 
We will see that this no longer remains true once the non-linearities of (6.67) are taken into 
account, so there is now an extra degree of freedom. 

A particularly nice way to study gravity hamiltonians is through the ADM formalism 
[97, 98]. A spacelike slicing of spacetime by hypersurfaces is chosen, and we change 
variables from components of the metric g^y to the spatial metric gij, the lapse iVj and the 
shift N , according to 

^700 = -N^ + g''N,N^, (6.68) 
9oi = Ni, (6.69) 
9ij = 9ij- (6.70) 

Here . . . are spatial indices, and g'^^ is the inverse of the spatial metric gij (not the ij 
components of inverse metric (7^^). 

The Einstein-Hilbert part of the action in these variables reads (see [99, 100] for detailed 
derivations and formulae) 



^ j d^'x^N [(")i? -K^ + K'^K,,] , (6.71) 
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where ^^''R is the curvature of the spatial metric gij. The quantity Kij is the extrinsic 
curvature of the spatial hypersurfaces, defined as 

Krj = ^ ig^J - V.iV, - V,iV,) , (6.72) 

where dot means time derivative, and the covariant derivatives are with respect to the spatial 
metric Qij. We then Legendre transform the spatial variables gij, defining the canonical 
momenta 

P'' = j^^=ir2V9{K''-Kg'^)^ (6.73) 

and writing the action in hamiltonian form 

2k'^L = (^j dH^x v"^m?j - H, (6.74) 
where the hamiltonian H is defined by 

H =(^j (fx p^^Qah^ ~L = j d'^x NC + NiC\ (6.75) 
and the quantities C and Q are 

C = ^ + K^- K'^K,,] , C = 2^Vj {K'^ - Kh'^) , (6.76) 

and here Kij should be thought of as a function of p^^ and g^j, obtained by inverting (6.73) 
for (jij and plugging into (6.72), 



V9\ 

All traces and index manipulations are performed with g^j and its inverse. 

For m = 0, the action is pure constraint, and the hamiltonian vanishes, a characteristic 
of diffeomorphism invariance. The shift and lapse Ni appear as Lagrange multipliers, 
enforcing the hamiltonian constraint C = and momentum constraints Ci = 0. It can be 
checked that these are first class constraints, generating the D diffeomorphism symmetries 
of the action. In D = 4, we have 12 phase space metric components, minus 4 constraints, 
minus 4 gauge symmetries, leaves 4 phase space degrees of freedom, the same counting as in 
the linear theory. The non-linear theory contains the same number of degrees of freedom as 
the linearized theory. 
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Now looking at the mass term, in ADM variables we have 



-2NH''hij + 2Ni [g'^ - 5'') Ni, (6.78) 



where hij = gij — 6ij. The action becomes 



S = ^ I d^'x p'^'gab -NC- NiC (6.79) 



2k' 



2 

m 



__ {h^^hi - K^h.i) + 25'^h,, - 2NH'^K, + 2N, [g^^ - 5'^) N,] . 

In the m 7^ case, the Fierz-Pauli term brings in contributions to the action that 
are quadratic in the lapse and shift (but still free of time derivatives). Thus the lapse and 
shift no longer serve as Lagrange multipliers, but rather as auxiliary fields, because their 
equations of motion can be algebraically solved to determine their values, 

C 



^.= ^{^"-S'rC'. (6.80) 



When these values are plugged back into (6.79), we have an action with no constraints 
or gauge symmetries at all, so all the phase space degrees of freedom are active. The resulting 
hamiltonian is 

(6.81) 

which is non-vanishing, unlike in GR. In 4 dimensions, we thus have 12 phase space degrees 
of freedom, or 6 real degrees of freedom. The linearized theory had only five degrees of 
freedom, and we have here a case where the non-linear theory contains more degrees of 
freedom than the linear theory. It should not necessarily be surprising that this can happen, 
because there is no reason non-linearities cannot change the constraint structure of a theory, 
or that kinetic terms cannot appear at higher order. 

As was argued in [32], the hamiltonian (6.81) is not bounded, and since the system is 
non-linear, it is not surprising that it has instabilities [81]. The nature of the instability, i.e. 
whether it is a ghost of a tachyon, what backgrounds it appears around, and its severity, is 
hard to see in the hamiltonian formalism. But in Section 8.2 we will see that this instability 
is a ghost, a scalar with a negative kinetic term, and that its mass around a given background 
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can be determined. It turns out that around flat space, the ghost degree of freedom is not 
excited because its mass is inflnite, but around non-trivial backgrounds its mass becomes 
finite. This ghostly extra degree of freedom is referred to as the Boulware-Deser ghost [32]. 

There is still the possibility that adding higher order interaction terms such as terms 
and higher, can remove the ghostly sixth degree of freedom. Boulware and Deser analyzed 
a large class of various mass terms, showing that the sixth degree of freedom remained [32] , 
but they did not consider the most general possible potential. This was addressed in [101], 
where the analysis was done perturbatively in powers of h. The lapse is expanded around 
its flat space values, = 1 + 6N. In this case, 6N plays the role of the Lagrange multiplier, 
and it is shown that at fourth order, interaction terms involving higher powers of 6N cannot 
be removed. It is concluded in [101] that the Boulware-Deser ghost is unavoidable, but this 
conclusion is too quick. It may be possible that there are fleld redeflnitions under which the 
lapse is made to appear linearly. Alternatively, it may be possible that after one solves for 
the shift using its equation of motion, then replaces into the action, the resulting action is 
linear in the lapse, even though it contained higher powers of the lapse before integrating out 
the shift. It is also possible that the lapse appears linearly in the full non-linear action, even 
though at any flnite order the action contains higher powers of the lapse. (For discussions 
and examples of these points, see [34, 102].) 

As it turns out, it is in fact possible to add appropriate interactions that eliminate the 
ghost [103]. In D dimensions, there is a — 2 parameter family of such interactions. We 
will study these in Section 9, where we will see that they also have the effect of raising the 
maximum energy cutoff at which massive gravity is valid as an effective fleld theory. ^'^ 

^■^Note that merely finding a ghost free interacting Lorentz invariant massive gravity theory is not hard - 
take for instance U{ri, h) — —2 [det ((5^^*^ + h^^ — /i] in (6.34), while letting the kinetic interactions be those 
of the linear graviton only. A hamiltonian analysis just like that of Section (2.1) shows that /loo and h^i 
both remain Lagrange multipliers. The problem is that this theory does not go to GR in the m — !■ limit, 
it goes to massless gravity. The real challenge is to construct a ghost free Lorentz invariant massive gravity 
that reduces to GR. 
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7 The non-linear Stiikelberg formalism 



In this section we will extend the Stiikelberg trick to full non-linear order. This will be 
a powerful tool with which to elucidate the non-linear dynamics of massive gravity. It will 
allow us to trace the breakdown in the linear expansion to strong coupling of the longitudinal 
mode. It will also tell us about quantum corrections, the scale of the effective field theory 
and where it breaks down, as well as the nature of the Boulware-Deser ghost and whether it 
lies within the effective theory or can be consistently ignored. 

7.1 Yang-Mills example 

We will first warm up with the spin 1 case, and we will set D = 4. The unique theory of 
interacting massless spin 1 particles is Yang-Mills theory [3]. Analogously to what we've 
done with gravity in Section 6.2, consider a non-abelian SU{N) gauge theory with gauge 
coupling g, and add a non-gauge invariant mass term with mass m for the gauge bosons, 
while leaving the kinetic structure unchanged from the massless case. 



The gauge fields are = —igAI^Ta, taking values in a Lie algebra with generators T^, with 
adjoint index a = 1, . . . , A^^ — 1. The generators satisfy the usual Lie algebra commutation 
and orthogonality relations [Ta,Tb] = i/^^'^Tc, Tr(TaTb) = ^6ab- The field strength is 
Ffj_iy = dfj,Ai^ — duAfj_ + [A^, Ay] = —igF'^^Ta. The theory (7.1) naively appears renormalizable, 
since there are no interaction terms with mass dimension greater than 4. But the propagators 
are those of a massive vector which do not go like ~ 1/p^, so naive power counting does not 
apply. 

In the absence of the mass term, the action is invariant under the gauge transformations 



where R = e"*""-^" G SU{N), and a°'{x) are gauge parameters. This reads infinitesimally 
(5y4^ = — ^dfj^a"' — fjj^"'A''^a'^. The field strength transforms covariantly — )■ RFfj^yR\ which 
reads infinitesimally SFj^^^ = fi,c"'(y^F^y 

The mass term breaks the gauge symmetry (7.2) (though it remains invariant under 
the global version), so we will restore it by introducing Stiikelberg fields. We pattern the 




(7.1) 



A^ -> RA^R^ + Rd^R\ 



(7.2) 
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introduction of the fields after the gauge symmetry we wish to restore, so we make the 
replacement 

^ UA^U^ + Ud^U\ (7.3) 

where 

U = e"^""^" G SU{N), (7.4) 

and the 7r"(x) are scalar Goldstone fields. The action now becomes gauge invariant under 
right gauge transformations^"', 

Af, RA^R^ + Rd^R\ U ^URl (7.5) 

The gauge kinetic term is invariant under this replacement, since it is gauge invariant, so 
the Goldstones appear only through the mass term 

2 2 
TD Ttl 

— Tr {A.A^ ^ -— Tr {D^U^D^U) , (7.6) 
9 9 

where D^U = d^U — U A^ is a covariant derivative, which transforms covariantly under right 
gauge transformations^'^, D^U {D^U)R^. We can go to the unitary gauge U = 1, and 
recover the massive vector action we started with, so the new action is equivalent. 

Expanding the terms in (7.6), we find kinetic terms for the vectors and scalars that 
require them to be canonically normalized as follows, 

A ~ gA, TT ~ — TT. (7.7) 

m 

Note that to lowest order in the fields, the non-linear Stiikelberg expansion (7.3) is the same 
as the linear one of Section 4.1. Thus the propagators all go like ~ ordinary power 

counting applies, and we can read off strong coupling scales from any non-renormalizable 
terms. 

For interactions, we have the usual normalizable Yang-Mills interaction terms with 
three and four fields, coming from the gauge kinetic term, 

~ g dA\ ~ g^A\ (7.8) 

^^Making the replacement — ?> WAfj,U — d^JJ would have led to left gauge transformations. 

^^The sigma model mass term (7.6) is invariant under SU{N)l x SU{N)ii global symmetry, U — ?► LUR\ 
of which the SU{N)ii part is gauged. The SU{N) subgroup L — R is realized linearly, and the rest is 
realized non-linearly. 
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From the mass term we find the non-renormahzable terms 

r^l^y'dH^, r^g(^y~'dAlT-, r^g^f^yAH'^. (7.9) 

Vm/ \m/ \m/ 

For g < 1, the lowest energy scale suppressing the non-renormalizable terms is ~ ^, 
which comes from the terms with only vr fields. The tree level amplitude for tttt — t- tttt 
scattering at energy E calculated from these terms goes like A ~ This amplitude 

becomes order one and unitarity is violated when E exceeds thus the Goldstones become 
strongly coupled at this energy, and this scale is the maximal cutoff for the theory, 

A=-. (7.10) 

9 

Note that when g is small this scale is parametrically larger than the vector masses m. 

We can take the decoupling limit which keeps this lowest scale fixed, while sending all 
the higher scales to infinity, 

g,m^O, A fixed. (7.11) 

The only terms that survives this limit are the scalar self-interactions (along with the free 
vector fields), 

5'decoupling = jd^x - A'Tt {d^U^d'^U) . (7.12) 

This is a limit which focuses in on the cutoff of the theory, ignoring all other scales^^. For 
this to be a valid limit, we should be looking at energies higher than the vector masses, 
and the coupling should be small. In this limit, the tt's becomes gauge invariant, but due 
to the way the Goldstones were introduced through traces of the combination UdfjU\ they 
retain a spontaneously broken SU{N)l x SU{N)ji global symmetry, U — )■ LUK^, of which 
the SU{N) subgroup L = R is realized linearly. 

Since we have an effective theory with cutoff A, there will be quantum corrections of 
all types compatible with the spontaneously broken SU{N)l x SU{N)^ global symmetry, 
suppressed by appropriate powers of the cutoff. For example we should find the operators 

~ Tr (^{d^U^d^uY^ , ~ Tr (d^U^d^U) , . . . (7.13) 

^^Note that the lowest scale is the only scale for which it is possible to take a decoupling limit. If we try 
to zoom in on a higher scale in a similar fashion, the terms with lower scales will diverge. 
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which in unitary gauge look hke the non-gauge invariant terms 

~Tr(A^), ~ Tr , . . . (7.14) 

Notice that the second operator in (7.14) modifies the gauge kinetic term in a non-gauge 
invariant way, and naively leads to ghosts. However, the mass of the ghost is rrig ~ m? / g"^ = 
A^, so it is safely at the cutoff. 

We might worry about the hierarchy between the small mass m and the high cutoff 
A ~ rn/g. If quantum corrections to the mass were to go like ~ A^, then the mass is 
pushed to the cutoff and there is a hierarchy problem which generally requires a solution 
in the form of fine-tuning of new physics at the cutoff. However, this does not happen 
here. There are only order one quantum corrections to mass ~ m^, coming from the 
generated operator — A^Tr {dfJJ^d^U^. Thus the small mass is technically natural, and can 
be consistently incorporated in the effective theory. 

This nice state of affairs is a consequence of the fact that gauge symmetry is restored 
in the limit as m — )■ 0, so that mass corrections must be proportional to the mass itself. For 
this to be true, it was important that there were no modifications to the kinetic structure 
of (7.1) not proportional to m, even though symmetry considerations would suggest that 
we are free to make such modifications. For example, suppose we try to calculate the mass 
correction to directly in unitary gauge by constructing Feynman diagrams with vertices 
read straight from (7.1). There are two interaction vertices ~ gdA^ and ~ g'^A'^ coming from 
the kinetic term. The mass term contributes no vertices but alters the propagator so that its 
high energy behavior is ~ At one loop, there are two IPl diagrams correcting the mass; 
one containing two cubic vertices and two propagators and one containing a single quartic 
vertex and a single propagator. Cutting off the loop at the momenta fcmax ~ A, the former 
diagram gives the largest naive correction 6m? ~ ^A^ ~ (The latter diagram gives the 
smaller correction 6m?' ~ ^A^ ~ A^.) This is above the cutoff, dangerously higher than the 
order one correction 6m'^ ~ m^ we found in the Goldstone formalism. 

What this means is that there must be a non-trivial cancellation of these leading 
divergences in unitary gauge, so that we recover the Goldstone result. This cancellation 
happens because the kinetic interactions of (7.1) are gauge invariant, implying that the 
dangerous k^k^ /m^ terms in the vector propagator do not contribute. Without these terms, 
the propagator goes like l//c^ and the estimate for the first diagram is 6m^ ~ S'^A^ ~ m^, in 
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agreement with the Goldstone prediction (the second diagram gives the smaller correction 
~ g^m^ log g). Non-parametrically altering the coefficients in the kinetic structure would 
spoil this cancellation and the resulting technical naturalness of the small mass (though such 
alterations could be done without spoiling technical naturalness if the alterations to the 
kinetic terms are suppressed by appropriate powers of m). For these reasons, it is desirable 
not to mess with the kinetic structure, and to introduce gauge symmetry breaking only 
through masses and potentials. The same will be true of massive gravity. 

This whole story, more than merely being a toy, can be thought of as a microcosm 
for the standard model. The fundamental particles seen so far in the electroweak sector are 
(the Higgs hasn't been seen as of this writing) spin 1/2 fermions and massive SU(2) spin 
1 gauge bosons (never mind the masslcss U(l) and complications of mixing). The gauge 
bosons masses are of order m ~ 10^ GeV, and the couplings g ~ 10^^. Their interactions 
at energies above m are well described by the above sigma model, up to an energy cutoff 
A ~ m/^f ~ 1 TeV. The reason for building the Large Hadron Collider is that something 
must happen at the scale A to UV complete the theory. 

If one demands that the UV completion be weakly coupled (as is suspected to be the 
case for the electroweak sector), one is led to introduce a new physical scalar, the Higgs, which 
unitarizes the amplitudes at energies above A. This UV completion is the standard model, 
a spontaneously broken gauge theory, where the Higgs has a mass and a perturbative 
quartic couphng A < 1, and gets a VEV v ~ ii/\f\ ~ A. The Higgs mass ~ \f\K sits 
somewhere between m ~ g'^v and the cutoff A. At the scale /i, the four point amplitude 
reaches the value ^ ~ A, the Higgs theory takes over, and the amplitudes cease growing with 
the energy, so that unitarity is not violated. Prom this perspective, studying the addition of 
a mass term to a gauge theory is not just an idle theoretical exercise. It leads one to uncover 
the Higgs mechanism and a new weakly coupled UV completion which is likely realized in 
nature. 

We will find an analogous story in the case of massive gravity. There is an effective 
field theory with a cutoff parametrically higher than the graviton mass, and the hierarchy is 
technically natural. The only missing part is the UV completion, which remains an unsolved 
problem. 



68 



7.2 Stiikelberg for gravity and the restoration of diffeomorphism 
invar iance 



We will now construct the gravitational analogue of the above. This method was brought to 
attention by [33, 104], but was in fact known previously from work in string theory [105, 106]. 

The full finite gauge transformation for gravity is (6.2), 

gt.u{x) ^ -Q^Q-;;9ap {f{x)) , (7.15) 

where f{x) is the arbitrary gauge function, which must be a diffeomorphism. In massive 
gravity this gauge invariance is broken only by the mass term. To restore it, we introduce a 
Stiickelberg field y(x), patterned after the gauge symmetry (7.15), and we apply it to the 
metric ^f^,,. 

The Einstein-Hilbert term yJ—gR will not change under this substitution, because it is 
gauge invariant, and the substitution looks like a gauge transformation with gauge parameter 
Y^{x), so no Y fields are introduced into the Einstein-Hilbert part of the action. 

The graviton mass term, however, will pick up dependence on s, in such a way that 
it will now be invariant under the following gauge transformation 

9A^) ^ -£il^9c., (fix)) , Y^\x) ^ {Y{x)r . (7.17) 

with /(x) the gauge function. This is because the combination Gfj_u is gauge mvariant (not 
covariant). To see this, first transform^^ g^i,, 

d.Y^d^Y^g^p {Y{x)) d^Y^d^Y^dJ^lYd^rirgx. UXyix))) , (7.21) 

^^The transformation of fields that depend on other fields is potentially tricky. To get it right, it is 
sometimes convenient to tease out the dependencies using delta functions. For example, suppose we have a 
scalar field (t>{x), which we know transforms according to (t>{x) — ^ 4'{f{x)). How should (j){Y{xy) transform? 
To make it clear, write 

c^[Y{x)) = j dy<P{v)6{y ~ Y{x)). (7.18) 
Now the field 4> appears with coordinate dependence, which we know how to deal with, 

dy(l>{f{ymy - Y{x)) = <j, {f{Y(x))) . (7.19) 
Going through an identical trick for the metric, which we know transforms as g^u{x) — >■ 'j^'j^gap (/(x)). 
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and then transform Y, 

^ d, [r\Y)Yd. [r\Y)Ydj'\f-^^y)dpr\f-^y)g,Ayi^)) 

= dp [rT\y9,y'dr [r'f \YduY^dj%-.(y)dpr\f-.^Y)gx. {y{x)) 

= (5^(5- d^Yfd,Y-gx, {Y{x)) = d^Y^d^Y'^gx, {Y{x)) . (7.22) 

We now expand Y about the identity, 

y"(x) + (7.23) 

The quantity G^uj, is expanded as 

= (5^ + 9^^") (5f + d,A^){g^^ + + \A^A^dxd,g^p + • • • ) 

+d^A''d,Af'g^p + d^A'^A^dpg^, + d^A'^A^df^g^^ + • • • (7.24) 

We now look at the infinitesimal transformation properties of ^f, Y , G, and Y, under 
infinitesimal general coordinate transformations generated hy f{x) — x + ^(x). The metric 
transforms in the usual way, 

= e^^A^^. + d^^gx^ + d^^^g^x. (7.25) 

The transformation law for the A's comes from the transformation of Y, 

Yf'ix) f-\Y{x))^' « F^(a;) - ^^{Y{x)), 
5Y^^-i^{Y), 

^ + A) = - A^'dae - ^A^A^dad^^" . (7.26) 

The A^^ are the Goldstone bosons that non-linearly carry the broken diffeomorphism invari- 
ance in massive gravity. The combination G^^, as we noted before, is gauge invariant 

5G^, = 0. (7.27) 

we find 

9a0 {Y{x)) ^ d^f\d0r\Ygxa {f{Y{x))) . (7.20) 
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We now have a recipe for Stiikelberg-ing the general massive gravity action of the form 
(6.34). We leave the Einstein-Hilbert term alone. In the mass term, we write all the /i^j^'s 
with lowered indices to get rid of the dependence on the absolute metric, and then we replace 
all occurrences of h^^^ with 

H,,{x) = G,,{x)-gl^X^)- (7.28) 

We then expand Gf^„ as in (7.24), and y as in (7.23). To linear order in hfj,i, = g^^ — gj^,) 
and Afj_, the expansion reads 

H^, = V + A, + V(°)/l^, (7.29) 

where indices on A are lowered with the background metric. This is exactly the Goldstone 
substitution we made in Section 5 in the linear case. 

In the case where the absolute metric is flat, gj^J = rj^y^ we have from (7.24), 

H^u = Ku + d^A, + d,A^ + d^A'^dyA^ + • • • (7.30) 

Here indices on A^ are lowered with rj^y and the ellipsis are terms quadratic and higher in 
the fields and containing at least one power of h. This takes into account the full non-linear 
gauge transformation. 

As in the linear case, we will usually want to do another scalar Stiikelberg replacement 
to introduce a f/(l) gauge symmetry, 

A^^A^ + d^<p. (7.31) 
Then the expansion for flat absolute metric takes the form 

= /i^,+a^A,+a,A^+2a^a,</.+a^A"a,A,+a^A"a,a„</.+a^a°</.a,A,+a^a"</.a,a,</.+- ■ ■ , 

(7.32) 

where again the ellipsis are terms quadratic and higher in the fields and containing at least 
one power of h. The gauge transformation laws are (7.26), (6.8), 

5h^y = d^^y + dy^^ + C^h^y, 
6A, = d^A-^^-A^^d^^.-^A'^A^do^d^C,---- , 
6(j) = -A. (7.33) 
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This method of Stiikelberg-ing can be extended to any number of gravitons and general 
coordinate invariances, as was done in [33, 107], in analogy with the gauge theory little Higgs 
models and dimensional de-construction [108, 109]. When multiple gravitons are present, all 
but one must become massive, since there are no non-trivial interactions between multiple 
massless gravitons [110], and these gravitons mimic the Kaluza-Klein spectrum of a discrete 
extra dimension. Other work in this area, including applications to bi-gravity and multi- 
gravity models, can be found in [111, 112, 113, 114, 50, 115]. 

7.3 Another way to Stiikelberg 

In the last section, we introduced gauge invariance and the Stiikelberg fields by replacing 
the metric g^^, with the gauge invariant object G^,^. This is well suited to the case where we 
have a potential arranged in the form (6.34), because all the background g^^^f^'^^s appearing 
in the contractions and determinant of the mass term do not need replacing. The drawback 
is that the Stukelberg expansion involves an infinite number of terms higher order in h^^. If 
we wish to keep track of the h^i,^s, this is not very convenient. 

Instead, we develop another method, which is to introduce the Stiikelberg fields through 
the background metric gj^J , and then allow g^^^, to transform covariantly. This method will 
be better suited to a potential arranged in the form (6.40), and will have the advantage that 
the Stiikelberg expansion contains no higher powers of /i^;^. 

We make the replacement 

^ {Y{x)) d,Y-d.Y^. (7.34) 

The y"(x) that are introduced are four fields, which despite the index a, are to transform 
as scalars under diffeomorphisms 

F"(x) ^ r"(/(x)), (7.35) 

or infinitesimally, 

= ^'^S^F". (7.36) 

This is to be contrasted with the transformation rule SY"" = S,'^ dvY'^ ~ {8^^°") Y^ which would 
hold if were a vector. Given this scalar transformation rule for y", the replaced g'^^} now 
transforms like a metric tensor. If we now assign the usual diffeomorpshim transformation law 
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to the metric g^u (so that it is now covariant), quantities hke g^^}g^^ and other contractions 
will transform as diffeomorphism scalars. We can take any action which is a scalar function 
of g'f^} and g^^, and introduce gauge invariance in this way^*^. 

This is convenient when we have a potential of the form (6.40). First we lower all 
indices on the /i^i/'s in the potential. Now the background metric g'^l^y appears only through 
h^v = g^v — gfvi so we replace all occurrences of h^y with 

H,.{x) = g,,{x) - g'^l (F(x)) d.Y^d.Y''. (7.37) 

Note that we need make no replacement on the g^^^^s used to contract the indices, or on the 
y/—g out front of the potential in (6.40). 

Expanding, 

= - A", (7.38) 

and using g^^ = gf} + h^^y, we have 

= V + 9?ld,A- + gfldyA" - g'^^jd.A^dyA^ + ■ ■ ■ , (7.39) 

where the ellipses denote terms that contain derivatives of g^l^} (and so vanish in the usual 
case of interest where g'f^} = rj^y). Note the difference in sign for the term quadratic in 
compared with (7.30). 

Under infinitesimal gauge transformations we have 

= -C + CduA"", (7.40) 

= V(,°)e. + Vl^^^M + A V' (7-41) 

where the covariant derivatives are with respect to g'^^} and the indices on are lowered 
using g^l^} . To linear order, the transformations are 

= -r, (7.42) 
= Vfe. + VW^M, (7.43) 



^^This is essentially the technique of spurion analysis, where a coupling constant is made to transform as 
a field. A quantity which is normally a background quantity, a coupling constant in the case of spurious, or 
the background gj^J in this case, is made to transform in some way that gives the action more symmetries. 
Note that this method of introducing gauge invariance can be carried out on any Lorentz invariant action, 
even one that does not contain a dynamical metric gf^^. For example, a plain old scalar field in flat space can 
be made diffeomorphism invariant in this way. This highlights the fact that general coordinate invariance is 
not the critical ingredient that leads one to a theory of gravity, since it can be made to hold in any theory. 
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which reproduces the hnear Stiikelberg expansion used in Section 5. 
In the case of a flat background, gjj^} = rj^^, the replacement is 

H^. = h^, + d^A, + d,A^ - d^A^d.A^, (7.44) 

with indices on A" lowered by ri^^y. Notice that this is the complete expression, there are no 
higher powers of h^y^ unlike (7.30). 

We will often follow this with the replacement — )■ + to extract the helicity 
mode. The full expansion thus reads 

H^^u = Ku + d^A, + d,A^ + 2d^d,<P - d^A'^d.Aa - d^A''d,d^<P - d^d''(j)d,A^ - d^d''cj,d,da(j>. 

(7.45) 

Under infinitesimal gauge transformations, 

5h^y = df,^y + dy^^ + C^hf,^, (7.46) 
6A^ = d^A-^^ + Cd,A^, (7.47) 
6(j) = -A. (7.48) 

Yet another way to introduce Stiikelberg fields is advocated in [116, 117, 118], in which 
they make the inverse metric g'^'^ covariant through the introduction of scalars g^'^{x) — ?► 
gOLji {Y~^{x)) daY^dfiY'^ . There have also been many studies, initiated by 't Hooft, of the so- 
called gravitational Higgs mechanism, which is also essentially a Stiikelberg-ing of different 
forms of massive gravity [119, 120, 121, 122, 123, 124, 125, 126, 127]. All of these are 
equivalent to the theories we study, as can be seen simply by going to unitary gauge [128]. 
At the end of the day, (6.34) is the most general Lorentz invariant graviton potential, and any 
Lorentz invariant massive gravity theory will have a unitary gauge with a potential which is 
equivalent to it for some choice of the coefficients Ci, C2, etc. 



8 Stiikelberg analysis of interacting massive gravity 

In this section, we will set D = 4 and apply the Stiikelberg analysis to the massive GR 
action (6.32) in the case of a fiat absolute metric. The mass term reads 

Smass = I dWt^ {Kc^Kp - K-h^p) ■ (8-1) 
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The Stiikelberg analysis instructs us to make the replacement (7.32), 

h^^u -> H^, = V + -9^^. + d,A^ + d^A'^d.A^ + 29^9,0 + d^d'^cpOM ■■■ . (8.2) 

The extra terms with h in the ellipsis will not be important for this theory, as we will see. 

At the linear level, this replacement is exactly the linear Stiikelberg expansion of Section 
4. We will have to canonically normalize the fields here to match the fields of the linear 
analysis. Using a hat to signify the canonically normalized fields with the same coefficients 
as used in Section 4 (although there we omitted the hats), we have 

1 - 1 - 1 

h = -Mph, A = -mMpA, = -m^Mp4>. (8.3) 

2 2 2 

We also get a whole slew of interaction terms, third order and higher in the fields, 
suppressed by various scales. We always assume m < Mp. cj) always appears with two 
derivatives, A always appears with one derivative, and h always appears with none, so a 
generic term, with rih powers of h^^, ua powers of A^ and n<^ powers of 0, reads 

~ m2M|/i"'^(M)"^(920)"^ ~ AJ"""~^"^"^"^/^"''(9i)"^(920)"^ (8.4) 

where the scale suppressing the term is 

a, = (aW-')"\ , = 3!!i±^!M±i!^. (8,5) 

n^ + nA + nh-2 

The larger A, the smaller the scale, since m < Mp. We have + + > 3, since we are 
only considering interaction terms. The term suppressed by the smallest scale is the cubic 
scalar term, n,^ = 3, ua = Uh = 0, which is suppressed by the scale A5 = (Mpm'^y^^, 



(av)' 



As = {Mpmy/\ (8.6) 



In terms of the canonically normalized fields (8.3), the gauge symmetries (7.33) read 

SA, = d,A-mi, + ^td^A,--^^A-PdM,---- , 
S(j) = -mk, (8.7) 
where we have rescaled A = ^A and C'' = ^C''- 
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Finally, note that since the scalar field always appears with at least two derivatives 
in the Stiikelberg replacement (8.2), the resulting action is automatically invariant under the 
global galilean symmetry 

(f){x) ^ c + b^x^, (8.8) 

where c and 6^ are constants. In addition, the action is automatically invariant under global 
shifts in — )■ + for constant c^. These symmetries are the gravitational analogs of 
the SU{N)l X SU{N)ii global symmetry of the Yang-Mills model in Section 7.1. It will 
persist even in limits where the gauge symmetries on and (p no longer act. 



8.1 Decoupling limit and breakdown of linearity 

As seen in Section 4.2, the propagators have all been made to go like ~ so normal 

power counting applies, and the lowest scale, A5, is the cutoff of the effective field theory. 
To focus in on the cutoff scale, we take the decoupling limit 

T 

m ^ 0, Mp ^ 00, T ^ 00, A5, — fixed. (8.9) 

Adp 

All interaction terms go to zero, except for the scalar cubic term (8.6) responsible for the 
strong coupling, which we may calculate using the replacement H^^, = 2d^d^(f) + d^d'^cpd^dafj) 
since we do not need the vector and tensor terms. As discussed in Section 4.2, we must 
also do the conformal transformation /i^^ = /i'^^ + m?(j)ri^u. This will diagonalize all the 
kinetic terms (except for various cross terms proportional to m which are eliminated with 
appropriate gauge fixing terms, as discussed in Section 4.2, and which go to zero anyway in 
the decoupling limit). 

After all this, the lagrangian for the scalar reads, up to a total derivative, 

S^ = j d'x - 3(90)2 ^ _! Y^^Y _ (□0)(5^5^^)2j + (8.10) 

The free graviton coupled to the source via j^h'^^T'^^ also survives the limit, as does the 
free decoupled vector. 

We can now understand the origin of the Vainshtein radius at which the linear expansion 

breaks down around heavy point sources. The scalar couples to the source through the trace, 

■j^(j)T. To linear order around a central source of mass M, we have 

Ml , , 

8.11 
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The non-linear term is suppressed relative to the linear term by the factor 

d^6 M 1 



Mp Af r5 ■ 



(8.12) 



rv-{i^] --fe) • (8.13) 



Non-linearities become important when this factor becomes of order one, which happens at 
the radius 

M 1 /GM^ 

When r < ry, linear perturbation theory breaks down and non-linear effects become impor- 
tant. This is exactly the Vainshtein radius found in Section 6.3 by directly calculating the 
second order correction to spherical solutions. 

In the decoupling limit, the gauge symmetries (8.7) reduce to their linear forms, 

5(t) = 0. (8.14) 

Even though (p is gauge invariant in the decoupling limit, the fact that it always comes with 
two derivatives means that the global galileon symmetry (8.8) is still present, as is the shift 
symmetry on A^. 



8.2 Ghosts 

Note that the lagrangian (8.10) is a higher derivative action, and its equations of motion are 
fourth order. This means that this lagrangian actually propagates two lagrangian degrees of 
freedom rather than one, since we need to specify twice as many initial conditions to uniquely 
solve the fourth order equations of motion [49], and by Ostrogradski's theorem [129, 130], 
one of these degrees of freedom is a ghost. The decoupling limit contains six degrees of 
freedom - two in the massless tensor, two in the free vector, and two in the scalar. This 
matches the number of degrees of freedom in the full theory as determined in Section 6.4, so 
the decoupling limit we have taken is smooth. The extra ghostly scalar degree of freedom is 
the Boulware-Deser ghost. Note that at linear order, the higher derivative scalar terms for 
the scalar are not visible, so the linear theory has only 5 degrees of freedom. 
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Following [101], let's consider the stability of the classical solutions to (8.10) around 
a massive point source. We have a classical background $(r), which is a solution of the 
(f) equation of motion, and we expand the lagrangian of (8.10) to quadratic order in the 
fluctuation ^9 = — $. The result is schematically 

£,~-(9^)^ + ^(9V)^. (8.15) 

There is a four-derivative contribution to the ip kinetic term, signaling that this theory 
propagates two linear degrees of freedom. As shown in Section 2 of [101], one is stable and 
massless, and the other is a ghost with a mass of order the scale appearing in front of the 
higher derivative terms. So in this case the ghost has an r-dependent mass 

'^ghostlO ~ • (8.16) 

This shows that around a flat background, or far from the source, the ghost mass goes 
to infinity and the ghost freezes, explaining why it was not seen in the linear theory. It 
is only around non-trivial backgrounds that it becomes active. Notice, however, that the 
backgrounds around which the ghost becomes active are perfectly nice, asymptotically fiat 
configurations sourced by compact objects like the Sun, and not disconnected in any way in 
field space (this is in contrast to the ghost in DGP, which occurs only around asymptotically 
de Sitter solutions). 

We are working in an effective field theory with a UV cutoff A5, therefore we should 
not worry about instabilities until the mass of the ghost drops below A5. This happens at 
the distance rghost where d^(^'' ~ A5. For a source of mass M, at distances r ^ ry the 
background field goes like ^(r) ~ ]^^5 so 



'"ghost is parametrically larger than the Vainshtein radius ry. 

As we will see in Section 8.4, the distance rghost is the same distance at which quantum 
effects become important. Whatever UV completion takes over should cure the ghost insta- 
bilities that become present at this scale, so we will be able to consistently ignore the ghost. 
We see already that we cannot trust the classical solution even in regions parametrically 
farther than the Vainshtein radius. The best we can do is make predictions outside rghost, 
and we will have more to say later about this. 
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8.3 Resolution of the vDVZ discontinuity and the Vainshtein mech- 
anism 



We are now in a position to see the mechanism by which non-hnearities can resolve the vDVZ 
discontinuity. This is known as the Vainshtein mechanism. It turns out to involve the ghost 
in a critical role. 

Far outside the Vainshtein radius, where the linear term of (8.10) dominates, the field 
has the usual Coulombic 1/r form. But inside the Vainshtein radius, where the cubic term 
dominates, it is easy to see by power counting that the field gets an r^/^ profile. 



At distances much below the Vainshtein radius, the ghost mass (8.16) becomes very 
small, and the ghost starts to mediate a long range force. Usually a scalar field mediates 
an attractive force, but due to the ghost's wrong sign kinetic term, the force mediated by 
it is repulsive. In fact, it cancels the attractive force due to the longitudinal mode, the 
force responsible for the vDVZ discontinuity, and so general relativity is restored inside the 
Vainshtein radius. 

We will now see this more explicitly. Following [131], some field re-definitions can 
be done on the scalar action (8.10), and the result is an action schematically of the form 
C = -{d^f + {d^ljf + ^ ^ A^^^- Here is the heahhy longitudinal mode, 

ip is the ghost mode, and the original scalar can be found from cj) = cj) — ip. Both are 
coupled gravitationally to the stress tensor. Note that the self-interactions appear in these 
variables as a peculiar non-analytic term (we can also see that the ghost mass around 
a background {ip) will be Ag''^/ (■?/') ^/^). The field is free and has the profile ~ 
everywhere, mediating an attractive force. 

The %l) field however has two competing terms, which becomes comparable at the Vain- 
shtein radius. The linear term dominates at radii smaller than the Vainshtein radius, so 
■0 ~ for r ^ ry. This profile generates a repulsive Coulomb force that exactly 

cancels the attractive force mediated by 0, so in sum there are no extra forces beyond 
gravity in this region. (The leading correction to the profile is found by treating the -0^/^ 
term as a perturbation, -0 ~ -00 + ip{i) + ■ ■ ■ , with ipQ ~ ^ , plugging in the equation 




(8.18) 
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of motion d'^ijj^i) + ip^^^^ = obtaining ^/^(i) ~ y-^j A5 r^/^, in agreement with 
(8.18).) The funny non-hnear term dominates at radii larger than the Vainshtein radius, so 
ip ~ AfrS r ry, and so the ghost profile is negligible in this region compared 

to the (p profile. Thus the ghost ceases to be active beyond the Vainshtein radius, and the 
longitudinal mode generates a fifth force. 

This is known as a screening mechanism, a mechanism for rendering a light scalar 
inactive at short distances through non-linearities (see the introduction and references in 
[132, 133], and in a different context [134]). 

One can think of this as a kind of classical version of a weakly coupled UV completion 
via a Higgs. Above the Vainshtein radius (low energies), there is only the long distance scalar, 
which starts to become non-linear (strongly coupled) around the Vainshtein radius, so one 
can think of this regime in terms of an effective field theory with cutoff the Vainshtein radius. 
Below the Vainshtein radius (high energies), a new degree of freedom, the ghost (analogous 
to the physical Higgs in the standard model), kicks in. Much below the Vainshtein radius, 
everything is again linear and weakly coupled, with the difference that there are now two 
active degrees of freedom, so one can think of this as a classical UV completion of the effective 
theory. 

Of course, this ghostly mechanism for restoring continuity with GR relies on an in- 
stability, which would become apparent were we to investigate small fluctuations beyond 
the gross-scale features described here. Furthermore, as we will see in the next section, the 
ghost issue is moot, since the classical mechanism described in this section occurs outside 
the regime of validity of the quantum effective theory and is swamped by unknown quantum 
corrections. 



8.4 Quantum corrections and the effective theory 

Quantum mechanically, massive gravity is an effective field theory, since there are non- 

renormalizable operators suppressed by the mass scale A5. The amplitude for tttt — t- tttt 

/ N 10 

scattering at energy E, coming from the cubic coupling in (8.10), goes like ^ ~ f . 
This amplitude should correspond to the scattering of longitudinal gravitons. The wave 
function of the longitudinal graviton (2.20) for a large boost is proportional to m~^, while 



80 



the largest term at high momentum in the graviton propagator (2.44) is proportional to m~^, 
so naive power counting would suggest that the amplitude at energies much larger than m 
goes like A 

yw-2 — However, as recognized in [^33], and calculated explicitly in [^135], there 
is a cancellation in the diagrams so that the result agrees with the result of the Goldstone 
description. We will encounter these kinds of cancellations again in loops, and part of the 
usefulness of the Goldstone description is that they are made manifest. 

The amplitude becomes order one and hence strongly coupled when ~ A5. Thus A5 
is the maximal cutoff of the theory. We expect to generate all operators compatible with the 
symmetries, suppressed by appropriate powers of the cutoff. In the unitary gauge, there are 
no symmetries, so we will generate all operators of the form 

Cp^.d'^hP. (8.19) 

We wish to determine the scales in the coefficient Cp g. 

After Stiikelberg-ing, the decoupling limit theory contains only the scalar 0, and the 
single coupling scale A5. In addition, there is the galileon symmetry — t- + c + c^x^. 
Quantum mechanically, we expect to generate in the quantum effective action all possible 
operators with this symmetry, suppressed by the appropriate power of the cutoff A5. The 
galileon symmetry forces each to carry at least two derivatives"*^^, so the general term we 
can have is 



(8.20) 



To compare to (8.19), we go back to the original normalization for the fields by replacing 
~ m?Mp(j) and recall that d^j,dv(j) comes from an h^y to find that in unitary gauge the 
coefficients Cp^q go like 

Cp,, ~ A-^P-^+^M^m^P = (mi6-^''-2f M^^-'^+^) . (8.21) 

This comparison is possible because the operations of taking the decoupling limit and com- 
puting quantum corrections should commute. 



^^Actually, there are a finite number of terms whiclr liave fewer than two derivatives per field, the so-called 
galileon terms [136] which change by a total derivative under the galileon symmetry [136]. However, there is 
a non-renormalization theorem that says these are not generated at any loop by quantum corrections [137], 
so we need not include them. We will encounter them later when we raise the cutoff to A3. 
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Notice that the term with p = 2, g = Oisa mass term, ~ h'^, corresponding to 

a mass correction = (^^)- This is down by a factor of m^/Ag from the tree level 
mass term. Thus a small mass graviton mass m ^ A5 is technically natural, and there is no 
quantum hierarchy problem associated with a small mass. This is in line with the general 
rule of thumb that a small term is technically natural if a symmetry emerges as the term 
is dialed to zero. In this case, it is the diffeomorphism symmetry of GR which is restored 
as the mass term goes to zero. The quantum mass correction will also generically ruin the 
Fierz-Pauli tuning, but its coefficient is small enough that ghost /tachyons associated to the 
tuning violation are postponed to the cutoff - indeed the resulting ghost mass, using the 
relations in the second paragraph of Section 2, is ~ A5. 

Similar to the spin 1 case in Section 7.1, it is important that there were no non- 
parametric modifications to the kinetic structure of the Einstein-Hilbert term, even though 
the lack of gauge symmetry would suggest that we are free to make such modifications. 
Suppose we try to calculate the mass correction directly in unitary gauge. The graviton mass 
term contributes no vertices but alters the propagator so that its high energy behavior is ~ ^ 
(the next leading terms go like ^ and then p-). At one loop, there are two IPI diagrams 
correcting the mass; one containing two cubic vertices -^d'^h? from the Einstein-Hilbert 
action and two propagators, and another containing a single quartic vertex -^d'^h* from 
the Einstein-Hilbert action and a single propagator. Cutting off the loop at the momenta 
^max ~ A5, the first diagram gives the largest naive correction ~ j^^^s A^^ ~ A5. (The 
second diagram gives a smaller correction.) This is at the cutoff, dangerously higher than 
the small correction ~ Tn?^ we found in the Goldstone formalism. 

This means that there must be a non-trivial cancellation of this leading divergence in 
unitary gauge, so that we recover the Goldstone result. This cancellation happens because 
the kinetic interactions of Einstein-Hilbert are gauge invariant, implying that the dangerous 
k^k^k'^kl^ /m^ terms in the graviton propagator do not contribute. Without these terms, the 
propagator goes like l/rn? and the estimate for the first diagram is 5m? ~ 4 Ac ~ rn?^, 
in agreement with the Goldstone prediction (again the second diagram again gives a smaller 
correction). Non-parametrically altering the coefficients in the kinetic structure would spoil 
this cancellation and the resulting technical naturalness of the small mass (though such 
alterations could be done without spoiling technical naturalness if the alterations to the 
kinetic terms are parametrically suppressed by appropriate powers of m). These kinds of 
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cancellations can be seen explicitly in the calculations of [135]. Some loop calculations for 
massive gravity have been done in [138, 139]. 

In summary, in unitary gauge the theory (6.32) in D = 4 is a natural effective field 
theory with a cutoff parametrically larger than the graviton mass, with the effective action 



S 



I d'x ^ \^gR - "^{hl - h')] + J2 cv,,d'h^. (8.22 



and a cutoff A5 = {m^Mpf^. 

We should take into account the effect that the unknown quantum operators have on 
the solution around a heavy source. Given that the linear field goes like ~ ]^^, the radius 
g at which the term (8.20) becomes comparable to the kinetic term (90)^ is 

This distance increases with p, and asymptotes to its highest value 

M V^^ 1 



Thus we cannot trust the classical solution at distances below tq, since quantum op- 
erators become important there. This distance is parametrically larger than the Vainshtein 
radius, where classical non-linearities become important. Unlike the case in GR, there is 
no intermediate regime where the linear approximation breaks down but quantum effects 
are still small, so there is no sense in which a non-linear solution to massive gravity can be 
trusted for making real predictions in light of quantum mechanics. 

In particular, the entire ghost screening mechanism of Section 8.3 is in the non-linear 
regime, and so it becomes swamped in quantum corrections. Thus there is no regime for 
which GR is a good approximation - the theory transitions directly from the linear classical 
regime with a long range fifth force scalar, to the full quantum regime. Note that it is the 
higher dimension operators that become important first, so there is no hope of finding the 
leading quantum corrections. Finally, the radius rg is the same as the radius rghost where 
the ghost mass drops below the cutoff, so it is consistent to ignore the ghost since it lies 
beyond the reach of the quantum effective theory. The various regions are shown in Figure 
3. Note that in the decoupling limit we are working in, the Schwarzschild radius (and the 
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radii associated to all scales larger than A5) are sent to zero, while the scale ~ :^ where 
Yukawa suppression takes hold is sent to infinity. 



Quantum 


Classical 


Non-linear (classically) 


Linear (classically) 


Ghost 






r 


/ M Y'^ 1 


/M\'/^ 1 






ry ~ 10^^ km 


rg ~ 10^^ km 



Figure 3: Regimes for massive gravity with cutoff A5 = (Mpm^)^/^, and some values witliin the 
solar system, for which A^^ ~ 10"*^^ km. Note that rq is a bit larger than the observable universe, 
i.e. this theory makes no observable predictions within its range of validity. 



9 The A3 theory 

We have seen that the theory (6.32) containing only the linear graviton mass term has some 
undesirable features, including a ghost instability and quantum corrections that become 
important before classical non-linearities can restore continuity with GR. In this section, 
we consider the higher order potential terms in (6.34) and ask whether they can alleviate 
these problems. It turns out that there is a special choice of potential that cures all these 
problems, at least in the decoupling limit. 

This choice also has the advantage of raising the cutoff. With only the Fierz-Pauli 
mass term, the strong coupling cutoff was set by the cubic scalar self coupling ~ a|^- "^^^ 
cutoff A5 = (Mpm^Y^^ is very low, and as we will see, generically any interaction term will 
have this cutoff. But by choosing this special tuning of the higher order interactions, we end 
up raising the cutoff to the higher scale A3 = {Mpui'^Y^^. 

It was recognized already in [33] , that if the scalar self-interactions could be eliminated, 
the cutoff would be raised to A3. This was studied more fully in [101], where the cancelation 
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was worked through and it was (mistakenly) concluded that ghosts would be unavoidable 
once the cutoff was raised. Motivated by constructions of massive gravity with auxiliary 
extra dimensions (see Section 10.3), this was revisited in [34, 102], where the decoupling 
limit lagrangian was calculated explicitly and was seen to be ghost free. In [103], the full 
theory was shown to be ghost free^°. 



9.1 Tuning interactions to raise the cutoff 

Looking back at the scales (8.5), the term suppressed by the smallest scale is the cubic scalar 
term, which is suppressed by the scale A5 = (Mpm"^)^/^, 

(91) 

The next highest scale is A4 = (Mpin^Y^'^, carried by a quartic scalar interaction, and a 
cubic term with a single vector and two scalars. 

The next highest is a quintic scalar, and so on. The only terms which carry a scale less than 
A3 = {Mpm'^Y^^ are terms with only scalars (9^0)", and terms with one vector and the rest 
scalars dA(d'^(j))'^. 

The scale A3 is carried by only the following terms 



(9.3) 



All other terms carry scales higher than A3. 

It turns out that we can arrange to cancel all of the scalar self couplings by appropriately 
choosing the coefficients of the higher order terms. We will work with the form of the 
potential in (6.40) where indices are raised with the full metric, and the Stiikelberg formalism 
of Section 7.3. We do because we eventually want to keep track of powers of h, so the form 
of the Stiikelberg replacement in Section 7.3 is simpler. We are interested only in scalar self 
interactions, so we may make the replacement (7.45) with the vector field set to zero, 

H^, ^ 2 d^d,(l) - d^da(l) d.d'^cj). (9.4) 
^'^The objections of [118] and [140] are addressed in [141] and [142] respectively. 
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The interaction terms are a function of the matrix of second derivatives n^j, = d^d^jcf). 
As reviewed in Appendix A, there is at each order in a single polynomial in II^j, which 
is a total derivative. By choosing the coefficients (6.40) correctly, we can arrange for the cf) 
terms to appear in these total derivative combinations. The total derivative combinations 
have at each order in cf) as many terms as there are terms in the potential of (6.40), so all the 
coefficients must be fixed, except for one at each order which becomes the overall coefficient 
of the total derivative combination. 

The choice of coefficients in the potential (6.40) which removes the scalar self interac- 
tions is, to fifth order [34], 

1 1 , , 

C1 = 2C3+-, C2 = -3C3--, (9.5) 

di = -6rf5 + 4(24c3 + 5), ^2 = 8^5-7(603 + 1), 
16 4 

d3 = 34-^(1203 + 1), d, = -6d, + ^cs, (9.6) 
16 4 



/l = M + 1^3 - + 24/7 , /2 = -1 - iic3 + 64 - 30/7 , 

/s = ic3 - 34 + 20/7 , U = -h- + 54 - 20/7 , (9.7) 

h = I.C3 - 3^5 + 15/7 , h = d5- 10/7 . 

At each order, there is a one-parameter family of choices that works to create a total deriva- 
tive. Here C3, d^ and fj are chosen to carry that parameter at order 3, 4 and 5 respectively. 
Note however that at order 5 and above (or D + 1 and above if we were doing this in D 
dimensions), there is one linear combination of all the terms, the characteristic polynomial of 
h mentioned below (6.39), that vanishes identically. This means that one of the coefficients 
is redundant, and we can in fact set d^ and its higher counterparts to any value we like 
without changing the theory. Thus there is only a two parameter family {D — 2 parameter 
in dimension D) of theories with no scalar self-interactions. This can be carried through at 
all orders, and at the end there will be no terms ~ (9^(/))". 

The only terms with interaction scales lower than A3 were the scalar self-interactions 
(9^0)"", and terms with one vector and the rest scalars dA{d'^4^)^. We have succeeded in 
eliminating the scalar self- interactions, but since these always came from combinations ( + 
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dfj_) the terms dA{d'^(f))'^ are automatically of the form d^A'^xju) , where the xjn) are the 
functions of d^^d^cj) described in Appendix A, which are identically conserved d^X^^l) = 
0. Thus, once the scalar self-interactions are eliminated, the dA^d'^cj))"' terms are all total 
derivatives and are also eliminated. 

Now the lowest interaction scale will be due to the terms in (9.3), 



(9.8) 



which are suppressed by the scale A3 = {Mpm'^Y^^, so the cutoff has been raised to A3, 
carried by the terms (9.8). 

The decoupling limit is now 

m 0, Mp 00, A3 fixed, (9.9) 

and the only terms which survive are those in (9.3). To find these terms we must now go 
back to the full Stiikelberg replacement (7.45), and we must also expand the inverse metric 
and determinant in the potential of (6.40) in powers of h. The h{d'^(j))"' terms, up to quintic 
order in the decoupling limit, and up to total derivatives are [34], 



-4X(L)(0) + ^^^X(^)(0) + ^^^^f^X(^(0)" 



^ K.T^^r (9.10) 



Here the xj^^ are the identically conserved combinations of d^dvip described in Appendix 
A. The ((9/1)2(920) terms can be found to cubic order in [143]. The terms with A's can in 
any case be consistently set to zero at the classical level, since they never appear linearly 
in the lagrangian, so we will focus only on the terms involving h and (f). Properties of this 
lagrangian, including its cosmological solutions, degravitation effects and phenomenology are 
studied in [144]. Spherical solutions are studied in [145]. The cosmology of a covariantized 
version is studied in [146]. 

In terms of the canonically normalized fields (8.3), the gauge symmetries (7.48) of the 
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full theory are 



6A^ = d^k-mi^ + ^td,A^, (9.11) 

IVl p 

2 

= d^i^ + d^i^ + —Cfhf,y, (9.12) 

IVlp 

5<p = -mk, (9.13) 

where we have rescaled A = '"'^^^ A and = In the decoupling limit (9.9), this gauge 

symmetry reduces to its linear form, 

SA^ = d^A, (9.14) 
5h^u = d^t + d,i^, (9.15) 
Scj) = 0. (9.16) 

The lagrangian (9.10) should be invariant under the decoupling limit gauge symmetries 
(9.16). Indeed, the identity d^'X^'u' = ensures that it is. The scalar is gauge invariant in 
the decoupling limit, but the fact that it always comes with two derivatives means that the 
global galileon symmetry (8.8) is still present, as is the shift symmetry on A^. 

Note that for the specific choice C3 = 1/6 and = —1/48, all the interaction terms 
disappear. This could mean that the theory becomes strongly coupled at some scale larger 
than A3, or there could be no lowest scale, since there are scales arbitrarily close to but above 
A3. In the later case, the theory would have no non-linear behavior, and so no mechanism 
to recover continuity with GR, and it would therefore be ruled out observationally. 



9.2 Exactness of the decoupling limit 

In [102], a nice trick was used to show that the decoupling limit lagrangian (9.10) is exact 
to all orders in the fields, that is, there are no further terms h{d'^(f))"' for n > 4. In fact, the 
method extends easily to any dimension, showing that in dimension D, there are no further 
terms h{d'^(f))"' for n > D. 

Define a new tensor, 

]C^{9, H) = n - - = E d^iH^^^, ^- = - il-2^t^.)2An - (9-17) 
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Here indices are raised using the full metric g^^t,, H^^ = g^'^Hau, and {H'^Y^ = H^^_^H"^^ ■ ■ ■ H""^-^ 
denotes the matrix product of n tensors. The tensor = g^j^a^^ is defined so that, given 
the Stiikelberg replacement with no vectors if^j/ = h^j^^ + 2 9^(9i^0 — dfj,da(t>di^d°'(j), we have 

IC^Ag,H) =d,d,<j). (9.18) 

Using /C as a new variable, the most general potential in (6.40) can be written as 

Wig, JC) = (/C^) - {icy + £i(/C3) + C2{lC){K) + c^{lCf + ■■■ (9.19) 

where brackets mean traces with respect to the full metric, and the tilde coefficients are 
arbitrary and can be related to those in (6.40) by expanding. Because of the property (9.18), 
this reorganization of the potential makes it easy to see what the scalar self-interactions look 
like - they are simply W{g,lC) = W{g,Il). Thus the A3 theory corresponds to choosing 
coefficients in (9.19) such that the /C terms appear in the total derivative combinations of 
Appendix A. Each total derivative combination can have an arbitrary overall coefficient, so 
the A3 theory corresponds to 

Ty(^7,/C) = 5^a„£™(/C), (9.20) 

n>2 

with arbitrary coefficients These coefficients correspond to the free coefficients of the A3 
theory, 0^2 = —2^, 0:3 = 2^C3, 04 = 2'^d5, etc. In (9.20), the sum is finite and stops ai n = D, 
since the total derivative combinations vanish for n > D. 

The decoupling limit interactions contain only one power of h, so the entire decoupling 
limit action is given by (note that in this section the fields are not canonically normalized) 

S = Jd^^^ {^^KuS^^^'^^Kp - "^h^^'X^^j , (9.21) 



where 



Using the relation 



we calculate 



^(^'^)L..^o = i(np^-ny, (9-23) 



(9.24) 
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where we have used the definitions 11 = d^d^cj), as well as 11°^ = rj^^, and 11"^^ = 0, and the 
xjn^^ are the identically conserved combinations of d^dycj) described in Appendix A. Thus we 
have 

^'^ = \Y. W"^ + ^4""'^) • (9-25) 



2 

n>2 



For D = A this agrees with (9.10), showing that (9.10) contains all the scalar and tensor 
terms of the decoupling limit. Some other re-summations are discussed in [147, 93]. 



9.3 The appearance of galileons and the absence of ghosts 

We can partially diagonalize the interaction terms in (9.10) by using the properties (A. 18). 
First, we perform the conformal transformation needed to diagonalize the linear terms, h^i^ — ?■ 
h^u + after which the lagrangian takes the form 

(9.26) 

Here the brackets are traces of Il^jy = dud^n and its powers (the notation is explained at the 
end of the Introduction). 

The cubic hcpcf) couplings can be eliminated with a field redefinition h^u — ^ h^^y + 
2(6c3-i) g^^g^^^ after which the lagrangian reads, 

40(6c3-l)(84 + C3: 



^^^^±^(90)^ ([n]^ - 3[n^][n] + 2\ti' 



(9.27) 

There is no local field redefinition that can eliminate the h(f)(j)(f) quartic mixing (there is a 
non-local redefinition that can do it), so this is as unmixed as the lagrangian can get while 
staying local. 
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The scalar self-interactions in (9.27) are given by the following four lagrangians, 

A = -\{dct>f - mm + 2[n3]) . (9.28) 

These are known as the galileon terms [136] (see also Section II of [137] for a summary of the 
galileons). They share two special properties: their equations of motion are purely second 
order (despite the appearance of higher derivative terms in the lagrangians), and they are 
invariant up to a total derivative under the galilean symmetry (8.8), 0(a;) — )■ 0(a;) + c + h^x^ . 
As shown in [136], the terms (9.28) are the only polynomial terms in four dimensions with 
these properties. 

The galileon was first discovered in studies of the DGP brane world model [35] (which 
we will explore in more detail in Section 10.2), for which the cubic galileon, £3, was found 
to describe the leading interactions of the brane bending mode [148, 149]. The rest of the 
galileons were then discovered in [136], by abstracting the properties of the cubic term away 
from DGP. They have some other very interesting properties, such as a non-renormalization 
theorem (see e.g. Section VI of [137]), and a connection to the Lovelock invariants through 
brane embedding [150]. Due to these unexpected and interesting properties, they have since 
taken on a life of their own. They have been generalized in many directions [151, 152, 153, 
154, 155, 156, 157], and are the subject of much recent activity (see for instance the > 100 
papers citing [136]). 

The fact that the equations are second order ensures that, unlike (8.10), no extra 
degrees of freedom propagate. In fact, as pointed out in [34], the properties (A. 17) of the 
tensors X^^ guarantee that there are no ghosts in the lagrangian (9.10) of the decoupling 
limit theory.^^ By going through a hamiltonian analysis similar to that of Section 2.1, we 
can see that /iqo and /loi remain Lagrange multipliers enforcing first class constraints (as they 
should since the lagrangian (9.10) is gauge invariant. In addition, the equations of motion 

This is contrary to [101], which claims that a ghost is stiU present at quartic order. As remarked however 
in [34] , they arrive at the incorrect decouphng hmit lagrangian, which can be traced to a minus sign mistake 
in their Equation 5, which should be as in (9.4). 
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remain second order, so the decoupling limit lagrangian (9.10) is free of the Boulware-Deser 
ghost and propagates 3 degrees of freedom around any background. 

Once the two degrees of freedom of the vector are included, and if there are no 
ghosts in the vector part or its interactions, the total number degrees of freedom goes to 5, 
the same as the linear massive graviton. The vector interactions were shown to be ghost free 
at cubic order in [143]. It was shown in [102] that the full theory beyond the decoupling limit, 
including all the fields, is ghost free, up to quartic order in the fields. This guarantees that 
any ghost must carry a mass scale larger than A3 and hence can be consistently excluded 
from the quantum theory. Finally, in [103] it was shown using the hamiltonian formalism 
that the full theory, including all modes and to all orders beyond the decoupling limit, carries 
5 degrees of freedom. The A3 theory is therefore free of the Boulware-Deser ghost, around 
any background. This can also been seen in the Stiikelberg language [141]. 



9.4 The A3 Vainshtein radius 

We can now derive the scale at which the linear expansion breaks down around heavy point 
sources in the A3 theory. To linear order around a central source of mass M, the fields still 
have their usual Coulomb form, 

Ih^^l (9.29) 

The non-linear terms in (9.10) or (9.27) are suppressed relative to the linear term by a 
different factor than in the A5 theory, 

~ t:^t^- 9.30 

Non-linearities become important when this factor becomes of order one, which happens at 
the radius 



J3) 
V 



MpJ A3 

This is parametrically larger than the Vainshtein radius found in the A5 theory. 

It is important that the decoupling limit lagrangian was ghost free. To see what would 
go wrong if there were a ghost, expand around some spherical background = $(r) + 
and similarly for h^^. The cubic coupling and quartic couphngs could possibly give fourth 
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order kinetic contributions of the schematic form, respectively, 



i3$(9V)^ (9.32) 



These would correspond to ghosts with r-dependent masses, 



2 Ai At 



or, given that the background fields go like $ ~ jf^-^, 



(9.33) 



'"i.-(r) ~ ^A|r. (f^)'A|rl (9.34) 
Thus the ghost mass sinks below the cutoff A3 at the radius 

As happened in the A5 theory, these radii are parametrically larger than the Vainshtein 
radius. This is a fatal instability which renders the whole non-linear region inaccessible, 

unless we lower the cutoff of the effective theory so that the ghost stays above it, in which 

(3) 

case unknown quantum corrections would also kick in at ~ '"ghost' 

swamping the entire non- 
linear Vainshtein region. 



9.5 The Vainshtein mechanism in the A3 theory 

In the A5 theory, the key to the resolution of the vDVZ discontinuity and recovery of GR 
was the activation of the Boulware-Deser ghost, which cancelled the force due to the lon- 
gitudinal mode. In the A3 theory, there is no ghost (at least in the decoupling limit), so 
there must be some other method by which the scalar screens itself to restore continuity 
with general relativity. This method uses non-linearities to enlarge the kinetic terms of the 
scalar, rendering its couplings small. 

To see how this works, consider the lagrangian in the form (9.26). Set = — C3/8, 
C3 = 5/36 to simplify coefficients, and ignore for a second the cubic hcpcj) coupling, so that 
we only have a cubic (p self- interaction governed by the galileon term £3, 

S = Jd'x - 3{d4>)' - ^(50)'n0 + ^<^T. (9.36) 
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This is the same lagrangian studied in [149] in the DGP context. 

Consider the static spherically symmetric solution, 0(r), around a point source of mass 
M, T ~ M6^{r). The solution transitions, at the Vainshtein radius ry"* = (^-^^^ x^, 
between a linear and non-linear regime. For r ^ the kinetic term in (9.36) dominates 
over the cubic term, linearities are unimportant, and we get the usual 1/r Coulomb behavior. 
For r ^ ry \ the cubic term is dominant, and we get a non-linear ^/r potential. 



2 / \ 

a3»,(3) ( \ ^ (3) 

^^3'\/ \ i-i) \ ' ^ 'y ) 

,3^(3)^/ 4" '\ .^.(3) 



{r)-{ \li (9.37) 



We can see the Vainshtein mechanism at work already by calculating the ratio of the 
fifth force due to the scalar to the force from ordinary newtonian gravity, 

3/2 



^\r)/Mp J-liayl r^r^^\ 



i^Newton M / {Ml,r^^ 



V 



1 (3) 

1 r ^ Ty . 



(9.38) 



There is a gravitational strength fifth force at distances much farther than the Vainshtein 
radius, but the force is suppressed at distances smaller than the Vainshtein radius. 

This suppression extends to all scalar interactions in the presence of the source. To 
see how this comes about, we study perturbations around a given background solution $(a;). 
Expanding 

= $ + V?, T = To + ST, (9.39) 

we have after using the identity {d'^ip)nLp = dy [d'^^pd^ip — ^r]'^'^ {d^py] on the quadratic parts 
and integrating by parts 

S^ = j d^x - Sidipf + 1^ (d^dy^ - ri^yD^) d'^ipd-'ip - j-^{d^yn^ + ^¥^5r. (9.40) 

Note that expanding the cubic term yields new contributions to the kinetic terms, with 
coefficients that depend on the background. Unlike the A5 lagrangian (8.10), no higher 
derivative kinetic terms are generated, so no extra degrees of freedom are propagated on any 
background. This is a property shared by all the galileon lagrangians (9.28) [158]. 

Around the solution (9.37), the coefficient of the kinetic term in (9.40) is 0{1) at 



(3) I T'^'^^ \ (3) 

distances r ^ Vy , but goes like I -y- J for distances r ^ Vy . Thus the kinetic term 
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3/2 



is enhanced at distances below the Vainshtein radius, which means that after canonical 
normalization the couplings of the fluctuations to the source are reduced. The fluctuations 
(y? effectively decouple near a large source, so the scalar force between two small test particles 
in the presence of a large source is reduced, and continuity with GR is restored. A more 
careful study of the Vainshtein screening in the A3 theory, including numerical solutions of 
the decoupling limit action, can be found in [145]. 



9.6 Quantum corrections in the A3 theory 

As in Section 8.4, we expect quantum mechanically the presence of all operators with at 
least two derivatives per 0, now suppressed by the cutoff A3 (we ignore for simphcity the 
scalar tensor interactions). 



(9.4i: 



These are in addition to the classical galileon terms in (9.27), which have fewer derivatives 
per (j), and are of the form 



(9.42) 



An analysis just like that of Section 8.4 shows that the terms (9.41) become important 
relative to the kinetic term at the radius r ~ f -^^^ j This is the same radius at which 

classical non-linear effects due to (9.42) become important and alter the solution from its 
Coulomb form. Thus we must instead compare the terms (9.41) to the classical non-linear 
galileon terms (9.42). We see that the terms (9.41) are all suppressed relative to the galileon 
terms (9.42) by powers of d/A-^, which is ~ ^ regardless of the non-linear solution. Thus, 
quantum effects do not become important until the radius 

~ (9-43) 

which is parametrically smaller than the Vainshtein radius (9.31). 

This behavior is much improved from that of the A5 theory, in which the Vainshtein 
region was swamped by quantum correction. Here, there is a parametrically large intermedi- 
ate classical region in which non-linearities are important but quantum effects are not, and 
in which the Vainshtein mechanism should screen the extra scalar. In this region, GR should 
be a good approximation. See Figure (4). 



95 



Quantum 
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Non-linear Linear 



r — 

rq ~ 10'^ km ry ~ 10^^ km 



Figure 4: Regimes for massive gravity with cutoff A3 = (Mpm^)^/^, and some values within the 
solar system. The values are much more reasonable than those of the A5 theory. 

As in the A5 theory, quantum corrections are generically expected to ruin the various 
classical tunings for the coefficients, but the tunings are still technically natural because the 
corrections are parametrically small. For example, cutting off loops by A3, we generate the 
operator ~ ^(□<^)^, which corrects the mass term. The canonically normalized is related 
to the original dimensionless metric by ~ -^dd(j), so the generated term corresponds in 

unitary gauge to Ag/i^ = M^iv? i^jf^ h^-, representing a mass correction 5m? ~ m? i^jt)- 
This mass correction is parametrically smaller than the mass itself and so the hierarchy 
m <^ A3 is technically natural. This correction also ruins the Fierz-Pauli tuning, but the 
pathology associated with the de-tuning of Fierz-Pauli, the ghost mass, is ~ srnf/m? ~ ^3' 
safely at the cutoff. 

We should mention another potential issue with the A3 theory. It was found in [f 36] that 
lagrangians of the galileon type inevitably have superluminal propagation around spherical 
background solutions. No matter what the choice of parameters in the lagrangian, if the 
solution is stable, then superluminality is always present at distances far enough from the 
source (see also [159]). It has been argued that such superluminality is a sign that the 
theory cannot be UV completed by a standard local Lorentz invariant theory [160], though 
others have argued that this is not a problem [161]. In addition, the analysis of [136] was 
for pure galileons only, and the scalar-tensor couplings of the massive gravity lagrangian can 
potentially change the story. These issues have been studied within massive gravity in [162]. 
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10 Massive gravity from extra dimensions 



So far, we have stuck to the effective field theorist's philosophy. We have explored the 
possibility of a massive graviton by simply writing down the most general mass term a 
graviton may have, remaining agnostic as to its origin. However, it is important to ask 
whether such a mass term has a top down construction or embedding into a wider structure, 
one which would determine the coefficients of all the various interactions. This goes back 
to the question of whether it is possible to UV complete (or UV extend) the effective field 
theory of a massive graviton. 

One way in which a massive graviton naturally arises is from higher dimensions. We 
will now study several of these higher dimensional scenarios, the Kaluza-Klein reduction, 
the DGP brane world model, and a model of a non-dynamical auxiliary extra dimension, 
showing in each case how massive gravitons emerge in a 4d description. 

10.1 Kaluza-Klein theory 

In the original Kaluza-Klein idea [163, 164] (see [165] for a review), gravity on a 5d space with 
a single compact direction is dimensionally reduced onto the four non-compact dimensions, 
where it is found that the lightest modes describe Einstein gravity, electromagnetism, and a 
massless scalar, all in interaction with each other. In almost all work on Kaluza-Klein theory 
(including the rather large subset going by the name of string compactifications) , only the 
lowest energy modes are considered. 

Beyond the lowest energy modes, there is an entire tower of massive fields. In the 
dimensional reduction of gravity, this tower will consist of massive gravitons. We will now 
review the dimensional reduction of pure 5d gravity down to four dimensions. We will 
work at the linear level, keeping all the massive modes, and we will see that the massive 
gravitons which arise are described by the 4d part of the 5d metric obeying precisely the 
Fierz-Pauli mass term (2.1). The Fierz-Pauli tuning of coefficients arises automatically from 
the dimensional reduction. In addition, the 5d components of the 5d metric become a tower 
of 4d scalars and a tower of 4d vectors. 

There is also the question of gauge symmetry. The 5d gravity action has 5d diffeomor- 
phism invariance. The result of the reduction, a tower of massive gravitons in 4d, has no 
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diffeomorphism symmetry, so where does this symmetry go? We will see that what comes 
out of the reduction is not the unitary gauge action (2.1), but rather the Stiikelberg-ed action 
(5.11). The 5d gauge symmetry becomes the 4d Stiikelberg gauge symmetry, and the towers 
of vectors and scalars become the Stiikelberg fields. 



We start with a massless graviton in 5 dimensions, with 5d Planck mass M> 



5; 

S = Ml j d'X-]^dcHABd''H^''+dAHBcd''H^''-dAH^''dBH+]^dAHd^H+^^HABT^^. 

(10.1) 

Here Has is the dimensionless 5d graviton, with indices A,B,... running over 5d spacetime. 
We divide spacetime into 4d coordinates x^, and a fifth coordinate y, so that X"^ = (x^,y). 
We compactify y so that it runs along a circle of circumference L, y G (0,/^). T^^ is the 
fixed external 5d stress tensor, which is conserved in 5d, dsT^^ = 0. 

Now we change variables by expanding in a Fourier series over the circle, 

oo 



71= — oo 
OO 



H^,yix,y) = ^ A^,„(x)e*' 

n=— oo 
oo 

Hyy{x,y) = (10-2) 

n=— oo 

Here n is an integer, oon = and we have the usual orthogonality relation dy (e*'^™*')* e^^"^ 

The coefficients in the Fourier expansion, 0„, and v4^^„, are the new variables 

and will become the 4d fields. Reality of the 5d fields imposes the conditions 

^liu,n = h^u-n-, A*^^n = ^t^-n, 4>n = 4>-n- (10.3) 

In addition, we decompose the 5d stress tensor in similar fashion, 

oo 

T^,(x,2/) = Y (10.4) 



n=— oo 
oo 



T,y{x,y) = ^M,n(x)e^""^ (10.5) 

n=— oo 
oo 

Tyy{x,y) = Y Jn{x)e'^"'. (10.6) 
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The fields t^^^n, j^i,n and j„, which satisfy reality conditions just like (10.3), 

will become the 4d sources. The equation for 5d stress tensor conservation, d^T^B = 0, 
when expanded out in components and in the Fourier series, implies 

9%,o = 0, d%,fl = 0, (10.8) 

J>.,n = —d%.,n, Jn = -^5^5V,n, n ^ 0. (10.9) 

Plugging the Fourier expansions into (10.1) and doing the y integral, we get the fol- 
lowing equivalent 4d action, 

S = LMl j d\ ^ Vo^^'^'^^/^afto - ^i^i.o + C {d^d,<P^ - r]^M 

12 1 

oo 

I \ ^ I,* c^iu,aPi _ , t"^ (\h |2 _ I /, |2\ _ I p |2 

~r / J '''fiu,n^ i''a/3,n \\ fJ-'^^'n-l \''"n\ J {-'^ fj,i',n\ 

n=l 

+ [2iUnA^^n {duhr - d^K) + K (d^d^K ' V^..O<P*J + c.c] 

n=l L 5 5 5 

(10.10) 

We have used the reahty conditions (10.3) and (10.7) to change the range of the sum. This 
action is exactly equivalent to (10.1), and describes all the 5d dynamics. We have not 
truncated anything or restricted the fields in any way, we have merely changed variables to 
ones that are more easily recognizable in 4d. 

From the prefactor we can read off the effective 4d Planck mass 

= LMl (10.11) 

We now study the fate of the gauge symmetry. The 5d action has the gauge symmetry 
5Hab = dA^B + ds'^A, for a gauge vector S^(X). Fourier decomposing the gauge parameter, 

oo 

71= — oo 
OO 

Ey{x,y) = Yl (10-13) 
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where the coefficients have reahty properties hke those of (10.3). The gauge transformations 
can be decomposed component by component to yield 

= + du^f^^n, (10.14) 

6A^^n = df,^n + iUJn^i,,n, (10.15) 
S(pf,u,n = ^iUn^n- (10.16) 

The zero mode of the 5d gauge parameter breaks up into a vector and a scalar, which 
become the linear diffeomorphism invariance o of the zero mode graviton h^j^^^Q, and the 
Maxwell gauge invariance of the zero mode vector A^^o (the zero mode scalar 0o is gauge in- 
variant). The n ^ modes, on the other hand, get transformations of exactly the Stiikelberg 
form (4.22). The 5d gauge symmetry has become the 4d Stiikelberg symmetry. 

In fact, the action (10.10) can be written solely in terms of the following gauge invariant 
combination for n 0, 

i 1 
h^,u,n H {d^,A^,n + duAf^^n) l^df,dy(t)n, (10.17) 

which is just the linear Stiikelberg replacement rule. The action (10.10) for the 7^ modes 
is precisely the complex version of our Stiikelberg action (5.11) for a massive graviton. The 
higher modes of the /i5 and 55 components of the 5d metric Hab have become the non- 
physical Stiikelberg fields, and are pure gauge. 

Fixing the unitary gauge A^ n = 0n = for n 7^ 0, and canonically normalizing, we 

have 



12 1 

H hiiyntn" H A,, n?'n H </'n7n 

00 

n=l 



(10.18) 



which shows that the theory (after the conformal transformation (4.24) for the zero modes), 
consists of a single real massless graviton, a single real massless vector, a single real massless 
scalar, and a tower of complex massive Fierz-Pauli gravitons with masses 

27rTi 

mn = u}n = ^—, n = 1,2,3,- ■■. (10.19) 
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These fields are all coupled to the various modes and components of the 5d stress tensor. In 
addition, it can be easily shown that the translation symmetry along y in the original theory 
becomes an internal U{1) rotating the phase of the massive gravitons. There are interesting 
issues that arise when one wishes to couple this U{1) to electromagnetism in the case of a 
single massive graviton [166]. 

To go beyond linear order, we would put the higher order in H interactions coming 
from the 5d Einstein-Hilbert term into the 5d action 10.1, make the same change of variables 
into Fourier components (10.2), then plug in and do the y integral. This will give a slew 
of interaction terms in 4d, involving all the modes interacting with each other. This should 
be a consistent, stable, ghost free theory of an infinite number of fully interacting massive 
gravitons, since it is equivalent to 5d Einstein gravity which we know to be consistent. There 
should be no strong coupling problems or low scale cutoffs, and the effective theory should be 
valid all the way up to the 5d Planck mass. All the 4d graviton modes should miraculously 
interact in such a way as to cancel out all the strong coupling effects we have uncovered for 
a single massive graviton [104]. It is possible to write these interactions for all the fields to 
all orders in closed form [167, 168, 169]. 

It turns out to be consistent to truncate the theory to only the zero modes (consistent 
in the sense that the processes of truncating and deriving the equations of motion commute). 
This leaves 4d Einstein-Hilbert self-interactions for the zero mode graviton, in addition to 
other interactions between the various zero mode fields. The ansatz that describes this 
truncation, which was the original Kaluza-Klein ansatz, is to take the 5d metric to be 
independent of y. It would be desirable to find a consistent truncation that involves only a 
single massive graviton (or a finite number) so that we could study the resulting consistent 
interactions. This does not appear to be possible in this simple model, but may be possible 
in compactifications involving more complicated manifolds or sets of fields. 

10.2 DGP and the resonance graviton 

The Dvali-Gabadadze-Porrati (DGP) model [35] is an extra-dimensional model which has 
spawned a great deal of interest (see the > 1300 papers citing [35]). It provides another, 
more novel realization of a graviton mass. Unlike the Kaluza-Klein scenario, in DGP the 
extra dimensions can be infinite in extent, though there must be a brane on which to confine 
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standard model matter (see [170] for lectures on large extra dimensions). By integrating out 
the extra dimensions, we can write an effective 4d action for this scenario which contains a 
momentum dependent mass term for the graviton. This provides an example of a graviton 
resonance, i.e. a continuum of massive gravitons. 

Another model which has revived a great deal of attention (> 4800 citations) is the 
Randall- Sundrum brane world [171], in which there is a brane floating in large warped extra 
dimensions. This model is not as interesting from the point of view of massive gravity at 
low energies, since the 4d spectrum is similar to ordinary Kaluza-Klein theory, containing 
ordinary Einstein gravity as a zero mode, and then massive gravitons as higher Kaluza-Klein 
modes. See [172] for a review on brane world gravity. 



The DGP action 

DGP is the model of a 3 + 1 dimensional brane (the 3-brane) floating in a 4 + 1 dimensional 
bulk spacetime. Gravity is dynamical in the bulk and the brane position is dynamical as 
well, and the action contains both 4d and 5d parts, 

(fX^^R{G) + ^ / d^x^R{g) + Sm- (10.20) 

Here X^, A, 5, ■ ■ ■ = 0, 1, 2, 3, 5 are the 5d bulk coordinates, Gab{X) is the 5d metric, and 
M5 is the 5d Planck mass, x^, /i, i/, . . . = 0, 1, 2, 3 are the 4d brane coordinates, g^u{x) is the 
4d metric which is given by inducing the 5d metric Gab onto the brane, and M4 is the 4d 
Planck mass. Sm is the matter action, which we imagine to be localized to the brane, 

Sm = j d'xCM{g,i^), (10.21) 

where ip{x) are the 4d matter fields. Due to the presence of a brane Einstein-Hilbert term, 
this scenario is also called brane induced gravity [173]. 

The dynamical variables are the 5d metric depending on the 5d coordinates, the em- 
bedding X"^{x) of the brane depending on the 4d coordinates, and the 4d matter fields 
depending on the 4d coordinates, 

Gab{X), X^{x), ^{x). (10.22) 

The 4d metric is not independent, but is fixed to be the pullback of the 5d metric, 

g^,{x) = d^X^d.X^GAB {X{x)) . (10.23) 
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Note that the dependence of the action on the enters only through the induced metric 

The action (10.20) has a lot of gauge symmetry. First, there are the reparametrizations 
of the brane given by infinitesimal vector fields ^'^(x), under which the X"^ transform as 
scalars and the matter fields transform as tensors (i.e. with a Lie derivative), 

5^X^ = ^f'd^X^, 5^tl) = C^ip. (10.24) 

Second, there are reparametrizations of the bulk given by infinitesimal vector fields H'^(X), 
under which Gab transforms as a tensor and the X^ shift, 

S^Gab = Va^b + Vb^a, feX^ = -H^(X). (10.25) 

The induced metric g^,y transforms as a tensor under 6^, and is invariant under 

We first proceed to fix some of this gauge symmetry. In particular, we will freeze 
the position of the brane. Note that the brane coordinate functions, X"^{x), are essentially 
Goldstone bosons since they shift under the bulk gauge symmetry, X^{x) — )■ X^{x) — 
S'^(X(x)). We can thus reach a sort of unitary gauge where the X^ are fixed to some 
specified values. We will set values so that the brane is the surface X^ = 0, and the brane 
coordinates x^ coincide with the coordinates X^, thus we set 

X'^(x) = x>', /i = 0,1,2,3, (10.28) 
X\x) = 0. (10.29) 

There are still residual gauge symmetries which leave this gauge choice invariant. Act- 
ing with the two gauge transformations 6^ and 5= on the gauge conditions and demanding 

^^To see invariance under S-e, transform 

SsGab - Ss{d^X^d,X''GAB{X{x))) 

= -d^E^d.X^'GAB iX{x)) - d^X^d^E^'GAB {X{x)) + d^X^d.X'' SeGab {X{x)) (10.26) 

then in transforming Gab, remember that both the function and the argument are changing, 

SeGab {X{x)) = CeGab {X{x)) - S^OcGab- (10.27) 

Putting aU this together, we find SeGab — 0. 
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that the change be zero, we find 



(X(x)) = 0. 



E^{X{x))+CduX^{x) ^ 

Xf(a;)=: 



X-^(x)=0 



(X(x)) 



E^{X{x))+e{x) 



(10.30) 



(X(x)) =e{x) 



(10.31) 



The residual gauge transformations are bulk gauge transformations that do not move points 
onto or off of the brane, but only move brane points to other brane points. Furthermore, 
the brane diffeomorphism invariance is no longer an independent invariance but is fixed to 
be the diffeomorphisms induced from the bulk. 

We now fix this gauge in the action (10.29), which is permissible since no equations of 
motion are lost. This means that the induced metric is now 



We split the action into two regions, region L to the left of the brane, and region R to 
the right of the brane, with outward pointing normals, as in Figure 5. We call the fifth 
coordinate X^ = y. The brane is at y = 0. 



0). 



(10.32) 



y 



L 




n' 



^J■ 

L 



R 



Figure 5: Splitting the DGP action 
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(fxdyy/^R{G) + j (fx £4 



[10.33) 



where £4 = ^\/—gR{g) + £m(5',V') is the 4d part of the lagrangian. To have a well 
defined variational principle, we must have Gibbons-Hawking terms on both sides [100], 
corresponding to the outward pointing normals. Adding these, the resulting action is 

c _ Ml 



( [ + [] d^xdyy/^R{G) + 2 (f d^x^^KL + 2 (f d^xy/^Kn 

\Jl JrJ Jl Jr 



+ / d'^x £4 , 



(10.34) 

where Kn, Kl are the extrinsic curvatures relative to the normals and ni respectively. 

We now go to spacelike ADM variables [97, 98] adapted to the brane (see [99, 100] 
for detailed derivations and formulae). The lapse and shift relative to y are N'^{x,y) and 
N{x,y), and the 4d metric is g^u{x,y). The 5d metric is 



Gab 



(10.35) 



The 4d extrinsic curvature is taken with respect to the positive pointing normal til, and is 
given by 



K„ 



2N 

where a prime means derivative with respect to y. The action is now^'^ 



(10.36) 



S 



Ml 



+ I d'^x £4. 



[10.39) 



It can be checked that a fiat brane living in fiat space is a solution to the equations 
of motion of this action. This is called the normal branch. There is another maximally 



^''The Ricci scalar and metric determinant are 



(10.37) 
(10.38) 



The total derivatives coming from 2\7a (n^Vsn'^ — n^K^ in the Einstein-Hilbert part of the action exactly 
cancel the Gibbons-Hawking terms. 
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symmetric solution with a flat 5d bulk, which contains a de Sitter brane with a 4d Hubble 
scale H ~ M^/M^. This is called the self-accelerating branch, and has caused much interest 
because the solution exists even though the brane and bulk cosmological constants vanish. 

Linear expansion 

To see the particle content of DGP, we will expand the action (10.39) to linear order around 
the flat space solution, and then integrate out the bulk to obtain an effective 4d action. We 
start by expanding the 5d graviton about flat space 

Gab = VAB + Hab. (10.40) 

We use the lapse, shift and 4d metric variables, with their expansions around flat space, 

9i,u = rj^y + hf,y, N^, = n^,, N = I + n. (10.41) 

We have the relations, to linear order in h^^,n^,n, 

H^, = h^„ H^5 = n^, H55 = 2n. (10.42) 

We will first expand the DGP action (10.39) to quadratic order in h^i^,n^,n. We will then 
solve the 5d equations of motion, subject to arbitrary boundary values on the brane and 
going to zero at infinity. We then plug this solution back into the action to obtain an 
effective 4d theory for the arbitrary brane boundary values. 

The 5d equations of motion away from the brane are simply the vacuum Einstein 
equations, which read, to linear order, 

- 2RAB{G)uno.r = ^^'^Hab + OaObH - d'^dAHnc - d^'dBHAc = 0. (10.43) 
We will solve (10.43) in the de Bonder gauge, 

d^'HAB - \dAH = 0. (10.44) 

With this, (10.43) is equivalent to 

U^^'^Hab = 0, (10.45) 
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along with the de Bonder gauge condition (10.44). In terms of the ADM variables, (10.45) 
becomes 



□V + 5' V = 0, (10.46) 

52„ 



□^M + = 0. (10-47) 

Un + din = 0, (10.48) 

where □ is the 4d laplacian. These have the following solutions in terms of boundary values 
h^^{x),nf,{x),n{x), 

Ku{x,y) = e-^^V(a^), (10.49) 

n^ix,y) = e-^^n^(x), (10.50) 

n{x,y) = e~^^n{x). (10.51) 

Here, the operator A is the formal square root of the 4d laplacian, 

A = v^^. (10.52) 

The A = n and A = 5 components of the gauge condition (10.44) are, respectively, 

d''h^^-^d^h + dyn^-d^n = 0, (10.53) 

d^n^ - ^dyh + dyu = 0. (10.54) 

For these to be satisfied everywhere, it is necessary and sufficient that the boundary fields 
satisfy the following at y = 0, 

d^h^^ - ^d^h - An^ - d^n = 0, 

d'^Uf, + ^Ah - An = 0. (10.55) 



These should be thought of as constraints determining some of the boundary variables in 

terms of the others^"*. We will at this point imagine that we have solved these constraints, 

and that the action is really a function of the independent variables. 

^^Note that we cannot think of them as determining n'', n in terms of /i^^. Acting with 9^ on the first 
equation, A on the second, and then adding, we find the equation 

df.d^hi"' - D/i = 0, (10.56) 

which is precisely the statement that the 4d hnearized curvature vanishes (which is, in turn, the hnearized 
hamiltonian constraint in general relativity). Thus, we must think of these constraints as determining some 
of the components of the metric. 



107 



The deDonder gauge is preserved by any 5d gauge transformation S"^ satisfying 

□ (^)H^ = 0. (10.57) 

The component must vanish at y = because the position of the brane is fixed. (10.57) 
then imphes that vanishes everywhere. The other components can have arbitrary values 
S^(x, 0) = ^^{x) on the brane, which are then extended into bulk in order to satisfy (10.57), 

E>^ix,y) = e-y^a^). (10.58) 

The residual gauge transformations acting on the boundary fields are then 

Sh^^ = d^^y + d^iu, 

Sn = 0. (10.59) 

The constraints (10.55) are invariant under these gauge transformations. The 4d effective 
action must and will be invariant under (10.59). 

The 5d part of the action reads 

^5 = ^ y d'^xdy N^g [R{g) + - K^.K^''] . (10.60) 

We want to expand this to quadratic order in h^y,n^j_,n and then plug in our solution. We 
will need the expansion of K^j_i, to first order, 

K^y = ^ {dyhf^i, - d^riy - d^n^) . (10.61) 

Expanding, we have (after much integration by parts in 4d) 

= / d^xdy nd^dyh^" - nOh + ^dxh^^d^h^^ - ^d^hd.h^" 

Now, in the last line, integrate by parts in y, picking up a boundary term at y = 0, and use 
(10.45) to kill the bulk part, 

Jj55 = I d'xdy nd^dyh^'' - nUh + V^'^'^" " \dMuh^'' 



-dyhdX + \{d^n^'f + dyh^.ud^n" + ]^n^Un^ 



+ I d^x -hidh+hif.ydyh^'r 
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We now insert the following term into the action 



^GF = - Y j [d^'HAB - ^dAHj . (10.62) 

The 5d equations of motion solve the de Bonder gauge condition, so this term contributes 
to the action (thought of as a function of the unconstrained variables) and we are free to 
add it. However, we write it in terms of the unconstrained 4d variables for now, 

-^'S'gf = / d'^xdy (d^h^^ - ]^d^h + dyU^ - d^^iA (10.63) 



1/. „ 1 



^5 

-^^[dX-^i^dyh + dynj . (10.64) 

Adding this to the previous 5d term, we find that after using the 5d Laplace equations, the 
entire action can be reduced to a boundary term at ?/ = 0, 

^ (^5 + Sgf) = I d'x --^K^Ah^' + ^-hAh - \iAn - \n^An^ 

5 'J 



Now a crucial point. We have been imagining solving the constraints (10.55) for the 
independent variables. But now, consider the action (10.65) as a function of the original 
variables h^^,, n^, n. Varying with respect to and ra, we recover precisely the constraints 
(10.55). Thus, we can re-introduce the solved variables as auxiliary fields, since the con- 
straints are then implied. The action now becomes a function of h^j,y,n^,n. 

Now add in the 4d part of the action. 



2 



J d^x V^Rig) + Sm + 2 (S, + Sgf) ■ (10.66) 

where Sm is the 4d matter action and the factor of 2 in front of the 5d parts results from 
taking into account both sides of the bulk (through boundary conditions at infinity we have 
thus implicitly imposed a Z2 symmetry). 

Expanded to quadratic order. 



1 1 



4 

+hAn + (-28^71 - d^h + 2d''h 



huuAh'"' + -hAh - nAn - n^An^" 
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where S'^'^'"'^ is the massless graviton kinetic operator (2.46), and 



m ^ ^ (10.68) 



is known as the DGP scale. 



The action (10.67) should now be compared to the massive gravity action in the form 
(4.44). It is invariant under the gauge transformations (10.59), under which plays the 
role of the vector Stiikelberg field, n plays the role of the gauge invariant auxiliary field. To 
get this into Fierz-Pauli form, first eliminate n as an auxiliary field by using its equation of 
motion. Then use (10.59) to fix the gauge = 0. The resulting action is 



2 r 



1 



-h.^S^^'^'^h^p - -m {h.^Ah'^'' - hAh) 



1 



+ ^h^,T^\ (10.69) 

which is of the Fierz-Pauli form, with an operator dependent mass term niA. 

One can go on to study interactions terms for DGP, and the longitudinal mode turns 
out to be governed by interactions which include the cubic galileon term ~ (90)^00 [148, 
149, 174], and are suppressed by the scale A3 = {Mim?Y/'^ (in fact this was where the 
galileons were first uncovered). In this sense, DGP is analogous to the nicer A3 theories of 
Section (9). The theory is free of ghosts and instabilities around solutions connected to flat 
space [149], but changing the asymptotics to the self-accelerating de Sitter brane solutions 
flips the sign of the kinetic term of the longitudinal mode, so there is a massless ghost around 
the self-accelerating branch [175]. This branch is completely unstable, which is bad news 
for doing cosmology on this branch. In addition to ghosts, there are other issues with other 
non-trivial branches, such as superluminal fluctuations [176], and uncontrolled singularities 
and tunneling [177]. 



Resonance gravitons 

The operator dependent mass term in (10.69) is known as a resonance mass, or soft mass 
[178, 179, 180, 181, 182, 183, 184, 134, 52, 55]. To see the particle content of this theory, 
we will decompose the propagator into a sum of massive gravity propagators. The linear 
Stiikelberg analysis, leading to the propagators (4.37), goes though identically, with the 
replacement — )■ rriA. The momentum part of the propagators now reads 

, ~' r^ ^ (10.70) 
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Setting z = —p'^, the propagator has a branch cut in the z plane from (0, oo), with disconti- 
nuity 

2m 

(10.71) 



A branch cut can be thought of as a string of simple poles, in the limit where the spacing 
between the poles and their residues both go to zero. The function f{z) = J^^dXp{X)^^ 
has a cut along the real axis everywhere that p is non-zero, with discontinuity —2inp{z). We 
can see this by noting 

/oo -1 -1 /'OO 
dXpiX) ^ ^ =/ dXp{X)\-27nSiz- X)]. 
z-X + ie z-X-ie 

Using all this, and the fact that analytic functions are determined by their poles and 
cuts, we can write the propagator in the spectral form 

°° -i m 
ds^— p(s), p(s) = ^— ^>0. (10.72) 



'p^^m^ Jo p^ + s ' 7r^/s{s + m'^) 

The spectral function is greater than zero, so this theory contains a continuum of ordinary 
(non-ghost, non-tachyon) gravitons, with masses ranging from to oo. This is what would 
be expected from dimensionally reducing a non-compact fifth dimension. The Kaluza-Klein 
tower has collapsed down into a Kaluza-Klein continuum. 

In the limit m — )• 0, where the action becomes purely four dimensional, the spectral 
function reduces to a delta function, 

p{s) -> 2S{s), (10.73) 

and the propagator reduces to —i/p^ representing a single massless graviton, vector and 
scalar, as can be seen from (4.37) (the extra factor of two is taken care of by noting that 
the integral is from to oo, so only half of the delta function actually gets counted). This 
theory therefore contains a vDVZ discontinuity. 

The potential of a point source of mass M sourced by this resonance graviton displays 
an interesting crossover behavior. Looking back at (3.11) with the momentum space replace- 
ment m? — )■ m^/jp', and using the relation = —hoo/M^ for the newtonian potential, we 
have 

-2M r d^p 



3M| J (27r 
2 M 



3 M|27r2r 





-1 




- m p 








© 
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(10.74) 



sin ( — ) ci ( — 1 ~'~ 2 f — J fvr — 2sif — 



M _|_ M 



^(^^) ^ ; 3 A/|47rr ^ 37r2Af| ro L ■ yroyj (10.77) 



M 



where 

si(x) = / — sint, ci(a;) = 7 + Inx + / — (cost — 1), (10.75) 
7 ~ 0.577. . . is the Euler-Masceroni constant, and the length scale Tq is 

ro = — = 10.76 

The potential interpolates between 4d ~ 1/r and 5d ~ 1/r^ behavior at the scale ro, 

7-l + ln(^^)] +C(r), r<ro 

3M|4^ ^ (;^) ' r>ro. 

Physically, we can think of gravity as being confined to the brane out to a distance ~ ro, 
at which point it starts to weaken and leak off the brane, becoming five dimensional. This 
is the behavior that is morally responsible for the self-accelerated solutions seen in DGP 
[21]. It has been suggested that corrections to the newtonian potential for r ^ ro may be 
observable in lunar laser ranging experiments [185, 186]. 

The resonance massive graviton can also be generalized away from DGP, by replacing 
the mass term with an arbitrary function of the laplacian [181, 183, 52], 

-> m^{D). (10.78) 

(See [187] for even further generalizations.) At large distances, where we want modifications 
to occur, the mass term will have a leading Taylor expansion, 

m\n) = L2("-i)n°, (10.79) 

with L a length scale and a a constant. In order to modify newtonian dynamics at large 
scales, (9 ^ the mass term should dominate over the two derivative kinetic terms, so we 
should have a < 1. Additionally, there is the constraint that the spectral function (10.72) 
should be positive definite, so that there are no ghosts. This puts a lower bound a > 
[183]. An analysis like that of Section 4.3 can be repeated with the more general mass term 
(10.79), and it turns out that degravitation can be made to work only for a < 1/2 [52]. 
DGP corresponds to a = 1/2, and so it just barely fails to degravitate, but by extending 
the DGP idea to higher co-dimension [188, 189] or to multi-brane cascading DGP models 
[190, 191, 192], a < 1/2 can be achieved and degravitation made to work [193]. Some N-body 
simulations of de-gravitation and DGP have been done in [194, 195, 196]. 
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10.3 Massive gravity from an auxiliary extra dimension 



There is a way of writing a mass term by using an integral over a non-dynamical auxiliary 
extra dimension [197, 198], and which provides a re-summation of some particular massive 
gravity theories. 

To see what this entails, consider the following perverse way of writing the mass term for 
an ordinary scalar field, (f){x). We let the scalar depend on an additional auxiliary parameter 
u, which takes values in the interval u G (0, 1). We then write the following action for the 
field u), 



S 



1 



1 



u=0 



/ du{du(t)f 



(10.80) 



The first term is a normal kinetic term, with derivatives over the x^ only, and the field there 
is 0(x) = (f){x,u = 0). The second term contains a u derivative only, and an integral over 
all the u values. This can be thought of as an action for the infinite number of fields (j){x, u) 
parametrized by m, where only the u = field has a kinetic term. 

We can get a 4d effective action for this theory by integrating out the u dimension, to 
obtain an equivalent theory for only 0(x). To do this, we solve the equations of motion for 
all the (j){x,u) with -u 7^ in terms of 0(x), and then plug back into the action, i.e. we can 
think of all the (j){x, u) with m 7^ as auxiliary fields and eliminate them via their equations 
of motion. To get a unique solution, we must set boundary conditions at m = 1, which we 
will choose to be a Dirichlet condition (j){x, 1) = 0. 

The equations of motion for -u 7^ and the boundary conditions are then, 
die 



(x,l) 



0, 



0, 0(a:,O) = 0(x 

with the solution 

(t){x, u) = (1 — u)(t>{x). 
Plugging back into the action, we find 

S = d''x (^- ^(50)2 - 

The integral over u has become an ordinary mass term for the scalar. 

We now apply a similar construction to the Einstein-Hilbert action, writing 

Mp ' 



10.81 



(10.82) 



(10.83) 



S = d X 



m 



du 



~9 {kl. - k^) 



(10.84) 
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Here the metric is g^^{x,u), it depends on the auxihary parameter u. k^^ = \dug^v is the 
extrinsic curvature into the auxihary dimensions, and indices in the mass term are raised and 
lowered with gn^{x,u) and its inverse. We choose boundary conditions so that the metric 
becomes flat at the end of the auxiliary dimension, 

gt^A^, 0) = g^,^{x), g^,y{x, 1) = r]^y. (10.85) 

The lagrangian (10.84) is invariant under 4d diffeomorphisms, but they are broken by the 
boundary conditions imposed at n = 1. This will be reflected in the effective 4d metric ob- 
tained from integrating out the extra dimension. The boundary condition metric essentially 
plays the role of the flxed metric used for the purely 4d constructions. 

We can now perturbatively solve the equations of motion in the auxiliary dimension 
and expand out the mass term. We write the metric at m = as g^y = rj^^ + h^^ and work 
in powers of h^^. The full w-dependent metric will have an expansion 

g,,{x, u) = T]^, + i/«(x, u) + H^^Jix, «) + ■■■, (10.86) 

where Hjn){x,u), Hjn){x,u), . . . are the terms containing respective powers of hf^i,. The full 
equations of motion for m 7^ derived from (10.84) are 

du [V^ ik'^" - kgn] = IV^g'" {kip - k') + 2^g {kk^^ - k^^k^ ■ (10.87) 

By plugging (10.86) into (10.87), we can collect like powers and solve order by order. To 
lowest order we have simply du (^H^J — rj^^H^^''^ = with solution 

H^^J{x,u) = il-u)h^,ix). (10.88) 

Plugging back into (10.84) we flnd exactly the Fierz-Pauli mass term at quadratic order, 

Mp 



S = d X 



2 

m 



-gR-'^ {hi -h^ + om) 



10.89) 



The order terms were calculated in [143], and it is found that their coefficients are in 
the right combination for raising the cutoff to A3, corresponding to C3 = 1/4, and it seemed 
likely that this theory was a way of re-summing one of the A3 theories. 

Recently, an exact solution for the equation (10.87) has been found, for an arbitrary 
boundary metric f^u^x) at m = 1 [199]. This solution allows us to flnd the 4d lagrangian 
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exactly, 

^ ~ J '^^•^ ' 

Here L is the matrix L^^, = f^^gxv, where f^^ is the inverse of the boundary condition metric 
at M = 1. This lagrangian shows exphcitly how the boundary condition metric becomes the 
fixed metric needed for the graviton mass term. Taking f^y = rjfj^y, then expanding to 
third order in h^^ = g^^ — rj^y reproduces the results of [143], but unfortunately the tuning 
necessary for the A3 theory is violated at fourth order. Thus this theory has a ghost, and a 
cutoff of A5. 

It turns out, however, that a different choice for the boundary metric /^j, can be made 
to restore the A3 tuning order by order, though it requires f^y to depend on the metric 
at M = if flat space is to be a solution in 4d [200]. They also present an argument to 
explain why with rj^i, boundary conditions the third order tuning comes out right while the 
fourth order and higher do not. Note that the combination fc^^ — k'^ of extrinsic curvatures 
appearing in (10.84) is the same combination that appears in the Gauss-Codazzi expansion of 
the flve dimensional curvature in terms of four dimensional quantities (10.37), which suggests 
that there may be some kind of hidden flve dimensional symmetry at work in this model, 
one which enforces the tuning of coefficients in the Fierz-Pauli terms. It is also possible 
to achieve the correct tuning by adding higher powers of k^j_y into the action (10.84) with 
the right coefficients [198, 200]. Some non-flat cosmological solutions of the auxiliary extra 
dimension model are studied in [197]. Extensions to higher co-dimension radially symmetric 
models, and to models including additional bulk Gauss-Bonnet terms, are considered in [200] . 

11 Massive gravity in three dimensions 

We have focused in this review on massive gravity in four dimensions. Extending to higher 
dimensions is a straightforward exercise that does not reveal much which is conceptually new 
beyond what exists in four dimensions (with some exceptions [201]). On the other hand, 

115 



-gR 



-g [detL] 



-1/2 



2[det L]"^/^cosh 



2\/3 



Tr[(lnL)^ 



^(TrlnL)^ 



(10.90) 



dropping down to three dimensions opens up new possibilities. These possibihties are due to 
pecuharities of three dimensions, and they would take another review to fully cover (see [202] 
for a short review). Here we will briefly mention the two main classes of fully interacting, 
covariant massive theories that have been studied in dimension three. 

11.1 New massive gravity 

The first is new massive gravity (NMG), proposed by Bergshoeff, Hohm, and Townsend 
[203]. NMG relies on adding terms which carry four derivatives of the metric, such as R^, 
R^iyR'^^, to the action in addition to the Einstein-Hilbert term. Since these terms are higher 
derivative, degrees of freedom in addition to the graviton will generally propagate. In four 
dimensions, adding such terms with generic coefficients will lead to the propagation around 
fiat space of a massive spin scalar and a massive spin 2 ghost, in addition to the massless 
graviton [204]. Judicious choices of coefficients exist which will project out the massive spin 
2 ghost or the scalar, but no choice will eliminate the massless graviton. There is no way to 
get an action that contains a healthy massive spin 2 and a healthy massless graviton. If the 
sign of the action is fiipped, the massive spin 2 can be made healthy, but then the massless 
graviton will be a ghost. 

In three dimensions, however, there is a loophole, because 3d Einstein gravity contains 
no propagating degrees of freedom (in three dimensions, the canonical analysis of Section 
2.1 shows that there are 3 first class constraints on the 3 spatial components of the metric 
and their canonical momenta, leaving no degrees of freedom). Thus, we can add curvature 
squared terms in a combination such that there is a massive spin 2 and no scalar, and then 
change the overall sign of the action so that the massive spin is healthy. 



Mp is the 3d Planck mass and m is the mass of the graviton. The specific factor of 3/8 in 
the higher derivative terms kills the spin mode. Note the wrong sign Einstein-Hilbert term 
- the massless spin 2 would be a ghost, but it does not matter because in three dimensions 
it carries no degrees of freedom anyway. The result is a theory which propagates only a 
massive spin 2 around fiat space, carrying two degrees of freedom^^. In contrast to the other 

^'^Recall that in three dimensions, spin for massive particles is hke hehcity for massless particles in four 
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massive gravity theories studied in this review, NMG requires no fixed metric in the action, 
so hke Einstein gravity, the theory is free of prior geometry. 

This theory has been recently analyzed with the Stiikelberg method [205]. Like the 
A3 theory in four dimensions, NMG theory carries a scale higher than that which would be 
generically expected (though in this case, it has been claimed that the theory is renormaliz- 
able [206]). The de-coupling limit is ghost free, and the longitudinal mode interacts through 
a cubic galileon interaction. In addition, the theory possesses no Boulware-Deser ghost to 
any order beyond the decoupling limit, and so it represents a completely consistent ghost 
free theory of a fully interacting massive graviton in three dimensions, one which is free of 
prior geometry. 



11.2 Topologically massive gravity 

The second type of massive gravity in three dimensions is topologically massive gravity 
(TMG) [207]. The action is 



1 / 2 



'11.21 



Here e^^'^ is the 3d epsilon tensor, equal to —^^^.'^^^ where e^^^ is the alternating symbol 
with e°^2 = /i is a mass parameter, and the term proportional to - is a Chern-Simons 
form, or secondary characteristic class. It is not a covariant tensor, but it changes into a 
total derivative under a diffeomorphism, so the action is indeed diffeomorphism invariant. 
The equations of motion are 

R^.. - Ir9>.u + -C^, = 0, (11.3) 

where C^^ = e^"^Va [R^u — \gi3uR) is the Cotton tensor. The Cotton tensor is symmetric, 
traceless, and covariantly conserved, and it vanishes if and only if the metric is conformally 
flat. 



dimensions, because the little group is so(2) in both cases. Thus a parity conserving theory of a massive 
spin s particle in three dimensions has two states, of so(2) charge ±s, which are related by parity. A parity 
violating theory, on the other hand, will carry only one degree of freedom for some particle, the right or left 
handed state. This will be the case for topologically massive gravity in Section 11.2. 
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Linearizing the equations of motion (11.3) around flat space, and plugging in a /i^J 
which is transverse and traceless gives 

(^,"5.^ + ^^.'"^a5.^) nhl^ = 0. (11.4) 

Acting on this with 5^5^ — j^e^"dx5^, we flnd 

(□ - /i^) D/iJJ = 0, (11.5) 

suggesting that the hnear curvature R^^ ~ D^^J propagates a particle with mass /i. The 
Chern-Simons term violates parity, and a more detailed analysis [207] shows that this model 
in fact propagates a single massive spin 2 degree of freedom, where the spin is left or right 
handed according to the sign of yU. This is in contrast to the parity conserving Fierz-Pauli 
mass term, which in three dimensions would propagate both a left and a right handed spin 
2. Finally, note the wrong sign Einstein-Hilbert term in (11.2), which is necessary so that 
the spin 2 particle is not a ghost. 

Adding a cosmological constant and putting this theory on AdS^ yields an interesting 
system with links to AdS/CFT. The AdS^ isometry group is 5*0(2,2) which locally decom- 
poses to SL{2,R) X SL{2,R), the group of left-moving and right-moving symmetries in a 
dual boundary 2d conformal fleld theory. This theory has ghosts or negative-energy modes 
for generic values of the parameters, but for a special choice, it was argued that half of the 
theory completely decouples, stability is restored, and the theory becomes chiral under only 
one of the SL{2, R) factors [208, 209]. For more on the subtleties, see [210, 211, 212, 213, 214]. 



12 Conclusions and future directions 

Massive gravity remains an active research area, one which may provide a viable solution to 
the cosmological constant naturalness problem. As we have seen, many interesting effects 
arise from the naive addition of a hard mass term to Einstein gravity. There is a well deflned 
effective field theory with a protected hierarchy between the cutoff and the graviton mass, 
and a screening mechanism which non-linearly hides the new degrees of freedom and restores 
continuity with GR in the massless limit. 

A massive graviton can screen a large cosmological constant, and a stable theory of 
massive gravity with a small protected mass offers a solution to the problem of quantum 
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corrections to the cosmological constant. It is a remarkable fact that the A3 theories of 
Section 9 exist and are ghost free, and that they are found simply by tuning some coefficients 
in the generic graviton potential. 

There are many interesting outstanding issues. One is the nature of the A3 theory. Is 
there a symmetry or a topological construction that explains the tunings of the coefficients 
necessary to achieve the A3 cutoff? Is there some construction free of prior geometry that 
would contain this theory? Is there an extra dimensional construction? 

There are also many questions related to the quantum properties of these theories. 
Apart from the order of magnitude estimates presented in this review and a few sporadic 
calculations, the detailed quantum properties of this theory and others remain relatively 
unexplored. The same goes for non-perturbative quantum properties, such as how a massive 
graviton would modify black hole thermodynamics, Hawking radiation or holography [215, 
216, 217, 218, 219, 220, 221]. 

The cutoff A3 is still rather low, however, so at best this theory in its current perturba- 
tive expansion can only provide a partial solution to the cosmological constant naturalness 
problem. Finding more natural constructions of these theories would go a long way towards 
solving the major issue, which is that of UV completion - is it possible to find a standard UV 
completion for a massive graviton, analogous to what the Higgs mechanism provides for a 
massive vector? Or is there some incontrovertible obstruction that forces any UV completion 
to violate Lorentz invariance, locality, or some other cherished property? If so, there may 
be a non-standard UV completion, or it could be that massive gravity really is inconsistent, 
in the sense that there really is no way whatsoever to UV complete it. There has been work 
on holography of massive gravitons in AdS/CFT [217, 218, 219, 220], which would provide 
a UV completion for theories in AdS space containing massive gravity. 

Even a partial UV completion, one that raises the cutoff to Mp, would be extremely 
important, as this is all that is required to offer a solution to the cosmological constant 
naturalness problem. One possibility is that the scale A3, while indicating a breakdown in 
perturbation theory, does not signal the activation of any new degrees of freedom, so that 
the theory is already self complete up to Mp. Since there are multiple parameters in the 
theory, it is likely that there is some other expansion, such as a small m expansion, which 
reorganizes the perturbation theory into one which yields perturbative access to scales above 
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A3. If this is true, it is important that (as is the case) the A3 theory is ghost free beyond the 
decouphng hmit. 

It should also be noted that massive gravitons already exist in nature, in the form of 
tensor mesons which carry spin 2. There is a nonet of them, which at low energies can be 
described in chiral perturbation theory as a multiplet of massive gravitons [222]. Here we 
know that these states find a UV completion in QCD, where they are simply excited states 
of bound quarks. 

We have focused in this review on theories with a vacuum that preserves Lorentz 
invariance, but there is a whole new world that opens up when ones allows for Lorentz 
violation. There exist theories that explicitly break Lorentz invariance, and theories such 
as the ghost condensate [223] which have Lorentz invariance spontaneously broken [224] 
by some non-Lorentz invariant background. In the former case, a systematic study of the 
possible mass terms and their degrees of freedom, generalizing the Fierz-Pauli analysis to 
the case where the mass term preserves only rotation invariance, is performed in [225]. For 
examples of the latter case, see [226, 227]. See also [228, 229], and [230, 231, 232] for reviews. 

There is still much to be learned about massive gravitons on curved spaces and cosmolo- 
gies (see [215, 233, 234, 235, 236] in addition to the references of Section 5). For instance, 
there are generalizations of the ghost-free higher cutoff A3 theory, both for arbitrary curved 
backgrounds [237], and for bi-gravity [238]. Is there a consistent fully interacting theory of 
the partially massless theories on de Sitter space? Are there consistent theories with cosmo- 
logical backgrounds, and in particular can they non-linearly realize the screening of a large 
bare cosmological constant while maintaining consistency with solar system constraints? 

Finally, a topic worthy of a separate review is the observable signatures that would 
be characteristic of a massive graviton. What would be the signatures of a cosmological 
constant screened by a graviton mass? For some examples of various proposed signatures, 
see [239, 240, 241, 242, 243, 244]. One surprising feature is the absence of any non-trivial 
spatially flat FRW solutions [245] (see also [246]). This means that the universe on scales 
larger than ~ 1/m must be inhomogenous or anisotropic, with homogeneity and isotropy 
restored on smaller scales by the Vainshtein mechanism. 

We will end this review by quoting the tantalizing current experimental limits on the 
mass of the graviton (under some hypotheses, of course) m < 7 x 10~^^ eV [247, 248], about 
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an order of magnitude above the Hubble scale, the value needed to theoretically explain the 
cosmological constant naturalness problem. 
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A Total derivative combinations 

Define the matrix of second derivatives 

n^. = d^d,(t>. (A.i) 

At every order in 0, there is a unique (up to overall constant) contraction of Li's that reduces 
to a total derivative^^, 

£™(n) = [n], (A.2) 

£r(n) = w-m. (A.3) 

£™(n) = [n]=^-3[n][n2] + 2[n=^], (a.4) 

£™(n) = [n]^-6[n2][n]2 + 8[n3][n] + 3[n2]2-6[n% (a.s) 



where the brackets are traces. is just the Fierz-Pauli term, and the others can be 

thought of as higher order generalizations of it. They are characteristic polynomials, terms 
in the expansion of the determinant in powers of H , 

det(i + n) = 1 + £™(n) + i£™(n) + l/:™(n) + \cT{ti) + ■ ■ ■ (a.6) 

The term £™(n) vanishes identically when n > D, with D the spacetime dimension, so 
there are only D non-trivial such combinations, those with n = 1, ■ ■ ■ , D. 



^^The proof of this fact is the same as the proof showing the uniqueness of the gahleons in [136]. See also 
[101]. 
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They can be written explicitly as 

£™(n) = J2 (-1)' v'^'^^'V'''-'''^ ■ ■ ■ v"-'^"-^ (n^...n^,,, • • • n^„,j . (a.t) 
p 

The sum is over all permutations of the u indices, with (—1)^ the sign of the permutation. 
They satisfy a recursion relation, 

cnn) = - [n-] £™„(n), (a.s) 

m=l 

with £j°(n) = 1. 

In addition, there are tensors x]j}) that we construct out of 11^;^ as foUows^^, 
The first few are 

= r^,u (A. 10) 

= mv,u-Ti,, (A.ll) 
= [[Tif-[Ti^])ri,u-2[U]Ii,, + 2Iil, (A.12) 

= ([nf-3[n][n2]+2[n=^])r^^,-3([n]^-[n2])n,, + 6[n]nj,-6nj, 
; (A.13) 



The following is an explicit expression 

{n — m)\ 



xt^ = E(-i)"7;73^n-£™^(n). (a.m) 



m=0 

They satisfy the recursion relation 



= -nn;xt~'^ + U'^'X^-^^. (A. 15) 



a/3 ^Mi'- 



Since £™(n) vanishes for n > D, X}Z vanishes for n > D. 
The xju^ satisfy the following important properties: 



^^Note that our definition of tlie X^"-* used here differs by a factor of 2 from that of [102] 
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They are symmetric and identically conserved, and are the only combinations of 11^ 
at each order with these properties: 



Qf^X^""^ = 

HP 



(A.16) 



• For spatial indices j and time index 0, 




(A. 17) 



Finally, we have the following relations involving the massless kinetic operator (2.46), 
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